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ABSTRACT 


The  Optical  Reduction  Programs  have  "been  written  to  obtain 
the  space  trajectory  of  a  re-entering  body  from  measurements  of 
position  on  two  photographic  plates  taken  of  a  re-entry  event 
from  two  separate  stations.  The  method  of  analysts  follows 
closely  that  developed  by  Whipple  and  Jacchta^,  Die  Programs  are 
a  modification  for  the  Lincoln  Laboratory  709^  computer  of  a 
program  original 3 y  written  at  the  Harvard  Observatory  for 
reduction  of  meteor  trails.  Die  computation  is  divided  into  two 
programs: 

1.  Die  Plate  Calibration  Program  calibrates  the  plate  by 
using  coordinates  and  plate  measurements  of  known  stars  as 
reference  points. 

2.  Die  Optical  Trajectory  Program  computes  range,  height, 
time,  direction  cosines,  and  distance  along  the  trail  from  the 
measured  points  on  the  trail. 

Diis  paper  discusses  these  two  programs  in  detail,  including 
program  listings,  flow  charts,  and  directions  for  running  the 
programs.  Die  mathematical  background  and  the  experimental  method 
for  obtaining  the  input  data  are  also  discussed. 
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GLOSSARY 


Sky  and  Plate  Coordinates 


Symbol 

Definition 

a,  A 

Altitude  and  azimuth  as  referred  to  the 
local  horizontal 

<5,  H 

Declination  and  hour  angle  as  referred  to 
the  equator 

6^  a 

Declination  and  ri^t  ascension  as  referred 
to  the  equator 

P,  ^ 

Celestial  latitude  and  longitude  as  referred 
to  the  ecliptic 

i,  ^ 

Direction  cosines  in  the  altitude-azimuth 
system 

1,  m,  n 

Direction  cosines  in  the  declination-hour 
angle  system 

\ 

Direction  cosines  in  the  declination-right 
ascension  system 

T 

Vernal  equinox  or  first  point  of  Aries 

€ 

Obliquity  of  the  ecliptic  (25  l/2*) 

1,  "n 

Standard  rectangular  plate  coordinates  with 
origin  at  the  optical  plate  center 

X,  Y 

Meas\ired  rectilinear  coordinates  on  the 
photographic  plate 

ff 

Angular  distance  fnom  the  optical  plate  center 

AX,  AY 

Corrections  to  measvued  rectilinear  coordinates 
on  the  photographic  plate 

Precession  direction  cosines  for  the  declination 
right  ascension  system 

vii 

Proper  motion  in  right  ascension  and  declination 
per  year 

a,b,c,d 

Plate  constants  in  four  constant  solution 

A,B,C,D 

Inverse  plate  constants  in  four  constant  solution 

a  ta  9 Si  jb  vb  tb 
x'  y  7/  rf  y  z 

Plate  constants  in  six  constant  solution  * 

Inverse  plate  constants  in  the  six  constant 
solution 

Time  and  Terrestrial  Coordinates 


S.T.  or  0 

Sidereal  time  at  observer 

Sidereal  time  at  Greenwich 

0 

o 

Sidereal  time  at  0000  Greenwich  time  at  Greenwich 

A0 

Correction  to  sidereal  time  reqviired  since  a 
solar  day  is  longer  than  a  sidereal  day 

U.T. 

Universal  or  Greenwich  time 

L 

Terrestrial  longitude 

Terrestrial  latitude  (geographic) 

Terrestrial  latitude  (geocentric) 

0 

Radius  of  earth  at  sea  level  at  observer' s 
latitude 

R 

Radius  of  earth  at  observer 

^AB'  ^AB 

Direction  cosines  in  equatorial  system  of  the 
vector  from  station  A  to  station  B 

^AB„^^AB^'^AB^ 

0  0  0 

Direction  cosines  in  the  equatorial  system  of 
the  vector  from  station  A  to  station  B  at 
sidereal  time  equal  to  zero  at  station  A 

X  V 

Zr  TtT  z 

Direction  cosines  in  the  equatorial  system 
of  the  zenithal  point  of  the  observer 

viii 


other  Symbols 


f  Lens  focal  length 

1  Length  of  the  trail  in  angular  measure 

Q  Angle  betveen  the  two  meteor  trails  seen  from 

the  two  stations 


Distance  from  station  A  to  the  plane  determined 
by  the  re-entry  object  and  station  B 

R. .  Range  from  station  A  to  the  point  i  on  the 

re-entry  object 

H  Hei^t  above  mean  sea  level 

D  Distance  along  the  trail  from  an  arbitrary  zero 

point 

P  Period  of  revolution  of  shutter 


N  Number  of  occultations  of  shutter  per  revolution 

t  Relative  time  from  arbitrary  zero  time 

^  Angle  that  blade  edge  makes  with  dash 


Subscripts 


s  Coordinates  of  stars 

c  Coordinates  of  plate  center 

q  Center  of  rotating  shutter 

Doiible  letter  subscripts  are  used  for  i>oints  eLLong  the  trail.  The  first 
letter  is  always  A  (for  station  A)  or  B  (for  station  B) 

i  A  point  on  the  trail  -  any  shutter  break  or  segnent 

b  Beginning  of  photographic  trail 

e  End  of  photographic  trail 


p 

G 

I, 

X,  Y 
R 

Plate  Calibration 
Computer  Symbol 

TITIE 

NN 

N 

xx,xy,xz 

1X,YY,YZ 

ZX,ZY,ZZ 

ACHR,ACMIN,ACSEC 

DCDEG,DCMIN,DCSEC 

XC,YC 

M 

DTK 

AC 

DC 


Pole  of  trail  considered  as  a  part  of  great  circle 
Greenwich 

Refer  to  J  and  standard  plate  axes 

Refer  to  X  and  Y  axes  on  plate 

Radiant  of  the  event 
Program 

Math  Symbol  Definition 


Title  of  event  -  any  ^2 
characters  desired 

Nimber  of  sets  of  precession 
constants 

n  Nmber  of  stars 


X  ,x  ,x 

x^  y'  z 


a 

c 


a 

c 


x  ,Y 
c^  c 


Ae 


a 

c 


a 

c 


Precession  constants 


Ri^t  ascension  of  plate 
center  in  hours,  minutes, 
seconds 

Declination  of  plate  center  in 
degrees,  minutes,  seconds 

Measured  coordinates  of  plate 
center 

Number  of  precession  set  to 
use  with  plate  center  or  star 

Correction  in  sidereal  time  to 
be  used  with  non- tracking  cameras 

Right  ascension  of  plate  center 
in  radians 

Declination  of  plate  center  in 
radians 


X 


ELC,EMC,ENC 

c 

Direction  cosines  of  unprocessed 
plate  center  in  declination 
right  ascension  system 

ELCP,EMCP,ENCP 

^  f|-L 

Direction  cosines  of  processed 
plate  center  in  declination 
right  ascension  system 

ELETA ,  EMETA,  ENETA 

^  T] 

Direction  cosines  of  Tj  axis 

ELEXI,EMEXI 

Direction  cosines  of  f  axis 

AHR,AMIN,ASEC 

a 

s 

Right  ascension  of  star  in  hours, 
minutes,  seconds 

A 

a 

s 

Ri^t  ascension  of  star  in 
radians 

DDEG,DMIN,DSEC 

6 

s 

Declination  of  star  in 
degrees,  minutes,  seconds 

D 

6 

s 

Declination  of  star  in  radians 

X,  y 

X,  Y 

Measured  coordinates  of  star 

ELS,EMS,ENS 

8 

Unprocessed  direction  cosines  of 
star  in  declination-ri^t 
ascension  system 

EL,EM,EN 

8 

Processed  direction  cosines  of 
star  in  declination-right 
ascension  system 

COSSIG 

a 

Angular  distance  of  star  from 

8 

plate  center 

EXI,ETA 

Standard  plate  coordinates 
of  star 

A,B,C,D 

a,b,c,(i 

Plate  constants  from  four 
constant  solution 

ACAP,  BCAP,  CCAP,  DCAP,  ECAP 

A,B,C,D,-A 

Inverse  plate  constants  from 
four  constant  solution 

DX,Dy 

AX,  AY 

Residuals  in  X  and  Y 

CF,SF 

cos  A0,  sin  A0 

Cosine  and  sine  of  A0 

xi 


ELC,EMC,ENC 

X  »  M  »  1/  f 

c  '  ^  c  Corrected  (for  A0)  direction 

EIEX,EMEX,ENEX 

cosines  of  the  plate  center 

1  ^  1  and  the  standard  axes  in  the 

EIET,EMET,ENET 

^  ,  ,  p  t  declination -ri^t  ascension 

T|  '  '  T)  system 

SIXCON 


SEXI2 

2_  2 
s 

SEXIET 

2 _ 

i  V 

SEXI 

2  _ 

1 

SEXEXI 

^xl 

SETA2 

2-^ 

V 

SETA 

2  _ 

SEXETA 

2  _ 

X  r\ 

SEX 

2 

X 

SEXEXI 

2  _ 
y  1 

SEXETA 

2  _ 
y 

SEX 

2 

y 

A 

Coefficient  matrix 

B 

ASEXI,BSEXI,CSEXI 

^  ”  Inverse  plate  constants 

ASETA,  BSETA,  CSETA 

,  for  six  constant  solution 

a  .0  .c 

Tj'  Tl^  T1 

xii 


Optical  Trajectory  Program 


TITIB 

—  - 

TITLE  1  (1) 

— 

TITIE  1  (2) 

— 

Constants 

CON 

p/n 

DTRAD 

V18O 

F2 

PI 

TT 

Title  of  event  -  any  72 
characters  desired 

Title  of  station  A 

Title  of  station  B 

Period  of  revolution/nvjnher 
of  occxaltations  per  revolution 

Conversion  of  degrees  to 
radians 

Pocal  length  of  camera  squared 

5.141592654 


Trail  Equation 


EN,N 

n 

SLOPE 

m 

SUMX 

2X 

SUMX2 

SUMXY 

SUMY 

X0 

^0 

YJZf 

Y 

0 

Nvmber  of  points  used  to 
get  trail  equation 

Slope  of  line  foxind  by  trail 
equation 

Sum  of  X  readings 

Sm  of  X  restdings  squared 

Svtn  of  the  product  of  the  X 
and  y  readings 

Sum  of  the  Y  readings 

First  X  reading  on  trail 

Average  Y  reading  on  trail 


xiii 


XTR 

AX 

Difference  between  actvial 
reading  of  point  and  X^, 

YTR 

ay 

Y^  for  points  on  trail” 

XBAR 

X 

Average  value  of  X  reading 

YEAR 

Y 

Average  value  of  Y  reading 

ZTRAIL 

— 

Y  intercept  of  trail  equation 

Relative  Coordinates 

EIEV 

— 

Elevation  of  station  above 
me&n  sea  level 

ELORH 

— 

Geographic  longitude  of 
station  in  hours,  minutes. 

ELONM 

— 

and  seconds  respectively 

ELONS 

— 

ELON 

L 

Geographic  longitude  of 
station  in  hours 

EIZ,EMZ,ENZ 

A 

z 

Direction  cosines  of  zenith 
of  station 

EXIABO,  ETAABO,  ZETAAB 

^AB'^^AB'^AB 

Direction  cosines  of  vector 

o  o  o 

from  station  A  to  station  B 
at  sidereal  time  zero  at 
station  A -and  at  the  time  of 
the  event,  respectively 

EXIAB,  ETAAB,  ZETAAB 

^AB  '^AB  '^AB 

Direction  cosines  of  vector 
from  station  A  to  station  B 
at  time  of  event 

GTHR,GD4,GTS 

— 

Given  event  Greenwich  time 
in  hours,  minutes,  and  seconds 

PHD,PHM,PHS 

Geographic  latitude  of 
station  in  degrees,  minutes, 
and  seconds 

m 

Geographic  latitude  of 
station  in  radians 

xlv 


CPHI,SPHI 

cos  sin  <j> 

Cosine  and  sine  of  latitude 

PHIP,PHIFD 

Geocentric  latitude  of 
station  in  radians  and 
degrees  respectively 

RAB 

^AB 

Range  from  station  A  to 
station  B 

RHO 

P 

Earth's  radius  at  the 
station 

THETJ^H,  THEO^,  THET0S 

e 

o 

Sidereal  time  at  0000  U.T. 
at  Greenwich  for  the  date 
of  the  event  in  hours, 
minutes,  and  seconds 

DTHETS 

Ae 

Correction  in  sidereal 
time  (seconds) 

THET, THETA 

e 

Sidereal  time  in  degrees 
or  radians 

CTH,STH 

cos  Q,  sin  0 

Cosine  and  sine  of  theta 

UT 

— 

Universal  time  of  event 

DL 

AL 

Difference  in  longitudes  of 
the  two  stations 

CDL 

cos  AL 

Cosine  and  sine  of  AL 

SDL 

sin  AL 

Calculation  of  Pole 

and  Radiant 

CAX, CBX,CCX 

CAY, CBY,CCY 

a  .L  ,c 

x'  x'  X 

a  ,c 

y'  y'  y 

Inverse  plate  constants  for 
six  constant  solution  (or 
four  constant  solution 
expanded  to  give  six  numbers) 

EXIBAR,ETABAR 

ItV 

Standard  plate  coordinates 
of  a  i)oint 

X 

X 

Measured  reading  of  X  on  plate 

Y 

Y 

Canptited  Y  of  point  reeid  at  X 

XV 


EIC,EMC,E1IC 


EIETA,EMETA,EIIETA 

EI£XI,EMEXI 

EL,EM,EN 

ELP,EMP,ENP 

DLONG 

1 

SINL 

sin  1 

smQ 

sin  Q 

EIR,EMR,ENR 

Direction  cosines  of  plate 
center  declination-rig^t 
ascension  system 

Direction  cosines  of  ^  axis 
declination-ri^t  ascension 
system 

Direction  cosines  of  X  axis 
declination-ri^t  ascension 
system 

Direction  cosines  of  any 
point  on  trail  declination- 
ri^t  ascension  system 

Direction  cosines  of  pole 
declination-ri^t  ascension 
system 

Trail  length  in  degrees 

Sine  of  trail  length 

Sine  of  angle  of  intersection 
between  the  two  trails  as 
seen  in  the  sky  from  stations 
A  and  B 

Direction  cosines  of  the 
radiant 


Calculation  of  Ranges,  Heights,  and  Distances 


XP 

YP 

COSZ, COSZl 


X  reading  of  any  point 

Y  Corresponding  Y  coc5)uted 

^  for  this  X  point 

cos  Z  Cosine  of  zenith  distance 

of  any  point 


D 


D 


Distance  along  the  trail 
from  an  arbitrary  zero 


EEL,iEM,PEN 


Direction  cosines  of  each 
point  in  declination-ri^t 
ascension  system 
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HEIGHT,  H 


H 
6] 

HI  h 

RANGE,  R  R 

NUM 

NDASH,ENDASH 

NWT 

EXIP,ETAP,ZETAP 


Height  above  mean  sea  level 
for  the  point 

Hei^t  correction 

Hei^t  above  the  tangent 
plane 

Range  from  the  station  to 
the  point 

Number  of  dashes 

Number  of  the  dash 
(numbering  from  1-NUM) 

Wei^t  or  reliability  of 
dash 

Direction  ccanponents  of 
point  declination-right 
ascension  system 


Calculation  of  relative  time 
XQ,IQ  X^,y^ 

XC,YC  X^,y^ 

DX,Dy 

ymPiS  X 

DITR 


CHEGA 


(O 


OCCUM?  N 

P  P 


Observed  coordinates  of  the 
center  of  rotation  of  the 
shutter 

Observed  coordinates  of  the 
projection  center 

Differences  in  two  centers 

X  reading  of  any  point 

Difference  between  Y 
computed  and  y  of  rotation 
center 

Angle  that  blade  makes  with 
dash 

Number  of  occultations  per 
revolution 

Period  of  revolution 
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SIGH,  IP  PLUS 


Code  describing  dii*ection 
of  i*evolution  of  shutter 


IJEXT 


Code  if  cards  are 
desired 


EIAZ,EMAZ,EHAZ 


Direction  cosines  of  point 
in  the  altitude -azimuth 
system 
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I.  INTRODUCTION 


The  optical  reduction  programs  for  the  709^  computer  have  been  written 
to  compute  the  space  trajectory  of  the  re-entering  bodies  of  the  Trail- 
blazer  experiments  fired  from  Wallops  Island,  Virginia.  When  a  body 
re-enters  the  atmosphere  it  is  heated  and  may  become  luminous;  thus,  it  may 
be  photographed.  Trails  are  recorded  photographically  fran  seve3ral  stations 
by  Lincoln  Laboratory  and  the  NASA  and  the  photographs  are  available  for 
analysis.  The  programs  discussed  in  this  report  compute  positions  along 
the  re-entry  trail  using  known  stars  as  calibrations.  Relative  time  may  be 
introduced  by  the  use  of  chopping  shutters  on  one  or  more  cameras,  so  that 
breaks  at  regular  time  intervals  are  produced  along  the  trail. 

The  mathematical  background  of  the  reduction  procedure  is  discussed 

in  Section  H  and  the  e3q)erlmental  techniques  in  Section  III.  The  programs 

themselves  are  documented  in  Sections  IV  and  V.  The  astrometric  reduction 

1  2 

is  modeled  after  the  technique  discussed  by  Whipple  and  Jacchla  and  McCrosky  . 
The  programs  are  adapted  from  programs  written  by  Dr.  McCrosky’ s  groxip  at  the 
Harvard  Observatory  for  meteor  analysis. 

A  camera  usually  maps  the  celestial  sphere  onto  a  flat  focal  plane. 

In  the  case  of  the  Suq>er-Schmidt  Cameras,  the  focal  plane  is  curved  and 
the  film  must  be  copied  by  projection  onto  a  flat  plate  for  analysis.  All 
measurements  for  this  reduction  are  made  on  the  flat  plates.  Positions  on 
the  plate  are  measured  as  discussed  in  Section  III. 

The  angular  direction  to  the  reference  stars  and  the  vehicle  trail 
are  detemined  with  respect  to  each  camera  station,  and  from  these  directions 
the  angular  bearing  of  measured  points  along  the  trail  can  be  found.  The 
trails  from  the  two  stations  may  look  quite  different  since  they  are 
photograjhed  at  different  distances  and  angles.  Also,  clouds  scmetlmes 
obscure  parts  of  the  trail.  Thus  the  beginning  and  end  points  of  the 
trails  recorded  from  the  two  stations  may  not  be  identical.  No  attempt 
has  been  made  to  synchronize  the  chopping  shutters  in  the  Wallops  Island 
camera  systems. 
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since  the  Trailblazer  vehicles  are  axially  symmetric,  the  only  force 
acting  to  change  the  direction  of  the  re-entry  from  a  straight  line  Is  the 
gravity  effect  on  a  freely  falling  body.  Dils  effect  Is  generally  small 
over  the  short  time  of  the  event  and  the  deviation  from  a  stral^t  line 
trajectory.  If  any,  vould  be  found  near  the  end  of  the  trail.  Also  the 
re-entry  body  of  the  Trailblazer  system  Is  fired  dovnward  and  re-enters 
about  10-12*  from  the  vertical  so  that  gravity  acts  mainly  along  the 
vehicle  trajectory.  Thus  the  re-entry  trail  Is  usviaHy  a  straight  line 
on  the  plate.  If  the  trail  Is  curved.  It  must  be  broken  Into  straight 
line  segments  and  the  program  must  be  run  separately  for  each  straight  line 
p€trt. 

Prom  the  known  positions  of  stars  near  the  trail  at  the  time  of  the 
event,  a  calibration  Is  performed  on  each  plate;  from  the  resvilts  of  this 
calibration,  the  euigular  position  of  points  along  the  trail  may  be 
determined  for  each  plate.  Since  the  trail  appears  as  a  straight  line,  the 
position  vectors  of  points  along  the  trail  as  seen  from  the  camera  station 
must  determine  a  unique  plane.  The  re-entry  must  lie  on  the  Intersection  of 
the  two  planes  determined  from  a  pair  of  photographs  taken  at  separate 
stations.  The  distance  from  the  camera  site  to  the  point  on  the  line  of 
Intersection  may  be  computed  Individually  for  each  point.  Therefore,  a 
set  of  slant  ranges  for  points  along  the  trail  may  be  found  for  each 
plate.  Height  ,  direction  cosines,  and  relative  distance  along  the  trail  may 
then  be  computed  for  each  measured  point  along  the  trail.  If  the  plates 
were  obtained  from  a  camera  with  a  chopping  shutter,  relative  time  may  be 
computed  for  each  dash. 

In  general  a  unique  time  for  any  given  point  on  the  plate  can  not  be 
found.  Also,  no  exact  correlation  exists  between  points  on  the  two  trails. 

If  a  common  point  such  as  a  flax^  can  be  Identified  on  both  plates,  the 
points  may  be  considered  to  occxir  at  the  seme  time.  The  height  of  the 
point  computed  from  both  statlcHis  should  be  Identical  (this  will  serve 
as  a  check  on  the  height  computations).  In  this  case,  also,  the  exact 


2 


sidereal  time  of  the  flsure  can  be  computed;  hovever,  since  in  general  flares 
need  not  occur,  this  computation  of  sidereal  time  is  not  included  in  the 
present  progrcm.  In  the  Project  PRESS  coded  shutters  have  been  built 
into  the  cameras  and  the  chopping  shutters  have  been  synchronized  so  that 
absolute  time  can  be  recorded  and  an  analysis  performed  using  unique  time 
points  from  two  or  more  stations.  The  Super- Schmidt  Cameras  used  in  the 
Trailblazer  pixjgram  are  not  equipped  to  record  absolute  time  since  suitable 
coded  shutters  can  not  be  easily  Installed.  It  should  be  noted  that  the 
re-entry  events  of  Project  HffiSS  are  much  longer  than  those  of  the  Trail- 
blazer  series;  thus  the  design  and  use  of  absolute  timing  equipnent  presents 
much  less  difficulty  in  Project  HIESS. 

The  present  set  of  programs  is  designed  to  compute  the  position  of 
the  body,  i.e.  hel^t,  range,  and  direction  cosines  of  each  point  along 
the  trail.  If  the  trail  is  chopped,  relative  times  can  be  computed  for 
the  dashes;  relative  distances  along  the  trail  from  an  arbitrary  point  on 
the  trail  are  also  computed.  These  programs  do  not  compute  velocity.  If 
the  distance  along  the  trail  as  a  function  of  time  is  known,  various  methods 
of  computing  velocity  might  be  used.  Harvard  has  obtained  good  results  in 
meteor  work  with  the  equation 

D  =  a  +  bt  +  c  exp  (kt) 

where  D  is  the  distance  along  the  trail;  t  is  the  elapsed  time  from  an 
arbitrary  reference  i)oint  on  the  trail  and  a,  b,  c  and  k  are  computed 
constants.  However,  this  fit  did  not  appear  particularly  suitable  for  the 
Trailblazer  work.  At  present  velocity  is  computed  from  the  resTolts  of  the 
optical  trajectory  program  by  hand  computation.  Since  this  method  is  not 
yet  sxiitable  for  machine  calcxilation,  the  computation  of  velocity  has  been 
omitted  from  the  present  program. 
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II.  MATHEMATICAL  BACKGROUND 


A.  Astronomical  Coordinate  Systems 

When  the  stars  are  observed  from  a  point  on  earth  they  appear  to  lie 
on  the  surface  of  a  sphere  with  the  observer  at  its  center.  Although  the 
eye  can  not  determine  the  distance  to  the  stars,  an  estimate  of  the  angles 
subtended  at  the  observer  by  any  pair  of  stars  can  easily  be  made.  Hiese 
directions  are  defined  in  terms  of  the  positions  on  the  surface  of  a 
sphere  -  the  celestial  sphere  -  in  vhich  the  straight  lines  Joining  the 
observer  to  the  stetrs  intersect  this  surface.  The  axis  of  the  earth  pierces 
the  sphere  at  two  points  called  "the  celestial  poles".  The  radius  of  the 
celestieil  sphere  is  considered  arbitrarily  large. 

The  position  of  any  point  on  the  surface  of  a  sphere  may  be  completely 
specified  by  reference  to  two  principal  great  circles  on  the  sphere,  or  by 
one  great  circle  and  a  point  on  that  great  circle.  Thus  a  point  on  the 
surface  of  the  earth  is  completely  specified  by  its  latitude,  measured  from 
the  eqmtor,  and  its  longitude,  measured  east  or  west  from  Greenwich. 

Several  methods  of  defining  position  on  the  celestial  sphere  are  in  common 
use. 


1.  Spherical  Coordinate  Systems 
a.  The  Altitude- Azimuth  System 

Let  0,  the  observer  on  the  surface  of  the  earth,  be  the  center  of 
the  celestial  sphere.  Z,  the  zenith,  is  the  point  directly  overhead 
(defined  by  a  plumb  line)  and  the  nadir  is  the  point  directly  underfoot. 
The  plane  throu^  0  perpendicular  to  the  direction  of  the  zenith  is  called 
the  plane  of  the  horizon  and,  as  shown  in  Ilg.  la,  cuts  the  celestial 
sphere  in  the  great  circle  NAS  dividing  the  celestial  sphere  into  two 
hemispheres,  the  upper  being  visible  and  the  lower  being  hidden  by  the 
earth  itself.  A  vertical  circle  is  a  great  circle  passing  throu^  the 
zenith  and  the  nadir. 
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Let  X  be  the  position  of  a  star  on  the  celestial  sphere  at  a  given 
moment.  Ihe  angle  AOX  or  the  great  circle  arc  AX  is  called  the  altitude,  a, 
of  the  star.  To  completely  specify  the  position  of  a  point  on  the 
celestial  sphere  the  particular  vertical  circle  on  which  it  lies  must  be 
specified.  Let  OP  be  parallel  to  the  earth' s  axis.  If  the  observer  is  in 
a  northern  latitude  the  position  P  is  called  the  north  celestial  pole.  Uie 
vertical  circle  throu^  the  pole,  ZPN  in  Pig.  la, is  defined  as  the  principal 
vertical;  it  is  the  circle  traced  on  the  celestial  sphere  by  the  observer' s 
meridian.  The  point  where  this  circle  cuts  the  horizon  is  called  the  north 
point,  N,  of  the  horizon  plane;  the  point  S  directly  opposite  N  on  the 
horizon  is  the  south  point , 

If  the  star  is  in  the  eastern* part  of  the  sphere,  as  is  X  in  Pig.  la, 
the  spherical  angle  PZX  or  the  great  circle  arc  NA  is  called  the  azimuth  A; 
if  the  star  is  in  the  western  part  of  the  sky,  as  is  Y,  the  azimuth  is  the 

-A 

angle  PZY  which  in  this  case  is  greater  than  180“.  Ihroxi^out  this  report 
azimuth  will  always  be  measured  from  north  in  an  easterly  direction  from 
0-560“.  Since  the  cuigle  POZ  in  Pig.  la  equals  the  angle  between  the  radius 
of  the  earth  which  passes  through  the  observer' s  position  and  the  earth' s 
axis,  POZ  is  seen  to  be  equal  to  the  co-latitude  of  the  observer. 

b.  The  Declination-Hour  Angle  System 

ISie  celestial  equator  is  the  great  circle  of  the  celestial  sphere 
perpendicular  to  the  polar  axis  and  is  thus  the  great  circle  in  which  the 
plane  of  the  earth' s  equator  cuts  the  celestial  sphere,  RWT  in  Pig.  lb. 

The  celestial  equator  and  the  horizon  intersect  at  two  points  E  and  W  which 
£ire  each  90*  from  N  and  S  and  eire  called  the  east  emd  west  points. 

Let  X  be  a  star  in  the  northern  hemisphere.  If  PXD  is  the  semi-great 
circle  throu^  X  and  the  pole  P,  the  sure  DX  is  defined  as  the  declination 
5  of  X,  i.e.  its  distance  in  degrees  north  (+)  or  south  (-)  of  the  celestial 
equator.  A  great  circle  passing  through  the  pole  of  the  celestial  sphere 

*If  the  observer  faces  north,  the  west  is  on  his  left  and  the  east  on  his 

rl^t. 
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is  called  a  meridian  or  an  hour  circle.  The  second  great  circle  required 

for  complete  specification  of  the  star’s  position  in  the  declination-hour 

angle  system  is  taken  as  the  observer's  meridian  (i.e.  that  meridisui 

passing  through  the  zenith)  PZRSQ  (L)  in  Fig.  lb.  At  any  moment  the 

I>osition  of  X  is  specified  by  the  angle  at  the  i)ole  between  the  observer’ s 

meridieui  and  the  meridian  EXQ  through  the  star  at  that  time.  This  angle 
/\ 

REX  or  ZEX  or  the  arc  RD  is  called  the  hoxjr  angle  H  and  is  measured  from 
the  observer' s  meridian  westward  0*  -360*  or  Oh  -24h. 

c.  The  Declination-Right  Ascension  System 

In  the  hour  angle-declination  system  only  one  coordinate,  the 
declination,  remains  constant  as  the  earth  rotates  whereas  the  hovur  angle 
increases  uniformly  from  Oh  -2Uh  during  a  day.  It  is  usually  more 
convenient  to  define  the  second  coordinate  relative  to  a  fixed  point  on 
the  equator  (althou^  all  the  stars  appear  to  change  position  during  the 
daily  rotation  of  the  earth,  their  relative  positions  with  respect  to  one 
another  remain  vinchanged over  long  periods  of  time).  The  point  chosen  as 
standard  is  the  vernal  equinox  or  the  first  point  of  Aries  denoted  by  t. 

/•S 

Then  the  arc  “T)  or  angle  Tpx.  is  called  the  ri^t  ascension  of  star  X, 
denoted  by  a,  and  is  measured  eastward  from  Oh  -24h  (or  from  0  -36O*) 

Fig.  lb.  Note:  right  ascension  is  meas\u:ed  in  the  opposite  direction 
to  hour  angle. 

d.  Ohe  Celestial  Latitude  -  Longitude  System 

Dae  earth  is  a  planet  rotating  around  the  sun  in  an  elliptical  orbit 
with  the  sun  at  one  focus  with  a  period  of  revolution  of  approximately  one 
year.  Dixring  the  year,  as  observed  from  earth,  the  sun  appears  to  make  a 
complete  circuit  of  the  sky  against  the  star  background.  The  plane  of  this 
orbit  is  called  the  plane  of  the  ecliptic  and  the  great  circle  in  which 
this  plane  cuts  the  celestial  siiiere  is  called  the  ecliptic.  This  plane 
is  found  to  be  inclined  at  an  angle  of  about  23  l/2*  to  the  celestial 
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equator;  this  inclination  is  called  the  obliquity  of  the  ecliptic,  denoted 
A 

by  €;  e  =  MTR  in  Pig.  Ic. 

Relative  to  the  earth  the  sun  appears  to  move  on  the  surface  of  the 
celestial  sphere  along  the  ecliptic  in  the  direction  Y^M,  and  twice  a  year 
its  position  on  the  celestial  sphere  coincides  with  the  intersection  of  the 
ecliptic  with  the  celestial  equator.  IDiis  position  t  at  which  the  sun’ s 
declination  changes  from  south  to  north  is  defined  as  the  vernal  equinox; 
it  is  the  reference  point  in  the  ri^t  ascension-declination  system. 

A  fourth  set  of  coordinates  can  be  derived  using  the  ecliptic  as  a 
fundamental  great  circle  and  the  vernal  equinox  t  as  a  principal  reference 
point.  In  Fig.  Ic  let  K  be  the  pole  of  the  ecliptic  and  KXA  a  great 
circle  arc  passing  throvigh  star  X  meeting  the  ecliptic  in  A.  The  arc  ta 
measiired  from  t  to  A  along  the  ecliptic  in  the  direction  of  the  sun' s 
anntial  motion,  i.e.  eastward,  is  cal  Ted  its  celestial  longitude  ^  and  is 
measured  fran  0  -360*.  The  arc  AX  is  cal  led  the  celestial  latitude  3  and 
is  measured  north  or  south  of  the  ecliptic. 

e.  Notes  on  Coordinate  Systems 

In  the  previous  definitions,  the  center  of  the  celestial  sphere  was 
taken  as  the  observer  on  the  surface  of  the  earth.  However,  since  the 
stars  are  at  distances  which  are  almost  infinitely  large  as  compared  with 
the  dimensions  of  the  earth,  no  appreciable  error  results  from  taking  the 
center  of  the  earth  as  the  center  of  the  celestial  sphere;  this  convention 
will  be  used  in  this  report  \dien  dealing  with  the  star  backgrovuid.  Although 
our  figures  have  been  drawn  for  the  noirthem  hemisphere,  similar  figures 
aijply  for  the  southern  and  the  definitions  of  the  coordinate  systems  apply 
to  both  hemispheres.  Ftor  an  observer  exactly  at  the  pole  the  definitions 

of  azimuth  and  meridian  break  down;  thus  in  the  following  study  it  will  be 
assumed  that  the  observer  is  not  at  a  pole. 

2.  Rectangular  Coordinate  Systems 

In  the  preceding  discussion,  points  on  the  celestial  si>here  were 
defined  in  tenns  of  sjiierical  coordinates;  however,  in  most  of  the  following 
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vork  it  is  much  more  convenient  to  use  a  rectangular  reference  system,  A 
ri^t  handed  coordinate  system  is  chosen  in  vhich  the  coordinates  x,  y,  z 
8u:e  measured  along  three  perpendicular  axes  such  that  the  rotation  of  the 
X  axis  into  the  y  axis  arovmd  the  z  axis  is  accomplished  hy  the  right  hand 
rule;  the  origin  of  the  system  is  taken  as  the  observer  or  the  earth' s 
center  (see  above).  Ihe  position  of  a  star  on  the  celestial  sphere  may 
be  represented  in  vector  form  by 

■where  i,  7,  and  k  are  unit  vectors  in  ■fche  x,  y,  z  directions  respectively, 
and  r^,  r^  and  r^  are  the  projections  of  r  on  "the  respective  axes.  The 
cosines  of  the  angles  between  the  positive  x,  y,  and  z  axes  and  the  vector 
r  measxired  in  the  direction  from  the  axis  toward  r  are  called  the  direction 
cosines  of  r.  We  -will  e:q)ress  all  coordinates  in  terns  of  direction 
cosines;  in  this  "way  the  magnitude  of  r  •will  not  appear  and  only  -the  directions 
■will  be  used. 


a.  Ihe  Altitude- Azimuth  System 

The  horizon  will  be  •taken  as  the  xy  plane,  •the  positive  x  axis  •will 
I)ass  ■throu^  the  east  point  of  the  horizon,  the  positive  y  suds  through 
the  north  point,  and  ■the  positive  z  axis  through  the  zeni^th. 

Let  S  in  Pig.  2a  be  any  star  on  the  celestial  sphere  •with  altitude 
a  and  azimuth  A.  !Iiien  using  V,  ^  to  represent  direction  cosines  in 
the  X,  y,  z  directions,  respectively 

6  =  cos  a  sin  A 
^  =  cos  a  cos  A 

S  =  sin  a  (2) 

b.  The  Declination- Hoxxr  Angle  System 

The  plane  of  the  equator  •will  be  taken  as  the  xy  plane,  "the  positive 
X  axis  •will  pass  -throio^  the  west  point  on  •the  equa^tor,  "the  positive  y 
axis  throu^  the  intersection  of  ■the  observer' s  meridian  and  •the  equator. 


8 


and  the  i>ositive  z  axis  throu^  the  celestieil  north  pole.  Thus,  in  the  form 
of  direction  cosines  m,  n  (Pig.  2b) 


i  =  cos  6  sin  H 

m  =  cos  <5  cos  H 

n  =  sin  <5, 


(3) 


c. 


The  Declination-Ri^t  Ascension  System 


The  plane  of  the  equator  vill  be  taken  as  the  xy  plane;  the  positive 
X  eixis  vlU  pass  through  the  vernal  equinox  t,  the  positive  y  axis  throvigh 
the  point  vhose  right  ascension  is  and  the  positive  z  axis  throu^ 

the  celestial  north  pole.  Thus  (Fig*  2c) 


k  =  cos  *5  cos  ^ 
H  =  cos  6  sin  ^ 
V  =  sin  <5 


■where  \  ^  are  direction  cosines  in  the  x,  y,  z  directions  respectively. 

B.  Precession  and  Proper  Motion 

Although  the  phenomenon  of  precession  ■was  observed  by  Hipparchus  over 
2000  years  ago,  the  first  correct  dynamical  explanation  'was  given  by  Ne'wton. 
As  is  ■well  known  the  earth  has  the  shape  of  an  oblate  spheroid  rather  than 
that  of  a  sphei^  due  to  an  equatorial  bulge  pixDduced  by  its  rotation  about 
its  polar  axis;  this  bulge  is  inclined  to  the  plane  of  the  ecliptic  by  about 
23  1/2®.  Thus  the  solar  and  limar  attractions  produce  moments  about  the 
center.  Since  the  earth  is  spinning  rapidly  on  its  axis  much  like  a  top, 
instead  of  turning  the  equatorial  plane  in-to  coincidence  ■with  the  ecliptic 
plane,  the  result  of  ■these  attractive  forces  is  to  cause  the  axis  of  the 
earth  to  process  about  the  axis  of  the  ecliptic.  Thus  the  earth' s  axis 
has  a  conical  motion  about  the  pole  of  the  ecliptic  so  that  the  pole  of 
the  equator  describes  a  small  circle  about  the  ecliptic  pole  and  the 
vernal  equinox  moves  backwards  along  the  ecliptic. 
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The  precession  of  the  equinoxes  causes  the  celestial  longitude  of  a  star 
to  increase  at  the  rate  of  about  50  seconds  of  arc  per  year  vhile  the 
celestial  latitude  shows  no  significant  change.  In  our  current  study 
we  are  interested  primarily  in  the  effect  of  precession  on  the  right 
ascension  and  declination  of  a  star.  When  the  coordinates  of  the  star  are 
expressed  in  direction  cosine  fom  precession  may  he  considered  the  result 
of  a  rotation  of  the  axes  of  the  system  with  the  origin  unchanged.  Let  us 
consider  the  effect  on  a  star  S  with  direction  cosines  ^  relative 

to  the  original  axes  X,  Y,  Z  of  a  rotation  of  these  axes  to  a  new  position 
X' ,  Y' ,  Z'  Tdiere  the  angles  between  OX' ,  OY' ,  OZ'  and  OX,  OY,  OZ  are  shown 
in  Table  I  where,  for  exan5>le,  3^  represents  the  angle  between  OX'  and  OZ. 

TABLE  I 


OX  OY  OZ 


OX’ 

X 

X 

x„ 

X 

y 

z 

OY' 

Y^ 

Y 

Y„ 

X 

y 

z 

OZ' 

z 

z 

Z 

X 

y 

z 

Thus  the  transfoimation  is  given  by  the  matrix  equation 


X  x'\ 

y  z 

Y  Y 

y  z 

\i 

z  z  , 

V  , 

"y  V  V  V 


(5) 


The  direction  cosines  ^  ,  |i',  v*  giving  the  position  of  a  star  at  any 
given  time  can  be  easily  obtained  from  the  direction  cosines  \  at  a 

standard  equinox  if  the  values  of  the  transfoimation  matrix  are  known. 
Usually  the  right  ascension  and  declination  of  a  star  are  tabulated  for  the 
eqxoinox  of  1950  in  star  catalogues;  the  values  of  the  transfoimation  matrix 

3 

can  be  found  . 
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In  1718  Halley  discovered  that  the  position  of  certain  bri^t  stars  had 

changed  appreciably  in  relation  to  the  general  star  backgroimd  since  the  time 

of  Hipparchus.  Dais  sviggested  that  these  stars  had  a  definite  space  velocity 

relative  to  the  sun  and  lay  at  finite  distances  frcm  the  earth,  Uie  proper 

motion  of  a  star  is  defined  as  its  apparent  angular  rate  of  motion  on  the 

celestial  sphere;  it  is  usually  measured  in  seconds  of  arc  per  year.  The 

proper  motion  is  usually  found  by  comparing  precise  observations  of  the 

right  ascension  and  declination  of  the  star  taken  many  years  apart.  The 

star  catalogues  give  the  components  ^  ^  the  proper  motion  in  ri^t  ascension, 

and  4  g  ,  the  proper  motion  in  declination,  per  year  or  per  century  for 

each  star  tabialated.  Since  the  components  of  proper  motion  depend  upon 

the  epoch  and  equinox  for  which  they  are  measured,  it  is  important  to 

apply  proper  motion  corrections  to  the  tabiilated  values  of  a  and  6  for  a 

star  before  precession  is  applied.  Thus,  for  exanqile,  suppose  the  ri^t 

ascension  and  declination  of  star  S  in  1935  referred  to  the  equinox  of 

1950*0  is  known  and  it  is  desired  to  conpute  the  coordinates  of  the  star 

in  I96U  referred  to  the  1964. 0  equinox.  Bie  simplest  method  is  to  first 

apply  proper  motion  corrections  frcm  the  tabulated  values  to  find  the  position 

of  S  in  1964  relative  to  the  equinox  of  1950.0  and  then  to  process  the  star 

to  the  eqinox  of  1964,  A  detailed  discussion  of  precession  and  proper 

k 

motion  may  be  fovind  in  Staart  . 

Another  effect  of  some  interest  is  nutation  which  is  a  sli^t  periodic 
oscillation  of  the  actual  pole  around  the  mean  pole  where  the  mean  pole  is 
processing  unifoimly.  Since  the  effects  of  this  nutation  are  small  no 
corrections  for  them  appear  in  this  work. 

C.  Time 

The  sidereal  time  at  any  place  on  earth  at  euiy  instant  is,  by  definition, 
the  hour  angle  of  the  vernal  equinox  at  that  instant.  When  t  is  on  the 
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observer' s  meridian  the  sidereal  time  is  0;  vhen  t  is  next  on  the  observer' s 
meridian  after  2kh  of  sidereal  time  or  one  complete  revolution  of  the  earth 
about  its  axis,  a  sidereal  day  has  passed.  Prom  Fig.  lb  it  is  obvious  that 

RT  =  RD  +  TD 

Now  RD  is  the  hour  angle  of  star  X  and  ol  is  its  right  ascension;  moreover 
RT  is  the  hour  angle  of  t.  Therefore 

Sidereal  time  (observer)  =  hour  angle  X  +  right  ascension  X 

S.T.  =  H  +  a 

Now  let  G  on  the  celestial  sphere  in  Fig.  lb  be  the  zenith  of  Greenwich; 

GPZ  is  the  longitude  of  the  observer,  GPX  the  hour  angle  of  X  from  Greenwich 
^  , 
and  ZPX  the  hour  angle  of  X  fran  the  observer' s  meridian  ^  .  Since 

GPX  =  ZK  -  =  ZK  -  longitude  of  observer 


obviously: 


K.A.  of  X  at  Greenwich  =  H.A.  of  X  at  ^  ±  longitude  of  Z 

Since  this  relation  is  general  it  also  holds  for  t.  Thus 

S.T.  at  Greenwich  =  S.T.  at  Z  ±  longitude  at  Z  (7) 

The  sign  of  Z  is  +  when  Z  is  west  of  Greenwich  and  -  when  Z  Is  east  of 
Greenwich. 

A  mean  solar  day  is  defined  as  the  Interval  between  two  successive  transits 
of  a  mean  sun*  across  the  observer' s  meridian  and  is  the  day  used  in  our 
system  of  civil  time.  A  sidereal  day  has  been  found  to  be  56.556s 
shorter  than  a  mean  solar  day. 


*The  mean  sun  is  a  fictitious  body  \diich  is  defined  to  move  in  the  celestial 
equator  at  a  uniform  rate  around  the  earth  so  that  it  completes  a  revolution 
in  the  same  time  as  the  true  sun  completes  a  circuit  of  the  ecliptic. 
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In  our  later  vork  it  vill  often  be  necessary  to  compute  the  sidereal 
time  at  a  point  on  earth  at  a  given  local  tine.  It  is  most  convenient  to 
compute  the  sidereal  time  at  Greenwich  first  and  then  to  use  equation  7 

5 

to  find  local  sidereal  time.  The  American  Ephemeris  and  Nauticeil  Alm^mc 
for  the  proper  year  contains  a  table  called  "Itaiversal  and  Sidereal  Times" 
which  gives  the  sidereal  time  for  each  day  of  the  year  for  Oh  Greenwich  time 
for  that  df^,  i.e.,  the  sidereal  time  at  midnight  at  Greenwich  for  the  given 
day;  each  day  starts  at  a  different  sidereal  time  since  the  length  of  a  day 
in  the  two  systems  is  different.  Uie  given  local  time  is  changed  into 
Universal  time  by  aidding  the  pixjper  correction  (+5  hours  for  E.S.T. ).  The 
sidereal  time  at  Greenwich  is  next  calculated  by 

Q„  =  e  +  U.T.  +  A0  (8) 

where  0^  is  the  sidereal  time  at  Greenwich,  0^  is  the  sidereal  time  at  Oh 
UT  taken  from  Ref.  5  and  is  the  correction  required  since  a  solar  day  is 
longer  than  a  sidereal  day. 

A0  =  >  56.556s  (U.T.h/24h)  =  236.556  (U.T.h/2Uh)  sec  (9) 

To  find  the  sidereal  time  at  station  A  from  its  value  at  Greenwich  at  the 
same  time  we  merely  use  equation  7  and  the  whole  equation  then  becomes 

0=0  +  UT  -  L.  +  A0  (10) 

A  o  A 

where  0^  is  the  desired  time  at  A  and  is  the  longitude  of  A  (west  in 
this  case). 

D.  Transformation  of  Coordinates 

In  our  later  work  it  will  be  necessary  to  transfom  between  the 
altitude  -  azimuth  and  the  declination  -  right  ascension  systems.  Let  us 
consider  the  problem  of  converting  from  the  ri^t  ascension-declination 
system  to  the  altitude-azimuth  system  in  direction  cosine  form.  We  are 
given  Q!  and  d  for  star  X  at  a  given  time  t  and  a  given  latitude  We  must 
first  transform  our  three  direction  cosines  \  M,  to  the  1,  m,  n  cosines 
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of  the  declination-hoxir  angle  system  (Fig.  3a). 

It  is  obvious  that  n=  ^  since  the  axis  is  identical  for  the  two  systems 
and  the  x,  y  plane  is  the  same  in  both  cases. 

0  =  a  +  ji  where  9  is  sidereal  time  at  the  observer 
so  that  the  transformation  may  be  seen  to  take  place  as  in  Fig.  3b 
where  the  rotation  of  axes  in  the  x,y  plane  is  shown. 

1  =  A  cos  (90  -  sin  (90  -  0)  =  >^  sin  0  -  u  cos  0 

m  =  X  sin  (90  -  0)  +  li  cos  (90  -  0)  =  X  cos  0  +  (i  sin  0  (ll) 

n  =  »^ 

Next,  we  must  transform  from  the  declination-hour  angle  system  to  the  altitude- 
azimuth  system  (Fig.  3c).  Ihe  y,z  -plane  is  camnon  to  both  systems  and  x^ 

=  -  Xg,  so  that  immediately 

I  =  p  CO80  -  X  sin0 

Also,  PZ  =  90  "  0  so  that  the  transformation  may  be  represented  in  the  plane 
as  in  Fig.  3d. 

5  =  M  COS0  -  sin0 

Tj  =  -  [m  cos  (90  -  0)  -  n  sin  (90  -  0)  ]  =  ^  cos  0  -  (>k  cos0  +  M  sin0)  sin  0 
=  m  sin  (90  -  J^)  +  n  cos  (90  -  JiJ)  =  (X  cos0  +  M  sin0)  cos  0  +  i/sin  0 

(12) 


E.  Astronomical  Riotography 

A  section  of  the  celestial  sphere  can  be  photographed  by  using  the 
general  configuration  of  Fig.  h.  Let  CO  be  the  opticeO.  axis  of  a  lens  with 
foceO.  length  f  and  assvnne  that  CO  passes  through  the  geometric  center  of  the 
plate.  Bie  focal  plane  of  the  lens,  FG,  is  perpendicular  to  OC.  Let  us 
consider  a  star  on  the  optical  eocis  at  C;  aJ]  the  rays  from  this  star  falling 
on  the  lens  will  be  brou^t  to  a  focus  at  0  and  an  image  formed  at  that 
point.  If  another  star  lies  on  the  celestial  sphere  at  S,  it  is  seen  to  form- 
an  image  at  R.  If  the  astronomical  coordinates  of  C  and  S  are  known,  since 
f  is  the  known  focal  length  of  the  lens  eind  o"  =  b,  we  can  find  the  relation 
between  the  linear  distance  d  on  the  photographic  plate  and  the  corresponding 


14 


angular  distance  o'  in  the  sky.  Thus 


d  =  f  tan  h  =  f  tan  o 


(15) 


We  will  now  define  a  system  of  axes  on  the  plate  through  the  plate 
center  such  that  the  position  of  a  star  on  the  plate  can  be  related  to  its 
astronomical  coordinates.  First,  the  direction  cosines  for  the  plate  center 
in  the  right  ascension-declination  system  are  found  by  applying  equation  4 
to  the  special  case  of  the  plate  center. 


=  cos 

6 

cos 

a 

c 

c 

c 

=  cos 

6 

sin 

a 

c 

c 

c 

p 

=  sin 

6 

Standard  coordinates  |  and  q  will  be  rectangular  coordinates  in  the  plane 
tangent  to  the  celestifiLL  sphere  at  C  with  their  origin  at  C  such  that  "rf 
Increases  northerly  along  the  hour  circle  through  C  and  |  increases  toward 
the  east  tangent  to  the  parallel  of  declination  of  C  (Fig.  5)* 

Having  found  the  direction  cosines  of  the  plate  center  in  terms  of  the  decli¬ 
nation  ri^t-ascension  system,  it  is  now  desired  to  determine  the  direction  cosines 
of  the  plate  axes  |  and  tj,  defined  above,  using  the  same  system.  Hie  angle 
OCP'  is  a  right  angle  since  q'  on  the  sphere  must  be  perpendicular  to  a  radius 
OC  of  the  sphere.  Therefore,  OP'C  must  equal  and  so  =*  cos  Also 

the  projection  of  tj  on  the  equatorial  plane  must  be  -  sin  so  that  the 
direction  cosines  of  tj  eire  given  by 


=  -  Sin  «c 

=  -  sin 
V  =  cos  ^ 

TJ  C 


COS  ot 
c 

sin  ^ 
c 


(15) 


Since  I  increases  along  a  tangent  to  a  parallel  of  declination,  its 
projection  on  the  z  axis  is  zero  so  that  =  0.  The  |  axis  is  directed 
eastward  from  the  point  on  the  sphere  representing  the  plate  center;  hence 
its  direction  is  +  90*  (J^g*  5)*  Hence: 
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\  =  -  Sin 

a.  =  +  cos  “ 

I  c 

=  0  (16) 

Nov  let  S  "be  a  star  on  the  celestial  sphere  with  coordinates  ^  ) 

s  s 

and  direction  cosines  in  our  normal  equatorial  system  (\  Let  us 

find  the  direction  cosines  of  S  relative  to  the  coordinate  system  of  the 
photographic  plate  (J,  "ri).  Bie  direction  cosines  of  the  standard  plate 
system*  s  axes  relative  to  the  equatorial  axes  are  shown  in  Fig.  6b.  Thus 
the  direction  cosines  of  S  in  the  plate- centered  system  are 


P|  -  Ys  *  “l  •'s 

^  ’  W  *  ^"8 

'  Vs  “o“s  *  ‘'o's 


It  is  also  true  that 


cos  ^  -  "K  7\  +  V  V 

s  c  s  *^c*^s  c  s 


(IT) 

(18) 


where  O’  is  the  angle  between  the  star  and  the  plate  center.  From  Fig.  4 

it  can  be  seen  that  the  star  located  on  the  celestial  sphere  at  S  is  seen  by 

the  camera  at  R,  as  if  it  were  located  on  the  tangent  plane  at  D.  This  type 

of  projection  is  called  a  gncmionic  projection  and  it  is  assumed  that  the 

camera  lens  is  located  at  the  center  of  the  celestial  sphere.  The  gnomonic 

projection  is  accomplished  by  dividing  equation  17  by  the  cosine  of  ^  . 

s 

One  then  obtains  the  pair  of  standard  plate  equations  for  the  star  S; 

'  (  \\  *  !  “="s 

■'is  '  (\’'s  *  *  W'>  /  '“"s  (19) 


F.  Plate  Constants 


Let  us  now  assume  that  we  have  a  photographic  plate  containing  a  re-entry 
trail  against  the  star  background  with  the  camera  driven  by  an  equatorial 
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moiint  so  that  the  stars  appear  as  points  on  the  plate.  The  trail  of  the 

re-entry  object  will  appear  as  a  broken  line  if  a  chopping  shutter  is  used. 

The  first  i)art  of  the  analysis  will  be  the  calibration  of  the  plate.  This 

will  allow  us  to  find  the  right  ascension  and  declination  of  any  point  on  the 

plate  (at  least  in  the  region  of  interest)  when  the  X  and  Y  measurements  on 

the  plate  are  known.  A  number  of  stars  (probably  about  6-15)  are  identified 

near  the  trail;  their  right  ascensions  and  declinations  are  found  from 

standard  tables  for  a  given  date.  After  corrections  etre  made  for  proper 

motion  and  precession  on  the  tabulated  coordinates,  we  have  exact  values 

of  ^  and  ^  for  each  star  at  the  time  of  the  event  *,  so  that  we  now  know 
s  s 

the  astronomical  coordinates  of  points  on  oxir  plate.  An  approximate  position, 

and  ^  for  the  plate  center  is  required  so  that  the  standard  plate 
c  c 

coordinates  and  can  be  calculated  for  each  star*  using  equation  19* 

Tvro  methods  of  relating  oi  ,  ^  of  the  star  and  the  X,Y  readings  on  the 
plate  are  in  common  use.  Hie  measvirements  are  taken  so  that  the  X-axis  is 
parallel  to  the  trail  of  the  re-entry  and  the  Y-axis  is  perpendicular  to 
the  X.  In  the  "four  constant  solution"  we  solve  the  two  equations 


af  +  bT|  +  c  =  X 

^S  's  ! 

bf  -  aTi"  +  d  =  Y 

^S  's  1 


s 


(20) 


s 


Two  stars  at  the  ends  of  the  trail  are  used  to  compute  the  four  constants 
a,b,c  and  d  and  the  standard  coordinates  of  any  point  on  the  trail  may  then 
be  found  from  the  inverse  transforms: 


(21) 


vhere: 


A  =  a  /  (a^  +  b^) 

C  =  -  (ac  +  bd)  /  (a^  +  b^) 


B  =  b  /  (a^  +  b^) 

D  =  (ad  -  be)  /  (a^  +  b^) 


^Proper  motion  is  corrected  to  the  nearest  month  and  precession  to  the 
nearest  1  January. 


IT 


The  "six  constant  solution"  requires  a  solution  of  the  equations 

a  £  +1311  +c  —  X;a  £  +c  =  Y  (22) 

X  ^s  X  's  X  s'  y  ^s  y  's  y  s  '  ' 


for  the  constants  a  ,b  ,c  ,a  ,b  ,  and  c  by  a  least  squares  method. 

3c  X  X  y  y  y 

standard  coordinates  ceui  then  be  found  by 


The 


^  *  “n  \  *  =T| 


•where: 


a.  =  yc 


I 


xy 


i 


=  - V‘ 


xy 


c,  =  (-b  c  +  b  c  )/c 

5  '  y  X  X  y''  xy 

c  =  a  b  -  a  b 
xy  X  y  y  X 


=  -  a  /c 

y  xy 

“  V"xy 

^®’y°x  “  ®’x‘^y^^°xy 


(23) 


The  six  constant  solution  is  more  complex.  However  it  has  the  ad'van'tage 
that  the  X  and  Y  measuring  axes  need  not  be  exactly  perpendicular  so  that 
it  can  be  used  to  correct  for  any  non-orthogonality  of  a  measuring  engine. 

It  will  also  correct  for  any  linear  distortions  on  the  plate. 

When  the  plate  cons'tants  have  been  found  a  process  of  second  order 
corrections  is  usually  applied.  Biese  corrections  make  allowances  for  errors 
caused  by  optical  distortions,  distortions  caused  by  -the  emulsion,  misplacement 
of  •the  plate  center,  errors  in  the  screws,  etc.  It  is  now  possible  to  compute 
theoretical  -values  of  X  and  Y  for  each  s-tar  —  i.e.,  the  position  where  a 
star  of  the  given  right  ascension  and  declination  should  lie  on  the  plate. 

From  these  -values  residuals  eire  compirted  from  -the  equations: 

^  Xg  =  Xg  (comp)  -  Xg  (obs)  ^  =  Y^  (comp)  -  Y^  (obs)  (23) 

The  observed  -values  of  X  and  Y  for  points  on  the  trail  are  usually  corrected 
by  the  equations: 


X  (corr)  =  X  (meas)  +  ^  (curve)  Y  (corr)  =  Y  (meas)  +  AY  (curve) 
where  AX  (curve)  and  AY  (curve)  are  -values  found  from  a  smooth  curve  through 
the  residuals  AX  and  AY.  13ie  measured  X  is  the  abscissa  for  the  correction 
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curves  and  the  AX  or  AY  is  the  ordinate.  By  examining  the  residueils  it  is 
usually  easy  to  spot  am  obviously  misidentified  or  misread  star  since  its 
residuals  will  be  much  larger  than  normal.  Often  one  or  more  stars  can  be 
eliminated  and  the  calibration  rerun. 

When  the  plate  calibration  has  been  completed,  it  is  possible  to  relate 
the  X,  Y  readings  of  any  point  on  the  plate  to  a  position  on  the  celestial 
sphere.  First  the  standard  plate  coordimtes  are  found  by  equation  21  or  23. 
Direction  cosines  of  the  point  can  then  be  obtained  directly  by; 

Tv  =  (\  +  +  Tv^I)/  (1  +  f  +  t) 

=  (^c  +  +  Mgl)/  (1  +  +  Ti^) 

=  (I^C  (1  +  1^+  Ti^) 

G.  T5ie  Trail  Equation 

It  is  assvnned  in  the  following  reduction  that  the  re-entering  vehicle 
travels  in  such  a  path  that  its  projection  on  the  plate  is  a  straight  line. 

If  this  is  not  the  case  the  reduction  must  be  done  separately  on  straight 
line  segments  of  the  trail.  It  Is  therefore  possible  to  ccanpute  a  straight 
line  of  the  foim 

Y  =  mX  +  b  (25) 

where  m  is  the  slope  of  the  line  and  b  the  y  intercept;  thus  only  X  need  be 
measured  directly  from  the  plate  after  the  equation  of  the  line  is  found. 
Reading  to  compute  the  trail  equation  is  done  by  dividing  the  trail  into  a 
suitable  number  of  equally  spaced  (usually  between  10  and  20)  sections  in 
X  and  reading  the  Y  corresponding  to  these  readings.  Then  the  X  and  Y 
measurements  are  related  by  a  least  squares  foimula 

m  =  (S  x^Y^  -  n  X  Y)/(2  X^^  -  n  X^)  b  =  Y  -  m  X  (26) 

where  n  is  the  number  of  points  read  and  Y  and  X  are  average  values. 

H.  Ihe  Pole  of  the  Trail 

As  discussed  above  the  vehicle  appears  to  travel  in  a  great  circle  if 
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the  trail  is  straight  along  the  plate.  If  the  trail  is  not  straight  the  pole 
is  computed  by  using  points  from  the  early  straight  line  i)art.  The  i)ole  of 
a  great  circle  is,  by  definition,  a  point  90“  from  each  point  on  the  great 
circle.  In  computing  the  pole  of  a  trail  we  will  use  the  beginning  and  end 
points  (a  middle  point  may  be  used  as  a  check).  It  is  readily  seen  that 
the  pole  of  a  great  circle  throu^  the  trail  can  be  fovind  from  the  cross 
product  of  the  beginning  and  end  points  on  the  trail.  Hence,  the  direction 
cosines  of  the  pole  given  by 


A 

sin 

f 

=  UP  .  pU 

P 

be  be 

P 

sin 

i 

11 

1 

p 

V 

sin 

i 

II 

1 

o'" 

(27) 

where  are  direction  cosines  of  the  beginning  point  of  the 

trail  and  A  ,  M  ,  those  ( 
e'  e'  e 

trail  (in  angiilar  measure). 


trail  and  A^,  those  of  the  end  point,  and  ^  is  the  length  of  the 


I.  Relative  Coordinates  of  the  Two  Stations 


It  is  necessary  to  use  two  photographic  plates  taken  from  two  separate 
stations  to  perfoim  a  complete  optical  reduction.  These  stations  will 
hereafter  be  called  station  A  and  station  B.  We  will  now  compute  the  relative 
positions  of  the  two  stations  with  respect  to  each  other.  We  will  assume  that 
we  know  their  geographic  latitudes  and  longitudes,  their  elevations  above  sea 
level  and  the  exact  time  of  the  event.  !Ihe  coordinate  system  will  be  that 
shown  in  Fig.  7*  We  must  first  find  the  geocentric  latitude*  of  the  two 
stations  since  we  eure  working  in  a  geocentric  system;  we  will  also  find  the 
earth' s  radii  and  R^  for  stations  A  and  B,  respectively,  using  the 
international  ellipsoid.  We  find^ 

^  -  11' 35* "6635  sin  2^  +  l."1731  sin  4^  -  0. "0026  sin  6^ 


* 

Geographic  latitude  is  measured  at  the  earth' s  surface  with  reference  to  the 
local  vertical;  geocentric  latitude  is  measured  from  the  center  of  the  earth. 
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P  =  a(0.9983200U7  +  O.OOI683U94  cos  20  -  0.0000035^9  cos  ^  +  O.OOOOOOOO8  cos  60) 
a  =  6378.388  Kin  =  20926.42796  Kft 

R  =  P  +  elevation  of  station  (28  contd) 


Letting  0^,  0^  be  the  geographic  latitudes,  the  earth' s  radii,  0^,  0^ 

the  geocentric  latitudes  and  the  longitudes  of  stations  A  and  B 

respectively,  ve  can  find  the  relative  rectangtOar  coordinates  in  the  right 
ascension-declination  system  of  station  B  with  respect  to  station  A  at 
sidereal  time  zero  at  station  A.  Note;  as  0^  =  0  the  x  axis  goes  through 
the  meridian  of  station  A. 


^3  MS  (1^  -  MS 

'  '*B  ^3  ('B  -  "'a’ 

‘aB  '  '*B  ^B  -  \  s*" 


and  obviously 

“ab  ^AB  AB  AB  AB 

00  o  o 

Next  we  must  compute  the  relative  coordinates  of  station  B  with  respect 

to  A  at  the  time  of  the  event  6.,  the  sidereal  time  of  the  event  at  station 

A' 

A.  In  general  0^  f  0^.  At  this  time  the  system  will  appear  as  in  Pig.  Tt>-  Thus 
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(50) 


We  also  make  use  of  the  direction  cosines  of  the  zenithal  points  for  the 
two  stations.  From  Pig.  7b 
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cos  0^ 

cos  0. 
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II 

cos  0_  cos 
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sin  0. 
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II 
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jd 
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sin  0^ 

^  = 
ZB 

sin  0g 

(31) 


where  0  is  now  the  geographic  latitude. 
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J.  The  Radiant 


The  radiant  of  a  meteor  or  other  moving  body  is  defined  as  the  direction 
from  which  it  appears  to  have  originated.  In  our  case  this  point  is  easily 
seen  to  be  the  point  at  which  the  projected  great  circles  defined  by  the 
trails  on  the  two  plates  intersect.  Die  angle  between  the  two  great  circles 
is  the  same  as  the  angle  between  their  poles;  the  radiant  is  the  cross 
product  of  their  poles.  Hence,  using  equation  27  for  stations  A  and  B 

\  sin  Q  =  ‘"gp  -  iAgp 

^  "  ^AP  ■  ^AP  ‘'bp 

‘'r  ^  ^  ^AP  ^BP  "  ^AP  ^P  (32) 

where  Q  is  the  angle  of  intersection  between  the  two  meteor  trails  as  seen 
from  station  A  and  station  B.  Die  radiant  chosen  must  lie  above  the  horizon; 
this  is  done  by  requiring  the  dot  product  Z*R  to  be  ixisitive  or  zero,  (Z 
being  a  vector  in  the  zenith  direction). 


K.  The  CoD5)utation  of  Range 


Since  the  trajectory  of  the  re-entering  vehicle  is  assumed  to  lie  on  a 
great  circle  only  one  possible  plane  exists  through  the  re-entry  track  and 
the  observing  station;  the  pole  for  that  great  circle  is  perpendicular  to 
this  plane.  The  distance  s  frcan  the  first  station  (a)  to  the  plane  determined 
by  the  missile  and  the  second  station  (B)  is  seen  to  be  (Fig.  8)  the  projection 
of  R^  (equation  29)  onto  a  perpendicular  to  the  plane,  i.e.,  onto  the  line 
through  the  pole.  Hence,  using  the  direction  cosines  of  the  pole  (equation 
27)  we  get  the  distance 


^A  "  ^AB  ^BP  ^AB  ^BP  ^  ^AB  ‘'bP 


(33) 


The  distance  or  range  R^^  of  station  A  to  any  point  on  the  plane  is  given  by 
R^^  =  sec  where  is  the  angle  between  P^  and  the  range  vector  R^^. 

Since  points  on  the  re-entering  track  are  just  specific  points  on  this 
plane  we  have  the  equation 
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cos  ■  ^Ai  \p  ^BP  A1  BP 


+  V 


SO 


(3^^) 


L.  Ihe  Computation  of  Height 


In  computing  height  above  sea  level,  the  height  h^^  above  a  plane 
throvigh  the  station  nomal  to  the  zenith  direction  is  first  computed.  The 
zenith  distance  for  any  point  on  the  re-entry  is  given  by 


cos 


“Ai  ^ZA  \±  ^  ^ZA  ^Ai  '  "ZA  '  Ai  (55) 

see  Fig.  9*  Thus  the  hei^t  above  the  plane  nonnal  to  the  zenith  is  found  by 
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(56) 


To  obtain  hei^t  above  sea  level,  the  elevation  of  the  station  H^,  and  a 
correction  6  h^^  for  the  earth' s  curvature  must  be  added.  As  in  Fig.  9 
the  total  height  of  the  vehicle  from  the  center  of  the  earth  is  given  by 
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and  thus  solving  for  6  h^^ 
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The  total  hei^t  above  mean  sea  level  is 
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where  H.  is  the  elevation  of  the  station. 
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(57) 


(38) 


M.  The  Calc\ilation  of  Distance  Along  the  Trail 


After  the  ranges  have  been  computed  for  any  two  jjoints  on  the  re-entry 
track  the  distance  between  them  c£m  easily  be  found.  Hie  ccmponents  of  the 
vector  in  the  astronomical  equatorial  system  are 
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Letting  the  first  point  seen  on  the  trail  he  called  ve  cam  find  the 
distance  from  the  beginning  of  the  tred.1  to  any  other  point  by 

“ai  =  ["ai  -  6ai)"  ("Al  -  *  (^Ai  -  (^0) 

N.  The  Computation  of  Relative  Time 


Except  in  special  cases  vhere  a  caramon  point  cam  be  identified  along 
the  trail  in  both  photographs,  it  is  impossible  to  find  exactly  simultaneous 
points  on  the  two  plates;  also,  no  absolute  time  is  available  for  any  point. 
However,  if  one  or  more  of  the  photographs  has  been  chopped,  a  relative  time 
(i.e.,  time  measured  from  am  arbitrary  zero  point)  can  be  computed  and  maikes 
possible  the  calculations  of  height  and  range  as  functions  of  time  and  also 
am  estimate  of  the  velocity  of  the  vehicle  as  the  time  derivative  of 
distance  ailong  the  trail. 

The  cameras  used  to  compute  relative  time  employ  a  shutter  which  opens 
and  closes  at  known  constamt  intervals,  thus  making  chops  on  the  plate. 

The  coordinates  of  the  mid-point  of  each  chop  aure  read;  these  points 
correspond  to  known  relative  times  since  the  chopping  rate  is  known. 

If  P  is  the  period  of  revolution  of  the  rotating  shutter  and  N  the 
nvmber  of  occultations,  all  equal  .ly  spaced,  per  revolution,  the  instant  t^ 
corresponding  to  the  i'th  break  is 


i(|) 


+  A  t. 


(41) 


where  A  t^,  is  a  function  of  the  position  of  the  re-entry  body  on  the  plate 
relative  to  the  center  of  rotation  of  the  shutter.  lOr  calculations  not 
involving  Super-Schmidt  cameras  it  will  be  assumed  that  A  t^  is  negligible. 

For  Sxiper- Schmidt  cameras  the  correction  is  often  fairly  lairge.  As  a 
first  approximation  let  us  assvnne  that  the  projection  center  (X^,  Y^)  very 
nearly  coincides  with  the  center  of  rotation  of  the  shutter  ^q^' 

practice  this  approximation  is  usually  valid  for  a  distance  between  the 
two  centers  less  than  one  cm.  In  Pig.  10  let  X^,  be  the  coordinates  of 
the  i'th  point  on  the  trail.  The  angle  that  the  blade  edge  makes  with  the 
i'th  dash  is  obviously  given  by 
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tan  Uj  =  (Xj  -  X  )  /  (Y^  -  Y  ) 


If  is  the  value  of  tc  at  the  origin  of  the  time  scale  we  have 
A  t. 


ik2) 


=  ±  (co^  -  u;^)  p/2  ^ 

where  the  +  sign  is  used  when  the  shutter  edge  chases  the  re-entry  body  as 
in  Fig.  10,  since  here  more  than  one  revolution  of  the  shutter  occurs 
between  chops  and  the  -  sign  is  used  when  the  edge  runs  toward  the  body. 

If  the  distance  between  the  center  of  rotation  of  the  shutter  and  the 
projection  center  is  too  large  for  equation  42  to  be  valid  an  exact  expression 
for  tan  has  been  found  by  Whipple  and  Jacchia^  -to  be 

u  (Xj  -  X  )  -  (Yj  -  V,)  1  /  U  di  -  \)  *  (X,  -  X  ) 

‘  ^  (W) 


tan  =  z 
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(Y  -  Y  )‘ 
'  q  c' 


and  f  is  the  focal  length  of  the  camera. 

0.  Ihe  Problem  of  Stationary  Cameras 


One  or  both  of  the  photographs  used  in  the  analysis  may  have  been  made 
by  a  stationary  camera.  In  this  case  the  star  image  consists  of  a  trail 
rather  than  a  point  on  the  plate.  Bie  ends  of  this  trail  corresjond  to  the 
opening. and  closing  of  the  camera.  Ihus  the  end  of  the  star  trail  where  the 
measurement  is  made  is  not  simultaneous  with  the  time  of  the  event.  The 
declinations  of  the  stars  are  unchanged  but  the  right  ascensions  are  shifted 


by  an  amount  A0  =  0, 


01  where  0.  is  the  sidereal  time  of  the  event  and 
A  A 


01  the  time  at  the  end  of  the  star  trail  where  the  plate  is  measured. 

Therefore  the  direction  cosines  ^Ai'*Al^  computed  by  equation  4 

must  be  changed  to  the  true  direction  cosines  by 
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It  is  easily  shown  that  if  only  the  direction  cosines  of  the  plate  center 
and  of  the  f  and  "ri  axes  are  corrected  by  equation  45  from  the  time  of  the 
star  trail' s  end  to  the  time  of  the  event,  the  direction  cosines  of 
each  point  on  the  re-entry  trail  computed  from  equation  2k  will  be  correct. 

HI.  THE  EXPERIMENTAL  OPTICAL  HIOGRAM 

Bie  optical  reduction  method  outlined  in  the  present  paper  for  detemnining 
the  space  trajectory  of  a  re-entering  body  requires  that  the  event  be 
photographed  simultaneously  from  two  stations  separated  by  a  relatively  long 
baseline.  If  velocity  is  to  be  computed,  at  least  one  of  the  cameras  must 
be  equipped  with  a  chopping  shutter  so  that  relative  time  may  be  obtained. 

The  Trailblazer  program  utilizes  two  Baker  Super-Schmidt  cameras,  one  at 
Arbuckle  Neck,  the  other  at  Eastville.  In  addition,  photographs  taken  by 
NASA  on  ballistic  cameras  from  their  Wallops  Island  and  Coquina  Beach 
stations  are  usually  available.  Figure  11  shows  the  relation  of  these  camera 
sites  to  each  other  and  to  the  launch  site. 

A.  The  Super-Schmidt  Camera 

Figure  12a  is  a  photograph  of  a  Baker  Super-Schmidt  camera  used  to 
obtain  the  re-entry  photographs.  Ihis  camera  has  an  aperture  of  12.25  inches, 
a  focal  length  of  8  inches,  an  effective  focal  ratio  of  O.85  and  a  field  of 
view  of  approximately  58  degrees.  It  is  supplied  with  an  interchangeable 
shutter  and  motor  capability,  so  that,  depending  upon  the  number  of  blades 
in  the  shutter  and  the  rpm  of  the  motor,  different  chopping  speeds  are 
possible.  The  surface  of  the  shutter  is  Immediately  in  front  of  the 
photographic  emulsion. 

A  schematic  diagram  of  the  optical  system  appears  in  Fig.  12b.  Two  part- 
hemispherical  correcting  shells,  concentric  with  the  center  of  curvature  of 
the  spherical  mirror,  provide  almost  complete  correction  of  spherical 
aberration.  Any  remaining  aberrations  and  achromatism  over  the  range  38OO- 

O 

7000  A  are  corrected  for  by  an  aspherical  correcting  plate  through  the  center 
of  curvatxire.  The  spherical  focal  surface,  concentric  with  the  rest  of  the 
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system,  lies  just  within  the  second  shell.  A  light  ray  entering  the  system 
passes  first  through  the  outer  shell,  the  correcting  plate,  and  the  second 
shell  before  it  is  reflected  at  the  mirror;  it  then  passes  throu^  the 
inner  shell  a  second  time  before  it  reaches  the  focal  surface. 

Since  the  focal  sxirface  of  the  Super-Schmidt  is  spherical,  flat  film 
must  be  molded  to  this  shape.  Ihe  film  is  held  in  the  focal  plane,  convex 
toward  the  mirror,  by  a  partial  vacuvnn  (about  20"  Hg)  between  the  film  and 
the  focal  plane.  A  special,  very  fast,  blue  sensitive  X-ray  film  is  usually 
employed  in  these  tests.  It  is  made  by  Eastman  Kodak  especially  for  the 
Harvard  Meteor  Project. 

Hie  optical  projection  in  the  Super- Schmidt  camera  is  very  nearly  a 
truly  concentric  sphericeLL  projection;  this  allows  the  star  images  to 
maintain  their  relative  positions  on  the  celestial  sphere  without  distortion. 
The  curved  film  is  projected  onto  a  flat  plate  by  a  very  nearly  gnomonic 
projection;  this  flat  plate  (a  positive)  is  then  used  for  the  analysis. 

Although  much  of  our  work  is  done  with  the  Staler- Schmidt  cameras  the 
reduction  methods  and  the  computation  program  discussed  here  will  work  with 
any  set  of  ballistic  plates. 

B.  Star  Identification 

At  each  station  the  re-entry  trail  and  the  field  stars  are  photographed 
simultaneously.  If  the  photograph  is  taken  by  a  camera  equipped  with  an 
equatorial  mount,  such  as  the  Super- Schmidt,  the  stars  appear  as  points  on 
the  plate  (Fig.  13a) i  if  the  camera  is  stationary,  the  stars  will  appear  as 
streaks  (Pig.  13b).  However,  if  the  time  of  opening  and  closing  of  the 
shutter  is  known  the  ends  of  these  streaks  are  identifiable  in  time  and  may 
be  used  in  the  analysis  as  discussed  in  Section  II-O. 

The  approximate  coordinates  of  the  plate  center  in  right  ascension  and 
declination  are  recorded  by  the  camera  operator;  the  direction  of  north  on 
the  photograph  is  also  indicated.  Thus  the  general  position  and  orientation 
of  the  photograph  are  known.  Ibr  Sviper- Schmidt  plates,  in  particulsir,  the 
use  of  Norton’ s  Star  Atlas  for  a  preliminary  identification  of  the  bri^ter 
stars  is  recommended  since  the  scale  of  its  star  maps  is  very  nearly  that  of 
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the  Super- Schmidt  film.  ^  overlaying  the  film  on  the  charts,  many  of  the 
stars  may  be  immediately  identified  and  the  same  stars  easily  found  on  the 
BD  or  CD  chart,  thus  outlining  the  area  of  interest.  Die  position  of  the 
trail  among  the  stars  is  then  found  on  the  BD  or  CD  charts  and  may  be 
sketched  in  lightly  for  future  reference. 

Stars  to  be  measured  should  be  chosen  at  uniform  intervals  along  the 
trail,  (X-direction)  and  alternately  on  opposite  sides  of  the  trail  (+  and 
-  Y).  In  order  to  minimize  centering  errors  the  stars  chosen  should  be  of 
neaorly  the  same  photographic  magnitudes  and  as  near  the  plate  limit  as 
possible.  Die  stars  to  be  measured  may  be  marked  on  the  star  charts  and 
their  BD  or  CD  numbers  as  well  as  their  approximate  coordinates  recorded 
from  the  BD  or  CD  tables.  When  the  BD  or  CD  number  of  a  star  is  known 
its  position  may  be  found  in  the  more  exact  star  catalogues. *  The  following 
information  should  be  recorded  for  each  star: 

Right  Ascension 

Declination 

Equinox  of  Reference 

Epoch  of  Observation 

Components  of  Proper  Motion  in  Right  Ascension  and  Declination. 

Die  stars  to  be  measured  may  be  marked  on  the  glass  side  of  the  plate. 
About  8  to  10  stars  are  usually  sufficient  for  a  plate  calibration;  at 
least  6  should  always  be  used. 

C.  Die  Position  of  the  Plate  Center 

In  addition  to  the  coordinates  of  the  selected  stars,  the  coordinates 
of  the  plate  center  must  be  estimated.  Pdut  stars  are  chosen  such  that 
two  lines  intersecting  at  the  center  may  be  drawn  throu^  the  stars  (i.e., 
each  two  stars  will  lie  on  a  line  through  the  center;  the  stars  should  be 
on  opposite  sides  of  the  center).  Die  four  stars  are  identified  on  the  star 
cheurts  and  the  lines  drawn  lightly  on  the  cheurts.  Die  point  of  intersection 
of  the  two  lines  on  the  star  chart  may  be  read  off  as  the  approximate  right 

*A  list  of  star  catalogues  will  be  found  in  Appendix  A. 
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ascension  and  declination  of  the  plate  center. 

D.  The  Measiirement  of  the  Plates 

The  plates  should  be  measured  on  a  measuring  engine  equipped  with 
two  mutually  perpendicvilar  precision  screws  so  that  the  X  and  Y  coordinates 
may  be  measured  simultaneously  in  a  plane  (Fig.  l4).  The  plate  shovild  be 
oriented  in  the  measuring  engine  so  that  the  trail  of  the  re-entry  is 
parallel  to  the  X  axis;  if  the  trail  is  not  completely  straight,  the  X  axis 
should  be  m€ide  parallel  to  the  straight  line  portion.  All  measiarements  of 
the  plate  must  be  made  without  removing  the  plate  from  the  machine.  To  avoid 
possible  backlash  effects,  all  measurements  are  made  with  the  screw  turning 
in  the  same  direction  (customarily  toward  increasing  scale  readings).  To 
check  the  consistency  of  the  readings  over  a  i>eriod  of  time,  two  ixjints  are 
selected  neaa*  the  ends  of  the  trail  and  read  at  occasional  intervals  during 
the  measuring  of  the  plates.  In  this  way  any  disturbances  of  the  machine  by 
temperature  changes  or  mechanical  vibration  will  show  up.  Readings  are 
made  to  the  nearest  micron. 

Various  methods  of  finding  the  coordinate  of  a  point  possible  — 
usually  an  average  of  three  settings  is  sufficient,  as  long  as  the  readings 
do  not  vary  by  more  than  a  few  microns.  It  is  sometimes  advisable  to  repeat 
a  series  of  measurements  with  the  plane  rotated  180*  to  check  for  possible 
systematic  or  accidental  errors.  Whatever  method  of  averaging  is  used,  only 
one  value  of  each  point  is  used  in  the  program.  Therefore  any  averaging 
or  statistical  analysis  on  the  readings  must  be  performed  prior  to  entry 
into  the  program. 

When  the  plate  has  been  adjusted  in  the  machine  so  that  the  X-axis  is 
parallel  to  the  re-entry  trail,  the  center  point  of  each  star  selected  for 
analysis  is  read  to  get  the  coordinates  in  X  and  Y  of  the  star.  Next,  the 
trail,  starting  from  the  first  distinct  point  and  ending  at  the  last 
distinct  point  (whether  or  not  these  points  are  seen  as  dashes)  is  divided 
up  into  from  about  10  to  20  equal  parts  (in  X).  The  Y  measure  ab  each  of 
these  equally  spaced  X's  is  read  for  computing  the  trail  equation. 
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If  the  re-entry  event  is  photographed  with  a  camera  equipped  with  a 
chopping  shutter  the  re-entry  trail  will  appear  as  a  fairly  straight  dashed 
line  --  the  dashes  being  caused  by  the  shutter  breaks  (Pigs.  IJa  and  15b). 

If  the  camera  does  not  have  a  chopping  shutter,  the  re-entry  trail  will 
appear  as  a  streak  among  the  stars  (Fig.  15c).  The  X-measurement  of  the 
middle  of  each  dash  should  be  taken  in  the  case  of  a  chopped  plate.  Ilhe  Y 
coordinate  may  then  be  computed  for  each  dash  by  the  trail  equation  so  that 
Y  need  not  be  read  for  the  dashes.  On  a  streak  photograph  any  point  of 
special  interest,  such  as  a  flare,  diould  be  measTired.  Also,  the  very  first 
and  very  last  points  visible  on  the  trail  must  be  measured.  At  least  five 
measurements  should  always  be  made  along  the  trail. 

A  diagramatic  sketch  of  the  required  measuranents  in  shown  in  Fig.  15* 
Ihis  diagram  is  only  theoretical  and  does  not  represent  a  particular  re¬ 
entry.  Stars  appear  over  the  entire  plate,  but  only  those  near  the  trail 
are  of  interest  here.  Stars  labeled  A-H  have  been  selected  for  measurement. 
These  stars  eure  seen  to  be  selected  on  alternate  sides  of  the  trail  and 
relatively  evenly  spaced  along  the  trail.  Star  F  is  seen  to  lie  on  the  trail 
itself.  Ihe  actual  trail  lies  between  the  two  lines  marked  T.B.  and  T.E. 
(beginning  and  end);  it  will  be  noticed  that  many  dashes  are  missing  at  the 
start  of  the  trail  and  a  few  near  the  end.  The  measurements  for  the  trail 
eqviation  will  be  made  at  eqxial  X  increments  between  TDB  and  TDE,  the  first 
and  last  distinct  points.  In  this  case  the  interval  is  divided  into  10 
equal  increments  (in  X)  labeled  Tl,  T2.  ..TU.  Notice  that  the  reading  is 
taken  along  the  trail  and  that  the  point  read  may  lie  between  two  dashes. 

A  Y  measure  is  taken  for  each  of  the  11  X*s.  The  first  point  T.B.  and  the 
last  point  T.E.  are  also  measured. 

The  center  of  each  dash  should  next  be  resid  in  X  starting  with  the  first 
distinct  dash,  D1  and  continuing  for  the  rest  of  the  distinct  dashes  (one  or 
more  may  be  missing  as  is  dash  l8  in  the  diagram).  The  flare  occurring  at 
dash  17  is  particularly  interesting  since  it  may  also  show  on  the  other  plate 
6uid  so  could  be  used  as  a  common  point.  The  plate  center  should  also  be 
measTired  in  X  and  Y. 
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E.  Canputations 


After  the  measurements  have  been  made,  the  next  step  is  a  plate 
calibration  vhich  will  relate  the  X,  Y  measurements  on  the  plate  to  positions 
in  the  sky.  Ihe  plate  calibration  program  is  run  using  as  input  the  measured 
X,  Y  of  the  known  stars  and  the  coordinates  of  the  stars  from  the  catalogues. 
The  output  consists  of  inverse  plate  constants,  and  direction  cosines  for 
the  plate  center  and  the  coordinate  axes  which  are  needed  as  input  to  the 
optical  trajectory  program.  When  the  results  of  the  plate  calibration  program 
are  available,  the  i*esidxials  (or  A  X  and  A  Y)  for  each  star  are  examined;  if 
any  of  them  are  extremely  large,  the  calibration  is  repeated  omitting  the 
offending  star.  When  the  calibration  is  satisfactory,  i.e.,  the  residuals 
are  very  small,  the  residuals  may  be  plotted  against  X  (Fig.  l6).  If  the 
residuals  are  within  the  limits  of  measurement  they  are  usually  ignored; 
in  other  cases  they  are  applied  to  each  measurement  along  the  trail  as 
discussed  in  Section  II-F. 

When  both  plates  have  been  calibrated  by  the  plate  calibration  program, 
the  input  cards  to  the  optical  trajectory  program  are  made  as  discussed  in 
Section  VI.  The  optical  trajectory  program  uses  the  results  of  the  calibration 
program,  the  geographic  coordinates  of  the  stations,  the  time  of  the  event, 
data  about  the  cameras,  the  measured  X,  Y  readings  taken  at  equal  X  increments 
for  computing  the  trail  equation  and  the  X  readings  for  each  measured  point 
along  the  trail  from  the  two  plates  and  has  as  output  such  meaningful  parameters 
as  range,  height,  direction  cosines,  and  distance  along  the  trail  for  each 
measured  point  along  the  trail. 

F.  The  Reduction  of  Trallblazer  Ik 

To  demonstrate  the  use  of  the  optical  trajectory  program  the  reduction 
of  Trailblazer  Ik  will  be  discussed.  Trallblazer  Ik  was  fired  frtm  Wallops 
Island  on  27  July  1962  at  21:50:00  E.S.T.  and  was  photographed  by  the  Super- 
Schmidt  camerab  at  Arbuckle  Neck  and  at  Eastville  at  21:56:21.583  (Figs.  17a 
and  17b).  The  plates  were  read  at  Heurvard  University  on  a  Mann  Measuring 
Engine  as  discussed  above. 
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Ihe  astronomic  coordinates  and  the  X,  Y  measurements  of  the  stars  are 
given  in  Table  II  for  the  SL  (Arbuckle  Neck)  camera.  Bie  coordinates  have 
alreeidy  been  corrected  for  proper  motion.  Care  must  be  taken  in  using  the  values 
of  proper  motion  given  in  the  star  catalogues,  since  some  catalogues  give 
proper  motion  as  A  and  some  as  A  o:  cos  <5  per  year  or  per  century.  Corrections 
in  proper  motion  are  made  to  the  nearest  month.  Also,  the  X,  Y  measures  have 
been  averaged  to  give  one  reading  for  each  point. 

TABLE  H 

Stars  Identified  and  Measured  for  the  SL  Hiotograph 
of  Trailblazer  IK  taken  at  Arbuckle  Neck 


Name 

Right 

Ascension 

Declination 

Equinox 

X 

Y 

mm 

mm 

A 

20h  6m 

13 -75^8 

-28*  35' 

13.95" 

1950 

173.1015 

20.411 

B 

20h  7m 

28.968s 

-29*  54' 

35.97" 

1950 

178.155 

20.738 

C 

20h  8m 

17.23  s 

1 

0 

0 

52.6  " 

1950 

182.086 

21.1255 

D 

20h  13m 

45.23  s 

-32“  45’ 

45.2  " 

1950 

189.9875 

18.8385 

E 

20h  llftn 

25.78  s 

-33"  50’ 

22.2  " 

1950 

194.171 

19.4495 

F 

20h  17m 

2.57  s 

-34"  44’ 

3^-3  " 

1950 

198.177 

18.4505 

Plate 

Center 

19h  45m 

H 

0 

CO 

H 

1 

1855 

133.205 

23.520 

The  first  step  is  to  perform  a  calibration  of  the  plate.  (We  will  discuss 
only  the  reduction  of  the  SL  data  as  the  SS  data  are  reduced  in  an  identical 
manner. )  Prom  Table  II  we  see  that  we  have  a  total  of  6  stars  all  of  which 
have  their  coordinates  referred  to  the  equinox  of  1950.  We  also  see  that 
the  coordinates  of  the  plate  center  are  referred  to  the  equinox  of  l855«  Bius 
two  sets  of  precession  constants  must  be  used;  the  first  will  precess  the 
coordinates  of  the  plate  center  from  and  the  second  will  precess 

the  coordinates  of  the  stars  from  1950-1963  (the  nearest  1  January).  The 
correct  precession  constants  cure  shown  in  Table  V.  A©  will  be  zero  since 
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the  St5)er-Schmidt  cameras  track  the  stars.  We  vill  run  one  four  constant 
solution,  using  stars  A  and  F  (or  1  and  6).  The  card  set-up  and  a  sample 
computer  print-out  are  given  in  Appendix  B. 

l^n  inspection  of  the  results  of  the  calibration  program,  it  is  seen 
that  the  residuals  are  very  small  and  within  the  accuracy  of  meas virement. 
!naerefore,  no  corrections  need  be  applied  to  the  measured  values  of  X  and 
Y.  Ihe  six  constant  solution  will  be  used  in  the  reduction.  Uie  next 
step  is  to  prepare  the  cards  for  the  Optical  Trajectory  Program. 

The  event  took  place  at  21:56:21.583  E.S.T.  on  27  July  I962.  Converting 
to  Greenwich  time,  we  have  02:56:21.583  G.M.T.  on  28  July  I962.  We  must 
therefore  find  the  0^  for  28  July  in  the  I962  Almanac  since  the  date  at 
Greenwich  must  be  used.  It  is  found  to  be  20h  20m  43.038s.  Uie  Greenwich 
time  and  geographic  coordinates  are  given  on  the  next  card. 

Prom  Table  HI  we  see  that  15  equally  spaced  points  were  read  along 
the  trail.  The  first  X  iioint  is  at  176.000mm  and  an  average  Y  on  the  trail 
is  about  20.l40mm.  These  are  subtracted  from  each  reading  eilong  the  trail 
and  the  differences  ar«  inserted  into  the  program.  This  set  of  cards  will 
be  used  to  compute  the  trail  equation  to  give  Y  as  a  function  of  X. 

The  next  two  cards  will  give  the  inverse  plate  constants  (frcm  Page  4 
of  our  print-out  for  the  six  constant  solution)  and  the  next  three  will 
give  direction  cosines  of  the  plate  center  and  standard  axes  (from  Page  5 
of  the  print-out).  The  names  of  the  variables  in  the  plate  calibration 
print-out  6u:e  identical  to  those  used  in  the  card  write-up  (Section  V). 

The  next  card  gives  the  X  coordinates  of  5  points  of  special  interest 
(the  first  iioint  must  be  the  physically  first  point  read  on  the  trail,  T.B. 
in  our  diagram  and  the  last  point  must  be  the  physically  last  point,  T.E. ). 
The  other  three  points  may  be  any  intermediate  points  such  as  flare  points 
or  other  points  of  special  interest. 

The  next  cards  will  give  information  about  the  dashes  as  diown  in 
Table  IV.  Table  IV  is  not  complete  but  enough  dashes  are  showito  give  an 
understanding  of  the  card  set  up. 
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TABLE  HI 


X,  Y  Readings  at  Equal  X  Increments  Used  in  Computing 
the  Trail  Equation  Data  from  the  SL  Haotograph 
of  Trailblazer  Ik  Taken  at  Arhuckle  Neck 


X 

Y 

AX 

AY 

mm 

mm 

mm 

mm 

176.0 

2O.IUO 

0 

0.000 

177.5 

20.141 

1.5 

0.001 

179.0 

20.148 

3.0 

0.008 

180.5 

20.142 

4.5 

0.002 

182.0 

20.146 

6.0 

0.006 

183.5 

20.142 

7.5 

0.002 

185.0 

20.139 

9.0 

-0.001 

186.5 

20.140 

10.5 

0.000 

188.0 

20.144 

12.0 

0.004 

189.5 

20.142 

13.5 

0.002 

191.0 

20.140 

15.0 

0.000 

192.5 

20.139 

16.5 

-0.001 

19U.0 

20.137 

18.0 

-0.003 

195.5 

20.140 

19.5 

0.000 

197.0 

20.136 

21.0 

-0.004 

X  =  175.861 
End  Point  Y  »  l$B,6Uk 


TABLE  TV 


X  Readixigs  for  Dashes  for  the  SL  Riotograph 
of  Trailblazer  Ik  taken  at  Arhuckle  Neck 


Dash  #  Wt. 

1  4 

2  4 

5  1 

4  2 

5  1 

6  2 


X 

mm 

175.946 

176.245 

176.534 

176.796 

177.057 

177.332 


87 


1 


198.567 
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TABLE  V 


Precession  Constants 


1855-1965  1950-1965 


0.99965576 

0.999991*98 

-0.02412862 

-0.00290555 

-0.010149490 

-0.00126516 

^x 

0.02412862 

0.00290555 

^y 

0.99970885 

0.99999578 

-0.00012664 

-0.00000185 

z 

X 

0.01049490 

0.00126516 

Z 

y 

-0.00012661 

-0.00000185 

z 

0.99994490 

0.99999920 

56 


Table  VI  gives  information  about  the  cameras  and  is  used  in  making 
out  the  remaining  cards.  A  similar  set  of  cards  must  be  made  for  the  SS 
camera*  Appendix  C  gives  a  listing  of  cards  used  as  input  to  the  trajectory 
program,  and  a  sample  computer  output.  Not  all  of  the  card  set  up  and 
print-out  are  shown,  since  Appendix  C  is  meant  to  show  a  sample  calcxalation 
and  not  give  useful  values  of  a  particular  shot. 

TABLE  VI 

Camera  Data  for  Trailblazer  Ik 


SL  Camera 

SS  Camera 

XQ 

133.205 

126.656 

IQ 

25.520 

10.533 

XC 

133.205 

126.656 

YC 

25.500 

10.533 

lj07^ 

i»07U 

Period 

0.1 

0.1 

Occult 

2 

2 

Direction 

Clockwise 

Counter 

IV.  THE  PLATE  CALIBRATION  HIOGRAM 

The  plate  calibration  program  is  designed  to  carry  out  the  computations 
necessary  to  relate  X,  Y  measurements  on  the  photographic  plate  to  astronomical 
coordinates  and  to  compute  residuals  for  the  correction  of  the  measurements 
of  points  along  the  trail.  The  plate  constants  are  computed  by  both  the 
four  and  the  six  constant  solutions  and  the  method  giving  the  best  (i.e. , 
smallest  or  most  consistent)  residuals  is  used  in  the  later  work. 
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The  input  to  the  program  includes  the  right  ascension  and  declination 
and  the  X,  Y  measurements  of  each  star  selected  and  of  the  plate  center. 
®ie  astronomical  coordinates  of  the  stars  are  corrected  for  proper  motion 
to  the  nearest  month  before  being  inserted  into  the  program;  however,  the 
corrections  for  precession  to  the  nearest  1  January  are  performed  within 
the  program  in  direction  cosine  form.  A  generalized  flow  chart  of  the 
plate  calibration  program  and  of  subroutine  SIXCON  appears  in  Jig.  l8. 

The  computer  outptrt  consists  of  the  following  print-outs: 

Page  1  Listing  of  the  inpxrt  data  for  the  plate  center  and  each 
star  in  easily  understandable  form 

Page  2  Direction  cosines  for  the  uncorrected  (A0  =  O)  plate  center 
and  the  coordinate  axes 

Direction  cosines  and  standard  coordinates  {^,  V,  D 
for  each  star 

Page  3* Plate  constants,  inverse  plate  constants,  and  residuals  for 
the  four  constant  solution  (each  four  constant  solution 
will  be  on  a  separate  page) 

Page  It*  Plate  constants  and  inverse  plate  constants  and  residuals 
for  the  six  constant  solution 

Page  ^  Corrected  direction  cosines  for  the  plate  center  and 

coordinate  axes  (Note:  since  in  our  case  A0  =  0,  the 
direction  cosines  on  Page  5  a^re  identical  with  those  on 
Page  2). 

Appendix  B  gives  a  print-oxit  sample  of  the  resvjlts  of  the  plate 
calibration  program. 

The  residuals  should  be  inspected  to  see  if  any  star  gives  either 
very  large  (probably  above  100  microns)  or  very  inconsistent  values.  Such 
a  star  should  be  eliminated  and  the  calibration  program  rerun.  When  the 
residvials  seem  to  be  small  and  lie  on  a  smooth  curve  (to  within  the 
accuracy  of  measxirement),  they  should  be  plotted  against  the  X  reading. 


If  more  than  one  fotir  constant  solution  is  run,  they  will  appear  on 
pages  3' r  etc. 
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Pig.  l6  is  a  sample  curve  of  the  residuals  for  a  theoretical  case.  A  smooth 
curve  is  dravn  through  the  residuals  and  this  curve  (one  for  AX  and  another 
for  ay)  is  used  to  correct  the  sets  of  X,  Y  readings  to  he  used  in  computing 
the  trail  equation  and  the  set  of  X  readings  measured  at  points  to  he  reduced 
(equation  23). 

A  listing  of  the  plate  calibration  program  appears  in  Appendix  D.  The 
six  constant  solution  is  performed  hy  subroutine  SIXCON  which  solves  the 
matrix  equation 

A  X  =  B 


where 


2 


A 


X 


Vi  ^ 


N 

J 


(146) 


using  library  routine  XSIMEQF.  Kiis  routine  solves  equation  lt6  for 
X  by  triangularization  of  the  A  matrix.  If  the  solution  is  successful, 
the  X  matrix  will  replace  the  original  A  matrix  in  storage.  In  the 
program  the  word  MATRIX  is  a  code  which  equals  1  for  a  successful  solution 
of  the  equation,  equals  2  for  overflow  or  underflow,  and  equals  3  for  a 
singular  matrix.  The  arguments  used  in  calling  the  subroutine  XSIMEQF 
together  with  their  values  for  this  example  are; 
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maximum  value  of  subscript  I  of  A  ( I,  j)  =10 
number  of  rows  in  matrix  A  =  5 
number  of  columns  in  matrix  B  =  2 
name  of  matrix  A  =  A 
name  of  matrix  B  =  B 
scale  factor  for  use  in  calculation  =  1 

one  dimensioneil  erasable  array  for  use  in  solution  =  WO  =  10 

A  more  detailed  description  of  the  subroutine  XSIMBQF  is  available. 

The  data  cards  are  set  up  as  follows: 

Data  Cards  for  the  Plate  Calibration  Program 

Card  1  Column  1  contains  a  "1"  and  columns  2-72  contain  a 

Hollerith  title  to  be  written  on  tape  AJ.  This  title  may  contain  any 
identifying  information  desired. 

Card  2  Contains  two  fixed  point  nunbers  (format  215) 

NN  =  number  of  sets  of  precession  constants  used  in  the  reduction 
in  colvmms  1-5  (format  15)  Up  to  5  sets  may  be  used 
N  =  nianber  of  stars  used  in  the  reduction  in  columns  6-10 
(format  15). 

The  next  groiq)  of  cards  will  consist  of  4  cards  for  each  of  the  NN 
sets  of  precession  constants. 

The  first  car^  of  each  precession  set  will  be  a  Hollerith  date  card 
constructed  as  follows: 

A  zero  in  column  1 

Year  of  equinox  of  table  in  columns  2-5 

Hollerith  hyphen  in  column  6 

Year  of  event  to  1  January  in  columns  7-10. 

The  next  3  cards  of  the  precession  set  will  each  have  3  precession 
constants  in  floating  point  form  (format  3E15.8)  as  follows: 

Precession  card  2 
in  columns  1-15 
X  in  columns  16-30 

y 

X  in  columns  31-^5 
z 


Arg^  = 
Argg  = 
Arg^  = 
Arg^  = 
Arg^  = 
Argg  = 
Arg^  = 


1*0 


Precession  card  5 
in  colvmms  1-15 
y  in  colunms  16-50 

y 

Y  in  columns  51-^5 
z 

Precession  card  4 
Z  in  columns  1-15 

X 

Z  in  columns  16-50 

y 

Z  in  columns  51-45 
z 

Each  precession  set  will  also  have  a  number  J  associated  with  it; 

J  being  assigned  according  to  the  placement  of  the  set  in  the  deck,  i.e., 

J  =  1  for  the  first  precession  set  read  in,  J  =  2  for  the  second  precession 
set,  etc.  The  precession  set  to  be  used  in  each  case  will  be  that  one 
corresponding  to  the  equinox  for  which  the  coordinates  of  the  star  (or 
plate  center)  are  recorded  in  the  standard  tables. 

The  next  card  will  have  10  numbers  (foiroats  2P4.0,  J8.4,  2P4.0, 

17.5#  2EL0.5,  15^  EI5.8)  which  give  the  coordin 
as  follows: 

Right  ascension  in  hours 

Right  ascension  in  minutes 

Right  ascension  in  seconds 

Declination  in  degrees 

Declination  in  minutes 

Declination  in  seconds 

X  reading  of  plate  center  in  mm 

y  reading  of  plate  center  in  mm 

J  =  number  of  precession  set  to  use 

A9  =  correction  of  plate  center  for 
stationary  cameras 

(Note:  AO  =  0  for  tracking  cameras.) 

The  next  N  cards  will  give  the  coordinates  of  the  N  stars  to  be  used 
in  the  analysis  and  will  have  9  numbers  per  card  (formats  2P4.0,  p8.4, 
2F4.0,  F7-3;  2P10.5,  I5)«  Corrections  in  ri^t  ascension  and  declination 


es  of  the  plate  center 

in  colxmins  1-4 
in  columns  5-8 
in  columns  ^-16 
in  colxmins  17-20 
in  columns  21-24 
in  colxmins  25-51 
in  columns  52-41 
in  colxmms  42-51 
in  colxmms  52-56 
in  columns  57-71 
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for  proper  motion  must  be  made  before  the  values  are  punched  on  the  cards. 
!I3ie  N  cards  wiH  have  the  following  format; 


Right  ascension  in  hours 

columns  1-4 

Right  ascension  in  minutes 

columns  5-8 

Right  ascension  in  seconds 

colimms  9“16 

Declination  in  degrees 

columns  17-20 

Declination  in  minutes 

columns  21-24 

Declination  in  seconds 

colvmins  25-31 

X  reading  of  plate  in  mm 

columns  32-41 

Y  reeding  of  plate  in  mm 

columns  42-51 

J  =  number  of  precession  set  to  use 

columns  52-5.6 

Each  star  will  have  a  number  L  associated  with  it  depending  on  its 
placement  in  the  deck,  i.e.,L  =  1  for  the  first  star  card  read  in,  L  =  2 
for  the  second  card,  to  L  =  N  for  the  last  star  card  read  in. 

The  next  card  will  have  one  nmber  MM  (foimat  15)  in  columns  l-5> 

MM  is  the  number  of  four  constant  solutions  to  be  run.  (Note;  it  may 
be  desired  to  try  different  combinations  of  stars  as  end  points  in  the 
four  constant  solution;  this  card  allows  for  this  possibility. ) 

The  last  MM  cards  will  each  have  2  numbers  (fonnat  215  )>  the  first 
nimber  is  the  L  for  the  star  used  as  the  beginning  point  in  the  four 
constant  solution  in  columns  1-5  and  the  second  number  is  the  L  for 
the  star  to  be  used  as  end  point  in  columns  6-10. 

V.  THE  OPnCAl  TRAJECTORY  HIOGERAM 

The  optical  trajectory  program  computes  such  physically  meaningful 
parameters  as  range,  hei^t,  direction  cosines,  distance  along  the  trail, 
and  relative  time  for  the  X,  Y  measurements  of  two  photographic  xlates 
of  the  re-entry  event.  Cards  may  be  made  for  a  future  velocity  calciilatlon; 
however,  velocity  is  not  computed  by  the  present  program. 

Input  to  the  program  includes  the  sets  of  X,  Y  readings  taken  at 
equal  X  increments  along  the  trail  to  be  used  to  compute  the  trail  equation, 
the  sets  of  X  readings  taken  at  selected  Intervals  along  the  trail  (such 
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as  the  centers  of  the  dashes)  vhere  a  i^uction  is  desired,  the  time  of  the 
event,  the  geographic  coordinates  of  the  camera  stations,  the  plate  calibration 
data  obtained  from  rtmning  the  calibration  program  for  the  two  plates,  and 
infoimation  about  the  cameras  used  to  photograph  the  event.  Only  one 
position  measurement  is  entered  for  each  point  read;  any  averaging  or 
residual  corrections  necessary  to  obtain  this  measurement  is  done  prior 
to  entry  into  the  program. 

A  generalized  flow  chart  of  the  optical  trajectory  program  and  of 
each  subroutine  appears  in  ilg.  19.  IHie  program  itself  is  given  in 
Appendix  E. 

The  computer  outpvrt  consists  of  the  following  print  oxits: 

1.  A  listing  of  the  input  data 

2.  Time  (Greenwich  and  sidereal)  and  geographic  coordinates  of 
the  two  stations;  direction  cosines  for  vector  from  station  1 
to  station  2  at  sidereal  time  zero  and  at  the  time  of  the  event 
Plate  csilibration  data,  trail  equation,  X,  Y,  eind  direction 
cosines  of  the  five  points  and  direction  cosines  of  the  pole 
for  each  station 

4.  Direction  cosines  of  the  radiant 

5.  Ranges  and  heists  for  the  five  points 

6.  Tables  of  X,  Y,  range,  cos  Z,  and  height  above  mean  sea  level 
for  each  dash  for  the  two  stations 

7.  Tables  of  relative  time,  range,  hei^t,  and  direction  cosines 

in  the  altitude-azimuth  system  for  each  dash  for  the  two  stations 

Note:  Certain  checks  have  been  performed  in  the  program  and  the 
small  errors  found  printed  out  for  inspection.  For  example,  checks  have 
been  made  to  determine  if  the  svm  of  the  squares  of  the  direction  cosines 
of  a  point  equaJsone, 

A  sample  output  is  shown  in  Appendix  C  for  the  Trailblazer  Ik  re-entry. 

A  card  output  containing  values  of  distance  along  the  trail,  relative 
time,  range,  hei^t,  direction  cosines  in  the  altitude-azimuth  system,  dash 
nunber  and  dash  wel^t  is  also  available. 
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The  dash  weight  gives  a  measvire  of  the  confidence  level  that  the 
person  read.ing  the  plate  put  on  that  dash.  !I3ie  weights  are  assigned  as 
follows : 


1  very  poor 

2  poor 

h  fair 

6  good 

8  excellent 

Ihe  weights  are  not  used  in  the  calcvilations  of  this  program  but  are 
included  in  the  card  output  for  possible  use  in  future  velocity  calculations, 
Hie  following  is  a  description  of  the  input  cards  necessary  for 
running  the  optical  trajectory  program. 

Card  1  Column  1  contains  a  ”1"  and  columns  2-72  contain  a  Hollerith 

title  to  be  written  on  tape  A5*  15iis  title  may  be  any  identifying 
information  desired. 

Card  2  will  contain  three  numbers  (formats  2P4.0,  18.4  ) 

ThetaO  in  hours  in  columns  1-4 

llietaO  in  minutes  in  coliznns  5-8 

IbetaO  in  seconds  in  columns  ^-l6 

ThetaO  is  the  sidereal  time  at  Greenwich  for  Oh  Universal  Time  for  the 
day  of  the  event  at  Greenwich  (the  day  will  be  defined  to  start  at  0000 
Universal  Time);  the  thetaO  is  found  in  the  table  of  "Universal  sind 
Sidereal  Times"  of  Pie  American  Ephemeris  and  Nautical  Almanac  for  the 
proper  year  -  thetaO  is  tabulated  separately  for  each  day  of  the  year. 

Ibllowing  cards  1  and  2  which  give  general  information  about  the 
event,  we  will  have  two  sets  of  cards,  one  for  each  station.  In  the 
following  discussion  the  making  of  these  cards  for  station  A  will  be  shown; 
those  for  station  B  are  made  in  exactly  the  same  way  and  immediately  follow 
those  of  station  A  in  the  deck. 

Card  3  will  have  a  zero  in  column  1  and  any  Hollerith  identification 
of  station  A  desired  in  columns  2-72. 

Card  4  will  have  10  numbers  giving  the  coordinates  of  station  A 
(format  5(2P4.0,  18.4),  m.7) 
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Event  time  (Greenwich) 

in  hours 

in  columns  1-4 

Event  time 

in  minutes 

in  columns  5-8 

Event  time 

in  seconds 

in  columns  9-l8 

Latitude  of  Station  A 

in  degrees 

in  columns  17-20 

Latitude 

in  minutes 

in  colxmms  21-24 

Latitude 

in  seconds 

in  columns  25-32 

Longitude  of  Station  A 

in  hours 

in  coltnnns  53-3^ 

Longitude 

in  minutes 

in  columns  37-40 

Longitude 

in  seconds 

in  columns  li-l-I<8 

Elevation  of  Station  A 

in  kilofeet 

in  columns  49-59 

Hie  next  group  of  cards  will  be  used  to  compute  the  trail  equation 
by  the  least  squares  method.  A  starting  and  are  chosen  and  the 
calculation  performed  on  the  differences  from  these  values  where 


A  X  = 
A  Y  = 


X  -  X 

i  o 

Y  -  Y  . 

i  o 


The  X  chosen  is  the  first  X  reading  on  the  trail,  the  Y  is  about 
o  o 

the  average  of  the  Y  readings  along  the  trail. 

Card  5  contains  a  single  number  N  which  is  the  number  of  i>oints  along 

the  trail  (format  15)  in  columns  1-5. 

Card  6  contains  two  numbers  (format  2  E15«8) 

X  in  columns  1-15 

o 

Y^  in  columns  16-50. 

Card  7  etc.  The  next  N/2  or  (N/2)  +  1  cards  contain  the  values  of 
A  Xj.  and  A  Y^  punched  two  to  a  card  (format  4  Ei5*8) 

in  columns  1-15 
in  columns  16-50 


A  1. 

X 


A  X, 
A  Y, 


i+l 

i+1 


in  columns  51-^5 
in  colrmns  t)6-6o 


The  next  five  cards  give  calibration  data  obtained  as  output  from 
the  plate  calibration  program;  each  card  (except  the  last  \rtiich  has  only 
two  numbers)  has  three  numbers  (format  5  E15*8)*  Hiese  values  are  given 
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as  output  fixm  the  plate  calibration  program  vhere  they  are  labeled  by 
the  same  symbol  as  in  this  description. 

Card  1  of  the  5 


AEXI 

in 

columns 

1-15 

BEXI 

in 

columns 

0 

1 

CEXI 

in 

columns 

51-45 

Card  2  of  the  5 

AETA 

in 

columns 

1-15 

BETA 

in 

columns 

16-30 

CETA 

in 

columns 

51-45 

Card  5  of  the  5 

LAMC 

in 

colmns 

1-15 

LANEXI 

in 

columns 

1 

VO 

H 

LAMETA 

in 

columns 

51-45 

Card  4  of  the  5 

MUC 

in 

columns 

1-15 

MUEXI 

in 

colxmms 

16-30 

MUETA 

in 

columns 

V>J 

H 

1 

Card  5  of  the  5 

NUC 

in 

colimms 

1-15 

NUETA 

in 

columns 

1 

H 

The  next  card  gives  5  values  of  Xj  the  first  X  on  this  card  must  be  the 
X  reading  at  the  physically  first  point  of  the  trail  and  the  fifth  X 
on  this  card  must  be  X  reading  at  the  physically  last  point  of  the  trail. 
Die  other  three  X’ s  should  be  ccomon  points  if  any  exist  or  any  other 
particularly  interesting  value  of  X.  Note;  this  card  must  always 
contain  5  X' s  whether  or  not  any  common  points  are  found;  also,  the  first 
and  last  X  on  the  card  must  be  the  physically  first  and  last  X  read  on 
the  event  (format  5FIO.5). 

The  next  card  will  contain  one  number  (format  I5)  which  will  be  the 
number  of  dashes  read  along  the  trail  plus  2  (for  the  end  points  of  the 
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trail  which  in  general  are  not  the  middle  of  dashes) 


NUM  =  number  of  dashes  +2  in  columns  1-5 


Ihe  next  grotqp  of  cards  will  give  the  infonnation  for  each  dash  with  three 
readings  to  a  card.  Ihe  first  reading  will  correspond  to  the  beginning 
of  the  trail  and  the  last  to  the  physically  last  point  on  the  trail 
(these  two  readings  will  not,  in  general,  be  read  from  the  middle  of  the 
dashes).  These  two  end  points  will  have  a  dash  number  of  0  and  a  weight 
of  1  (format  3  (215  ,  HO. 5))* 


Dash  number 
Weight  of  dash 
X  reading  of  dash 
Dash  nunber 
Weight  of  dash 
X  reading  of  dash 
Dash  number 
Weight  of  dash 
X  reading  of  dash 


in  columns  1-5 
in  columns  6-10 
in  columns  11-20 


in  colvmms  21-25 
in  columns  26-30 
in  columns  31-^ 
in  colvnnns  UI-U5 
in  columns  U6-5O 
in  columns  51-60 


The  next  and  last  group  of  the  set  gives  information  about  the 
cameras  used. 

Card  1  of  the  group  contains  fovur  numbers  (format  4  EI5.8) 


Observed  X  coordinate  of 
shutter  center  XQ 


in  columns  1-15 


Observed  Y  coordinate  of 
shutter  center  XQ 


in  columns  I6-3O 
in  columns  31-^5 
in  columns  *46-60 


X  projection  center  XC 
Y  projection  center  YC 


The  next  card  will  contain  3  nmbers  (format  3  EI5.8) 
=  focal  length  squared  in  columns  1-15 


Period  of  revolution  of 
shutter 


in  colvanns  16-30 


Number  of  occultations 
per  revolution 


in  columns  31-^5 


Hie  last  card  for  station  A  will  have  two  numbers  (format  215) 

Bie  first  number  will  be  in  columns  1-5 

=  1  if  card  output  is  desired 
=0  if  cards  are  not  desired 

Die  second  number  will  be  in  coltrains  6-10 

=  1  if  covmterclockwise  revolution  of  St5>er- Schmidt  camera 
=  -1  if  clockwise  revolution  of  Super-Schmidt  camera 
=0  if  non-Super- Schmidt  camera 

Biis  completes  the  input  cairds  for  station  A.  A  second  set  in 
the  same  formats  starting  with  card  5  will  be  required  for  station  B. 
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3-21-6094 

ZENITH  (Z) 


Fig.  la.  The  altitude-azimuth  system. 
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3-21-6095 


Z 


Fig.  lb.  The  declination-hour  angle  and  the  declination-right 
ascension  systems. 
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3-21-6097 
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3-21-6098 
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Fig.  215.  The  declination-hour  angle  system. 
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3-21-6099 


Pig.  2c.  Bie  declination-ri^t  ascension  system. 
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3-21-6100 


Fig.  3a-  Transformation  from  the  ri^t  ascension-declination  system 
to  the  hour  angle-declination  system. 


3-21-6101 


Pig.  3'b.  Transfoimation  from  the  ri^t  ascension-declination  system 
to  the  hour  angle-declination  system. 
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3-21-6102 


Fig.  5c.  Transfonnation  from  the  hour  angle-declination  system  to  the 
altitude-azimuth  system. 


3-21-6103 


Fig.  3cL.  Transfoimation  from  the  hour  angle-declination  system  to  the 
altitude-azimuth  system. 
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3-21-6104 


Pig.  U,  Diagram  of  the  optical  system  used  for  photographing  celestial 
objects. 
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3-21-6105 


z 


P* 


PLg.  5*  Standard  plate  coordinate  system  with  center  at  the  optical 
center  of  the  plate. 
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3-21-6106 


Pig.  6a.  Transformation  ftcoi  the  decllnatlon-rig^t  ascension  system 
to  the  standard  plate  coordinate  system 


Pig.  615.  Transformation  from  the  decllnatlon-rl^t  ascension  system 
to  the  standard  plate  coordinate  system. 
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3-21-6108 


Pig.  7a.  Relative  coordinates  of  station  B  with  respect  to  station  A 
at  sidereal  time  zero  at  station  A. 
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3-21-6109 


Pig.  7b.  Relative  coordinates  of  station  B  with  respect  to  station  A 
at  sidereal  time  of  the  event. 
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3-21-6110 


Fig.  8.  ConstzTictlon  to  show  the  plane  through  the  missile  trail  and 
station  B  and  the  range  vector  to  measured  points  from  station  A. 
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3-21-6111 


Fig.  9»  Construction  to  show  the  method  of  computing  height  above 
mean  sea  level. 
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3-21-6112 


x^g.  10.  The  effect  of  the  rotation  of  the  shutter  on  points  on  the 
photographic  plate. 
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3-312-7062(1) 


Pig.  11,  Camera  locations  in  relation  to  the  launch  site  and  the 
predicted  trajectory  of  the  missile. 
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P^^36-430 


Fig.  12a.  The  Harvard  Super-Schmidt  camera  at  Arbuckle  Neck,  Virginia. 
Bie  small  camera  on  the  ri^t  is  a  baby  Schmidt. 
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Fig.  12b.  Bie  Sviper- Schmidt  optical  system  (f!ran  "Meteor  Science  and 
Engineering,"  Donald  McKinley). 
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Pig.  13a.  Eastville  Super-Schmidt  photograph  of  Trailhlazer  33c  sixth 
stage  re-entry.  Chopping  rate,  I5  chops  per  second. 
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C312-665 


Pig.  13b.  Coquina  Beach  K  2k  aerial  camera  photograph  of  Trailblazer  Ic 
sixth  stage  re-entry.  No  chopping  shutter. 
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C21-265 


Pig.  13c.  Eastvllle  WIID  BC  4  photograph  of  Trallblazer  Ih  sixth  stage 
re-entry.  No  chopping  diutter 
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Jig.  ih  A  2  screw  rotatable  horizontal  measxirlng  engine  of  a  type  , 
suitable  for  reading  ballistic  plates. 
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Pig.  l6a.  Sample  curves  showing  residuals  in  X  measurements. 


3-21-6115 


Pig.  16b.  Sample  cxnrves  showing  residuals  in  Y  measurements 
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-21-6386 
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PLg.  17a.  Arbuckle  Neck  Super-Schmidt  photograph  of  Xrallblazer  Ik  sixth 
stage  re-entry.  Chopping  rate,  20  chops  per  second. 
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Fig.  17b.  Eastvllle  Sxxper-Schmldt  photograph  of  Trallblazer  Ik  sixth 
stage  re-entry.  Chopping  rate,  20  chops  per  second. 
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Pig.  18,  Plow  chart  of  the  plate  calibration  program. 
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Main  Program 


Fig.  19a.  Flow  chart  of  the  optical  trajectory  program. 
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19b.  Plow  chart  of  the  optical  trajectory  program. 
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APPENDIX  A 


The  following  charts  and  catalogues  may  be  used  to  obtain  the 
astronomical,  coordinates  of  the  stars  chosen  for  the  plate  calibration. 


Charts 


Bonner  Durchmusterung  (B  D)  -  23“ 

Cordoba  Resultados  (C  D)  -25*  -  -90“ 

Norton' s  Star  Atlas  very  preliminary  identification  of 
brightest  stars 

Preliminary  Catalogues 

Bonner  Durchmusterung  (B  D) 

Schumfeld' s  Southern  Extension  of  Bonner  Durchmusterung  (BD) 

Cordoba  Durchmusterung  (c  D)  used  mainly  -25“  -  -52“ 

Cape  Photographic  Durchmusterung  (C  P  D)  used  meiinly 
-52*  to  -90“ 

Intermediate  Catalogues 

Henry  Draper  Catalogue  intermediate  catalogues  used  with  the 
Boss  General  Catalogue 

Final  or  Precision  CateJ-ogues 
AGKg  with  EBL  Nord  (for  proper  motion) 

(Zweiter  Katalog  Der  Astronomischen  Gesellschaft) 

Yale  Zone  Catalogues 

(Transactions  of  the  Astronomical  Observatory  of 
Yale  University) 

Cape  Zone  Catalogues  appearing  in  Cape  Annals 
Cape  Zone  Catalogue  with  Cape  Proper  Motion  Catalogue  or 
Cape  Catalogue  of  Paint  Stars 
La  Plata  Catalogues  (or  Cordoba  Catalogues) 

Boss  General  Catalogue 


-90“  -  -  2“ 

+90“  -  +85“ 
+60“  -  +50“ 
+50“  -  -50“ 

_50"  -  -40“ 
-52"  -  -64* 

-40"  -  -52“ 
-64*  -  -82“ 
-90“  -  +90* 


-  2“  -  +90“ 

-  2*  -  -25* 
.23“  -  -90“ 

-40*  -  -90“ 
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ITRAILBLAZER  IK  SL  PLATE 
2  6 
01855-1963 

0.99965377E+00-0.24128620E-01-0.10A94900E-01 
0.24128620E-01+0.99970885E+00-0.12665000E-03 
0.10494900E-01-0. 12661000E-03+0. 99994491 E+00 
01950-1963 

0. 99999498 E+00-0.29055300E-02-0. 1263 1600E-02 
0.29055300E-02+0.99999578E+00-0.ie400000E-05 
0.12631600E-02-0.18400000E-05+0.99999920E+00 

19.  +45. +00. 0000-18. -17. -00. 000  1  PC 

20.  +06. +13. 7540-28. -35. -13. 950+173. 10150+020. 41100 
20.+07.+28.9680-29.-54.-35.970+178. 15500+020. 73800 
20. +08. +17. 2300-30. -56. -52. 600+182. 08600+021. 12550 
20. +13. +45. 2300-32. -45. -45. 200+189. 98750+01 8. 83850 
20. +14. +25. 7800-33. -50. -22. 200+194. 17100+019.44950 
20. +17. +02. 5700-34. -44. -34. 300+198. 17700+018.45050 

1 

1  6 


82 


rsj  (Nj  (M  ru  rg  (M 


Qk 


1TRAILBLA2ER  IK  OPTICAL  TRAJECTORY  SL-SS 
20.+20.+43.0380 
OSL  PLATE  ARBUCKLE  NECIC 

02. +56. +21 .583  +37 .+51 .+23. 2660+05 *+02 .+02 .78 30+00* 0000000 


15 

0.17600000E+03+0.20140000E+02 

O.OOOOOOOOE+00+O.OOOOOOOOE+00+0.15000000E+01+0.10000000E-02 
0.30000000E+01+0.80000000E-02+0.45000000E+01+0.20000000E-02 
0.60000000E+01+0.60000000E-02+0.75000000E+01 +0*200000 OOE-02 

o.9oooooooe+oi-o*ioooooooe-o2+o.io5oooooe+o2+o.ooooooooe+oo 

0*12000000E+02+0*40000000E-02+0.13500000E+02+0.20000000E-02 

0.15000000E+02+O.OOOOOOOOE+00+0.16500000E+02-0.10000000E-02 

0.18000000E+02-0.30000000E-02+0*19500000E+02+O.OOOOOOOOE+00 

0.21000000E+02-0*40000000E-02 

0. 1162 7049E-02-0. 474431 60E-02-0.42743869E-01 

0*47272654E-02-0. 11645432E-02+0.65499187E+00 

0.44365523E+00+0. 8845 1 413 E+00+0. 144238 30E+00 

0.841 17500E+00+0. 46651 330E+00-0.27347735E+00 

0.30918370E+00+0.95100233E+00  . 

175.86100+187*74500+194*74700+197.73400+198*64400 

89 

0  1+175*86100  1  4+175.94600  2  4+176.24500 

3  1+176.53400  4  2+176*79600  5  1+177*05700 

6  2+177.33200  7  4+177.59900  8  4+177.90000 

9  3+178.16400  10  4+178.44300  11  4+178.72800 


TITLE 


COOR  A 
N  TR  A 
XO  YO  A 
TRl  A 
TR2  A 
TR3  A 
TR4  A 
TR5  A 
TR6  A 
TR7  A 
TR8  A 
EXIS  A 
ETAS  A 
LAMB  A 
MU  A 

NU  A 

5  X5  A 

DASMl  A 
0ASH2  A 
DASH3  A 
DASH4  A 


DATA  MISSING 


87  1+198*56700  0  1+198*64400 

0.13320500E+03+0. 23520000E+02+0. 13320500E+03+0* 235000 OOE+02 
0.40740000E+04+0.10000000E+00+0.2000000CE+01 
0-  1 

OSS  PLATE  EASTVILLE 

02. +56. +2 1*583  +37 *+20 *+46. 4300+05. +03 .+36 *7650+00* 0000000 
15 

0.18000000E+03+0.44885000E+02 

O.OOOOOOOOE+OO+O.2OOO0O0OE-O1+0.15OOOOOOE+O1+O.4OOOO0O0E-O2 

0.30000000E+01-0* 10000000E-02+0.45000000E+01+0.10000000E-02 

0.60000000E+01-0*50000000E-02+0*75000000E+01-0*80000000E-02 

0.90000000E+01-0*90000000E-02+0.10500000E+02-0*10000000E-02 

0.12000000E+02-0.50000000E-02+0.13500000E+02-0.13000000E-01 

0.15000000E+02-0.70000000E-02+0.16500000E+02-0.60000000E-02 

0. 18000000E+02-0*30000000E-02+0*19500000E+02+0*50000000E-02 

0.20100000E+02-0*10000000E-02 

0.25734744E-02-0.41302408E-02-0.28662772  +00 

0.41605029E-02-0.25337586E-02+0# 54358602  +00 

0. 7591 0489E+00+0.65095623E+00+0# 39516512  -02 

0. 6509474 lE+00+0. 75911 5 17E+00-0. 33886189  -02 

0.52056018E-02+0*99998645E+00 

176. 84700+186. 71 700+199* 26900+202. 44800+202#52500 
51 

0  1+176.84700  1  1+186*71700  2  1+186.97500 

3  1+187.36800  4  1+187*66000  5  1+187*96200 


DASH30  A 

CAM2  A 

TITLE  B 
COOR  B 
N  TR  B 
XO  YO  B 
TRl  B 
TR2  D 
TR3  B 
TR4  B 
TR5  B 
TR6  B 
TR7  B 
TR8  B 
EXIS  B 
ETAS  B 
LAMB  6 
MU  B 
NU  B 
5  XS  B 

DASHl  B 
DASH2  B 


DATA  MISSING 


52  1  +  202*22000  53  1  +  202*44800  0  l  +  202.«7<(nn 

0.12665600E+03+0*10333000E+02+0.12665600E+03+0*10333000E+02 

0*40740000E+04+0. lOOOOOOOE+OO+O .20000000E+01 
0+  1 


DASH17  B 
CAMl  B 
CAH2  B 
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87 


fRROR«  8.  1.6626451E-09  -3.7252983E-09 
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•  WILSON.FLORENCE  PLATE  CALIBRATION  PROGRAM  OF  FEB.  18.1965 

•  LIST 

•  LABEL 

•  SYMBOL  TABLE 

CPCAL  plate  CALIBRATION  PROGRAM 

DIMENSION  STATEMENTS 

DIMENSION  XXISI.XYISI tXZ(5>.YX(5) .YY(5)«YZ(5) .ZXIS) .ZY(5).ZZ(5). 
1X(2S).Y(Z5).EL(25>.EM(25).EN(25».C0SSIG(25).EXI(25) .ETA(25).TITLE( 
212) 


COMMON  STATEMENTS 
COMMON  TITLE.EXI tCTA.X.Y.N 

FORMAT  STATEMENTS 

20  FORMAT! 12A6) 

21  FORMATdOIS) 

22  FORMAT (72H 

1  ) 

23  FORMAT! 3E15. 8) 

24  F0RMAT!1H  9F13.8) 

25  FORMAT !2F4.0>F8.4.2F4.0.F7.3.2F10.5.15.E15.8) 

26  FORMAT !17H0RT.  ASC.  CENTER-F4.0 .2H  HF4.0.2H  MF8,4.2H  S/13H  DEC.  CE 

1NTER>F4.0.2H  DF4.0.2H  MF7.3.2H  S/4H  XC-F10.5.7H  YC-F10.5/8H  DTH 

2ETA-1PE15.7) 

27  FORMAT!72HO  RIGHT  ASCENSION  DECLINATION  X 

1  Y  /37H  HR  MIN  SEC  DEG  MIN  SEC) 

28  FORMATIIH  F4.0.F5.0.F9.4.F8.0f F5 .0.F8.3.F13.3.F13.3 ) 

29  FORMAT I17H0RT.  ASC.  CENTER-F4.0.2H  HF4.0.2H  MF8.4,2H  S/13H  DEC.  CE 

1NTER-F4.0.2H  DF4.0.2H  MF7,3.2H  S/6H0LAMC>1PE15.7.10H  MUC-lPEl 

25.7.10H  NUC-1PE15.7/8H  LAMETA-1PE15.7.10H  MUETA-1PE15.7. 10 

3H  NUETA-1PE15.7/8H  LAMEX1-1PE15.7.10H  MUEX1-1PE15.7.10H  N 

4UEXI>1PE15.7) 

30  FORMAT (104HONO.  LAMOA  MU  NU 

1  COSSIG  ETA  XI) 

31  FORMATII4.1P6E18.7) 

32  FORMATI3H  A-1PE1S.7.10H  B-1PE15.7.10H  C>1PE15.7.10H 

1  0«1PE15.7/6H  ACAP-1PE15.7.10H  BCAP-1PE15,7,10H  CCAP 

2-1PE15.7.10H  DCAP-1PE15.7/7H  STARS  15. 2H  15) 

33  FORMATI72HONO  X!S)  Y!S)  OX  DY 

1  ) 

34  FORMAT! I3.2F10.3.2F10.5) 

35  F0RMAT!1H  1P3E15.7) 

36  FORMAT !32HOCORRECT ED  COORDINATES  OF  CENTER) 

37  FORMATI6H0LAMC-1PE15.7.10H  MUC-1PE15.7.10H  NUC-1PE15.7/8 

IH  LAMETA-1PE15.7,10H  MUETA-1PE15.7.10H  NUETA-1PE15.7/8H  LAME 
2X1>1PE15.7.10H  MUEXI>1PE15.7»10H  NUEX 1 • 1PE15.7 ) 

38  FORMATI6H0AEXI-1PE15.7.8H  BEX1-1PE15.7 .8H  CEXI •1PE15.7/6H  AETA 

1-1PE15.7.8H  BETA-1PE15.7.8H  CETA-1PE15.7) 

DEFINE  PI  AND  DTRAD-DEGREES  TO  RADIANS 
PI-3.141592654 
OTRAD-3. 141592654/180. 

read  and  write  title 

1  READ  INPUT  TAPE  2. 20. TITLE 
WRITE  OUTPUT  TAPE  3. 20. TITLE 
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C  NUMBER  OF  SETS  OF  PRECESSION  CONSTANTS  AND  NUMBER  OF  STARS 

READ  INPUT  TAPE  2t21fNN.N 

READ  IN  FOUR  CAROS  FOR  EACH  PRECESSION  SET 
CARD  ONE  GIVES  OATES  FOR  PRECESSION 

CAROS  TWO-FOUR  GIVE  PRECESSION  CONSTANTS—THREE  PER  CARO 
00  100  I-lfNN 
READ  INPUT  TAPE  2f22 

READ  INPUT  TAPE  2 <23 tXXI 1 1 tXYl 1 1 tXZ I  1 1 t YX( 1) • YY 1 1 1 « YZ ( It tZXI I ) tZY ( 
llltZZin 

WRITE  OUTPUT  TAPE  3.22 

100  WRITE  OUTPUT  TAPE  3.24 .XXI I ) .XYI It .XZ ( 1 1 .YX( 1 1 . YYI I ) . YZ ( It .ZXI 1 1 .Z 
lYlII.ZZdt 

PLATE  CENTER  CALCULATIONS 
READ  AND  WRITE  PLATE  CENTER  DATA 

READ  INPUT  TAPE  2.2S.ACHR.ACMIN.ACSEC.DCDE6.DCMIN.DCSEC.XC>YC.M,0T 
IH 

WRITE  OUTPUT  TAPE  3.26.ACHR.ACMIN.ACSEC.DCDEG.DCMIN.DCSEC.XC.YC.DT 
IH 

COMPUTE  RIGHT  ASCENSION  AND  DECLINATION  OF  PLATE  CENTER 
IN  DEGREES  AND  RADIANS 
AC-15.* (ACHR+(ACMIN>ACSEC/60. 1/60.) 

0C-DCDEG*IDCMIN+DCSEC/60.)/60. 

AC-0TRA0*AC 
0C-0TRA0*0C 

COMPUTE  UNPRECESSED  DIRECTION  COSINES  OF  PLATE  CENTER 
EQ.  14 
IF(OCt60t61>62 

60  DC— DC 
ENC— SINFCDCI 
GO  TO  63 

61  ENC-0. 

GO  TO  63 

62  ENC-SINFIDC) 

63  ELC-COSFIOC)*COSF(AC) 

EMC-COSF ( DC  I *6 1 NF ( AC ) 

PRECESSION  OF  PLATE  CENTER 

M  IS  THE  PRECESSION  SET  TO  USE  FOR  PLATE  CENTER 

I-M 

PRECESS  PLATE  CENTER  DIRECTION  COSINES  BY  EO.  5 
ELCP-XX ( I ) *ELC+XY  1 1 ) *EMC*XZ  C I ) *ENC 
EMCP- YX ( 1 1 *ELC+YY  < I ) *EMC+ YZ ( 1 1 •ENC 
ENCP-Zxm*ELC*ZYI  I  l*EMC*ZZ<  I  l*ENC 

COMPUTE  DIRECTION  COSINES  OF  XI  AND  ETA  AXES — EO.  15.16 
ENETA-SORTFI l.-ENCP*ENCP ) 

ELETA— ENCP*ELCP/ENETA 
EMETA— ENCP*EMCP/ENETA 
ELEXI— EMCP/ENETA 
EMEXI-ELCP/ENETA 
ENEXI-O. 

WRITE  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3.27 
C 
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CALCULATIONS  FOR  STARS 
PERFORM  LOOP  FOR  EACH  STAR 
DO  101  J-1>N 

READ  IN  DATA  FOR  CURRENT  STAR 

READ  INPUT  TAPE  2»29»AHR*AHINf ASEC.DDEG«DMIN.DSEC»X( Jl .TC Ji >M 
WRITE  OUTPUT  TAPE  3»2S»AHR t AMIN*ASEC tDDEG.DHIN.OSEC.XIJ I t Y IJ) 
COMPUTE  UNPRECESSED  RIGHT  ASCENSION  AND  DECLINATION 
OF  STAR  IN  DEGREES  AND  RADIANS 

65  A-15.»IAHR'f ( AHINtASEC/60*>/60») 

66  D-DDEG-f (DMIN'^OSEC/60.i/60. 

A-DTRAO«A 

D-DTRA0*D 

COMPUTE  UNPRECESSED  DIRECTION  COSINES  OF  STAR  — EO.  4 
IFID)67f6«t69 

67  0— D 

ENS— SINFIOI 
GO  TO  70 

68  ENS>0* 

60  TO  70 

69  ENS-SINF(0l 

70  ELS»C0SF(O»*COSFIA» 

EMS-COSF(0>*SINFIA( 

PRECESSION  OF  CURRENT  STAR 
M'PRECESSION  SET  TO  USE  WITH  CURRENT  STAR 

I-M 

PRECESS  DIRECTION  COSINES  OF  CURRENT  STAR — EO.  5 
ELI  J»»ELS*XX(  n+EMS*xr(  n4ENS*X2(  I ) 

EM<J»-ELS*YXI I »+EMS»YY( I l+ENS*Y2( I » 

EN(J)"ELS«2X( I l+EMS«ZYI I l+ENS»ZZ( I » 

COMPUTE  STANDARD  COORDINATES  OF  CURRENT  STAR— EO.  18.19 
C0SSIGIJI-EL(J»*ELCP4EM(J)*EMCP+EN( J>«ENCP 
EXl(J)-IELIJI*ELEXI«EM(Ji«EMEXI >/COSSIC<J) 

101  ETAI JJ»(ELIJI*ELETA+EM( J»*EMETA+EN( JI«ENETA)/C0SSI6IJ» 

WRITE  OUT  DATA  COMPUTED  FOR  PLATE  CENTER 
WRITE  OUTPUT  TAPE  3.20. TITLE 

WRITE  OUTPUT  TAPE  3 t29.ACHR .ACMIN .ACSEC.OCDEG.DCMIN.DCSEC .ELCP.EMC 
IP.ENCP.ELETA.EMETA.ENETA.ELEXI .EMEXI .ENEXI 
WRITE  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3.30 

WRITE  OUT  COMPUTED  DIRECTION  COSINES  AND  STANDARD 
COORDINATES  FOR  EACH  STAR 
DO  102  J-l.N 

102  WRITE  OUTPUT  TAPE  3.31.J.EL(J) .EMU) .ENIJI .COSSIGIJ) .ETAIJt .EXI <J) 

NUMBER  OF  FOUR  CONSTANT  SOLUTIONS  TO  PERFORM 

read  input  tape  2. 21. mm 

LOOP  FOR  EACH  FOUR  CONSTANT  SOLUTION 
DO  103  J-l.MM 

READ  IN  STARS  TO  USE  FOR  END  POINTS  IN  SOLUTION 
READ  INPUT  TAPE  2.21.L.LL 

PERFORM  FOUR  CONSTANT  SOLUTION 

COMPUTE  A.B.C.D  FOR  THESE  TWO  STARS— EO.  20 
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OEN-<ETA<LL)-ETA(L)  )**2+CEXI  (LD-EXKL)  )**2 
IF<DEN-l.E-10)103»104f lOA 

104  A-CX<LL)-X(L) )»CEXI < LL )-EX I ( L ) ) -C Y C LL )-Y ( L) ) • C ETA { LL )-ETA ( L) ) 
A-A/OEN 

B-<EXHLL)-EXHL)  )  •  « Y<  LL  )-Y  ( L  )  ) -►  <  ETA  ( LL)-ETA  (  L)  )  *  (  X  ( LL)-X  IL  )  ) 
B-B/DEN 

C-X(LL)-A*EXI CLL)-B*ETA(LL) 

0«YCLL)-B*EX1 (LL)+A*ETACLL) 

0EN-A*A^B*B 

C  IF  OEN  TOO  SMALL  IGNORE  THIS  STAR 

IF(OEN-1.E-10)103*105»105 

COMPUTE  INVERSE  PLATE  CONSTANTS— EO.  21 

105  ACAP«A/0EN 
BCAP-B/DEN 
CCAP— ( A*C*»-B»0 ) /OEN 
DCAP-(A#0-B*C)/0EN 
ECAP— ACAP 

WRITE  OUT  PLATE  CONSTANTS  AND  INVERSE  PLATE  CONSTANTS 
FROM  FOUR  CONSTANT  SOLUTION 
WRITE  OUTPUT  TAPE  3*20tTITLE 

WRITE  OUTPUT  TAPE  3 *32 t A *B t C *0* ACAP . BCAP *CCAP *DCAP *L *LL 
WRITE  OUTPUT  TAPE  3 *38  * ACAP *BCAP *CCAP *BCAP tECAP *0CAP 

WRITE  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3*33 

WRITE  INVERSE  PLATE  CONSTANTS  ON  CARDS  FOR  USE  IN 
TRAJECTORY  PROGRAM 
WRITE  OUTPUT  TAPE  14*20*TITLE 

WRITE  OUTPUT  TAPE  14* 35 * ACAP*BCAP *CCAP *BCAP*ECAP *DCAP 

COMPUTE  RESIDUALS  FOR  EACH  STAR  FOR  FOUR  CONSTANT  WORK 
DO  106  I-1*N 

COMPUTE  OX  AND  DY  FOR  CURRENT  STAR— EQ«  23 
0X-A»EXI  C  1  )>B*ETA(  I  )-f*C-X(  I  ) 

0Y«B*EXI C I )-A*ETA( I )+0-Y( I ) 

WRITE  OUT  OX  AND  DY  FOR  CURRENT  STAR 

106  WRITE  OUTPUT  TAPE  3*34* I *X( I ) *Y( I ) *DX*DY 
103  CONTINUE 


PERFORM  A  SIX  CONSTANT  (LEAST  SQUARES)  SOLUTION 
CALL  SIXCON 

COMPUTE  CORRECTED  DIRECTION  COSINES  OF  PLATE  CENTER 
AND  AXES  — EQ.  45 

OTH  READ  IN  AS  SEC  OF  TIME  MUST  BE  CONVERTED  TO  RADIANS 
0TH-0TH*PI/43200* 

SF-SINF(DTH) 

CF-COSF(DTH) 

ELC-ELCP*CF-EMCP»SF 

EMC«EMCP»CF+ELCP*SF 

ELEX-ELEXI»CF-EMEXI*SF 

EMEX»EMEXI*CF+ELEXI»SF 

elet-eleta*cf-emeta»sf 

EMET-EMETA«CF4ELETA#5F 


Sk 


n  n 


WRITE  OUT  CORRECTED  DIRECTION  COSINES 
WRITE  OUTPUT  TAPE  3.20»TITLE 
WRITE  OUTPUT  TAPE  3.36 

WRITE  OUTPUT  TAPE  3.37 .ELC. EMC. ENCP »ELET »EMET .ENETA.ELEX.EMEX.ENEX 
II 

C  WRITE  CORRECTED  DIRECTION  COSINES  ON  CAROS 

121  WRITE  OUTPUT  TAPE  1A.35.ELC.ELEX.ELET.EMC.EMEX.EMET.ENCP.ENETA 
C  DO  NEXT  SET  OF  DATA  IF  ANY 

60  TO  1 
END 
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♦ 

* 


* 


LIST 

LABEL 

SYMBOL  TABLE 
SUBROUTINE  SIXCON 


CSCON 


SUBROUTINE  FOR  SIX  CONSTANT  SOLUTION 


DIMENSION  STATEMENTS 

DIMENSION  TITLE(12) fEXI (25) tETA(25)fX(25)fY(25)fA(10*10)*B(10*l0) * 
ITWO(IO) 

COMMON  TITLE*EXI*ETA*X»Y.N 


FORMAT  STATEMENTS 


20  FORMAT (12A6) 

21  FORMAT(6HOAEXI-lPE15.7feH  BEXI -1PE15#7 .8H  CEXI -1PE15*7/6H  AETA 
l-lPE15#7feH  BETA-lPE15.7f 8H  CETA-1PE15.7) 

22  FORMATC4HOAX-lPE15.7t8H  BX«1PE15.7*8H  CX-1PE15.7/4H  AY-IPE 

115*7.8H  BY-1PE15.7.8H  CY«1PE15.7) 

23  FORMATdH  1P3E15.7/1H  IP3E15.7) 

33  FORMAT C72HONO  X<S)  Y(S)  OX  DY 

1  ) 

34  FORMAT! l3t2F10«3f2Fl0«5) 

35  FORMATC8HOMATRIX-15) 

SET  ALL  SUMS  TO  ZERO  AT  START 

OHE-1, 

SEXI2-0. 

SEXIET-0* 

SEXI-0* 

SEXEXI-0. 

SETA2-0« 

SETA-0* 

SEXETA-0* 

SEX-0* 

SEYEXI-0* 

SEYETA-0* 

SEY-0. 


EN-N 


N-NUMBER  OF  STARS 
PERFORM  LOOP  FOR  EACH  STAR 
DO  100  J-ltN 

SUM  OF  XI  SQUARED 
SEXI2-SEXI2^EXI CJ)*EXI ( J) 

C  SUM  OF  XI  TIMES  ETA 

SEXIET-SEXIET-^EXl  (J)«ETA(J) 

C  SUM  OF  XI 

SEXI-SEXI-fEXKJ) 

C  SUM  OF  X  TIMES  XI 

SEXEXI-SEXEXI^X(J)«EXI (J) 

C  SUM  OF  ETA  SQUARED 

SETA2-SETA2+ETA(J)«ETA(J) 

C  SUM  OF  ETA 

SETA-SETA^ETA(J) 

C  SUM  OF  X  TIMES  ETA 
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c 


SEXETA-SEXETA+X(  J»»ETAU) 

SUM  OF  X 
SEX-SEX-fXU) 

C  SUM  OF  Y  TIMES  XI 

SEYEXl»SEYEXl-*-y(  J)*Exn  J) 

C  SUM  OF  Y  TIMES  ETA 

SE YET A-SEY  ETA+Y I J ) *ET A ( J I 
C  SUM  OF  Y 

100  SEY-SEY+YIJ) 

SET  UP  MATRICES  AS  IN  EQ*  46 
Adtll-SEXia 
A( 1.2)-SEXIET 
A(1>3)-SEXI 
AI2*1)>SEXIET 
A(2«2>-SETA2 
A(2«3)-SETA 
A(3.1I-SEXI 
A13t2»-SETA 
AI3.3)-EN 
Bd.lCSEXEXI 
B(2.U-SEXETA 
B(3«1)-SEX 
Bd*2>-SEYEXI 
Bdt2>-SEYETA 
200  B(3«2»-SEY 

SOLVE  MATRIX  EQUATION  BY  LIBRARY  ROUTINE  XSIMEOF 
HATRIX-XSIMEOFdO«3«2«A.B«ONE«TWO> 

WRITE  OUTPUT  TAPE  3t20tTITLE 

PRINT  OUT  CODE  WORD-1  IF  SOLUTION  SUCCESSFUL 
OTHERWISE  IS  GREATER  THAN  1 
WRITE  OUTPUT  TAPE  3f35. MATRIX 

WRITE  OUT  PLATE  CONSTANTS  FROM  SIX  CONSTANT  SOLUTION 
WRITE  OUTPUT  TAPE  3 •22*Ad  d  I  •  A( 2 1 1 1  .A(  3 1 1  >  •  Ad.2 1  .A( 2 .2  I  .A(  3>2  ) 

COMPUTE  INVERSE  PLATE  CONSTANTS— EQ»  23 
CXY-ACl  tU»AI2.2)-Ad»2  )*A(2f  1» 

ASEXI-A(2«2)/CXY 
BSEXI— A«2tl)/CXY 

CSEXI-(-A(2f2)*A(3.1)*A<2.1>»A{3.2n/CXV 
ASETA—Ad.2>/CXY 
BSETA-AdtD/CXY 

CSETA-(Ad«2>«AI3tl)-Ad«l)«A(3»2»  )/CXY 

WRITE  OUT  INVERSE  PLATE  CONSTANTS 
WRITE  OUTPUT  TAPE  3»21 »ASEXI t BSEXI tCSEXI t ASETA.BSETAtCSETA 
PUT  INVERSE  PLATE  CONSTANTS  ON  CAROS 
WRITE  OUTPUT  TAPE  14,23.ASEXI f BSEXI .CSEXI .ASETA.BSETA.CSETA 

WRITE  OUT  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3.33 

COMPUTE  AND  WRITE  OUT  RESIDUALS  FOR  EACH  STAR 
COMPUTED  BY  SIX  CONSTANT  SOLUTION 
DO  101  J-l.N 


97 


n  n 


OX-A(lf  1)»EXI  IJ)+A(2tl  )«ETA<  JI-fAOtD-Xi  J) 
0Y-A(1*2)«EXI (J)«A<2>2I»ETA( J)4-A(3t2)-r(  J) 
iOl  WRITE  OUTPUT  TAPE  3.3A , J.X C J) , r fJ I .DX.DY 

RETURN  TO  MAIN  PROGRAM 

RETURN 

END 
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WILSONtFLORENCE  OPTICAL  TRAJECTORY  ASSEMBLY  OF  2/17/65 

LIST 

LABEL 

SYMBOL  TABLE 


CTRAJ 


OPTICAL  TRAJECTORY  REDUCTION  PROGRAM 


DIMENSION  STATEMENTS 


DIMENSION  CAXI2)>CAY(2)»CBX(2)»CBY(2).CCXI2ltCCYI2)«COSZ(250>2)tCO 
lSZl(5t2>»CPHI(2)»CPHIPI2)tCTH(2)»D(250>2).DLONG(2).EL(5>2)tELAZ(25 
20>2)»ELC(2)tELETA(2)»ELEV(2)>ELCXI 12 ) tELON ( 2 ) »ELONH( 2 ) t ELONMt 2 ) tEL 
30NSI2)»ELP(2) tELZ(2).EM(5»2)»EMAZ(250t2)»EMC(2)»EMETA(2)>EMEXI(2)> 
4EHP(2) »EMZ(2) tCN(5>2l (ENAZIZSCtZ) •ENC( 2 1 >ENETA< 2 1 . ENP( 2 ) .ENZ( 2 ) tET 
5AP(250t2).EXIP(250t2) *F2(2) .6THRI2I .0TM(2I .GTS(2).H{5.2) *HEIGHT(25 
60»2) .IFEXT(2) »IFPLUSI2)»NDASH(250«2) tNUMI2) »NWT(250<2) tOCCULT(2) ,P 
7l2>»PEL(250t2)tPEM(250>2)fPENI250»2) tPHD(2)>PHI(2).PHIP(2)fPHM(2) t 
8PHS(2)tR(5t2) .RANGE(250f2) .RH0(2) tSINL(2)>SL0PE(2).SPHI (2) tSPHIP(2 
9) tSTH(2)«THETA(2l tTlHE(250t2l t T ITLE< 12  I < TITLEl ( 12»2 ) 

DIMENSION  XI5t2).XBAR(2)«XC(2l •XP(250t2l tXQ(2)»XTR(30) tY(5>2)«YBAR 
H2I *YCt2)tVP«250»2) .YQ(2).yTR(30) tZETAPI 250< 2 ) .ZTRAI L ( 2 ) 


COMMON  STATEMENTS 


COMMON  CAX»CAY»CBX»CBY»CCX.CCY»COL.COSZtCOSZl»CPHI.CPHIP»CTH>D>DL. 
IDLONGtDTRAO.ELtELAZtELCfELETAfELCV.ELEXI tELON •ELONHtELONMtELONSt EL 
2PtELRtELZtEMtEMAZtEMCtEHETAtEMEXI tEMPtEMRtEHZtENtENAZtENCtENETAtEN 
3PtENRtENZtETAABtETAAB0tETAPtEXIABtEXIAB0tEXIPtF2tGTHRtGTMtGTStHtHE 
AIGHTtIFEXTt IFPLUStJtKtNtNDASHtNUHtNWTtOCCULTtPtPELtPEMiPENtPHOtPHI 
StPHIPtPHMtPHStPI tRtRANGEtRHOtSDLtSINLtSINOtSLOPEtSPHI tSPHIPtSTHtTC 
6ORRtTHETAtTHET0HtTHET0MtTHET0StTIMEtTITLEiTITLEI tX .XBAR tXC tXPtXO t X 
TTRtXOtYtYBARtYCtYPfYQtYTRtYOtZETAABtZETAPtZTRAIL 


FORMAT  STATEMENTS 


20  FORMAT(12A6) 

21  FORMAT (101 5) 

22  FORMATI3(2F4.0tF8.4) tFll.7) 

23  FORMAT ( 7H0THET0-F4. OtlHHFA. Of 1HMF8, At IHS) 

24  FORMAT(6H  TIME( I 1 t 2H)>F4tO t 1HHF4.0 t 1HMF8 .4t2HS  /5H  PHI ( 1 1 t2H)>>F4. 0 
ltlHOF4«OtlHMF8.4tlHS/6H  LONG( 1 1 t2H ) -FA tO 1 1HHF4 .0 1 1HMF8 .4 1 1HS/6H  EL 
2EVC llt2H)«F11.7t4H  KFT) 

25  FORMAT(4E15.8) 

26  FORMAT(3E15.8) 

27  F0RMAT(7F10.5) 

28  FORMAT(2I5tF10.5t2I5tF10.5t2I5tF10.5) 

29  FORMAT (38H0ACCUMULAT0R  OVERFLOW  AT  END  OF  DASHPT) 

30  FORMAT (30H0DI VIDE  CHECK  AT  END  OF  DASHPT) 

31  FORMAT (27HONO  CHECKS  AT  END  OF  DASHPT) 

32  FORMAT) 21HOACCOMULA TOR  OVERFLOW) 

33  FORMAT) 13H0DIVIOE  CHECK) 

35  FORMADIH  4E15.8) 

36  FORMATdH  3E15.8) 

37  FORMATdH  7F10.5) 

38  FORMATdH  213tF10.5t2I 5tF10.5t2l5tF10.5) 


TURN  OFF  INDICATORS  AT  START  FOR  LATER  TESTING 


1  IF  ACCUMULATOR  OVERFLOW  lOtlO 
10  IF  DIVIDE  CHECK  lltll 


C 
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DEFINE  PI  AND  DTRAD>DEGREES  TO  RADIANS 
11  DTRAD-1.74532925E-2 
PI-3.1A159265A 

READ  AND  WRITE  TITLE 
READ  INPUT  TAPE  2«20tTITLE 
WRITE  OUTPUT  TAPE  3*20<TITLE 

READ  IN  SIDEREAL  TIME  AT  0  HOURS  2  FOR  DATE  OF  EVENT 
READ  INPUT  TAPE  2t22tTHET0H.THET0H.THET0S 
WRITE  OUTPUT  TAPE  3*23 .THETOH. THETOM.THETOS 

LOOP  FOR  READING  IN  DATA  FROH  EACH  OF  TWO  STATIONS 
DO  100  J-1*? 

j.j 

CURRENT  STATION  TITLE 

READ  INPUT  TAPE  2*20. I TI TLEl i JL* J » t JL*1 • 12 ) 

WRITE  OUTPUT  TAPE  3*20* (TITLE1IJL*J) »JL*1»12) 

DATA  FOR  COMPUTING  RELATIVE  COORDINATES 

GT-GIVEN  TIME  (GREENWICH) - PH-LATITUDE  OF  STATION 

ELON-LONGITUDE  OF  STATION - ELEV-ELEVATION  OF  STATION 

READ  INPUT  TAPE  2*22*GTHR(J) .6TM( J) (GTSI J) .PHO( J) .PHM( J) .PHS( J) (EL 
10NH(J).EL0NM(J)*EL0NS(J)*ELEV( J) 

WRITE  OUTPUT  TAPE  3*24* J.GTHR( J) *GTM( J) .GTS( J) tUtPHO) J) tPHNC J) .PHS 
IIJ) *J*ELONH( J) *ELONH(J) *ELONS< J) *J*ELEV( J) 

DATA  FOR  COMPUTING  TRAIL  EQUATION 

N'NUMBER  OF  POINTS  ALONG  TRAIL — XO*rO>APPROXIMATE  X*Y 

ALONG  TRAIL - XTR*VTR-AOOlTION  TO  X0*y0  TO  GIVE  EXACT 

READING 

READ  INPUT  TAPE  2*21*N 
READ  INPUT  TAPE  2*29*X0*Y0 

READ  INPUT  TAPE  2 *25* ( XTR ( JJ ) * YTR( JJ) * JJ>1 *N) 

WRITE  OUTPUT  TAPE  3*35*X0*Y0 

WRITE  OUTPUT  TAPE  3*35 * (XTR (JJ) *YTR (JJ) *JJ>1 *N ) 

COMPUTE  TRAIL  EQUATION  FOR  STATION  BY  SUBROUTINE  TRAIL 
CALL  TRAIL 

READ  IN  DATA  TO  USE  IN  COMPUTING  POLE  OF  EVENT 
DATA  INCLUDES  INVERSE  PLATE  CONSTANTS  AND  DIRECTION 
COSINES  OF  PLATE  CENTER  AND  COORDINATE  AXES 
DATA  OBTAINED  FROM  PLATE  CALIBRATION  PROGRAM 
*2  READ  INPUT  TAPE  2*26*CAX ( J ) *CBX ( J ) *CCX( J) *CAY( J) *CBY( J) *CCY( J) *ELC 
1(J)*ELEXI( J)*ELETA(J)*EMC(J)*EMEXI( J)*EMETA(J) *ENC( J) *ENETA{ J) 
WRITE  OUTPUT  TAPE  3*36*CAX( J ) *CBX< J ) *CCX ( J) *CAY ( J) *CBY ( J ) *CCY( J ) >E 
1LC(J)*ELEXI(J)*ELETA(J)*EHC( J|*EMEXI(J)*EMETA(J)*ENC(J)*ENETA(J) 

READ  IN  X  VALUES  FOR  FIVE  POINTS  OF  SPECIAL  INTEREST 
READ  INPUT  TAPE  2*27*(X(K*J) *K-1*5) 

WRITE  OUTPUT  TAPE  3*37*(X(X*J)  *)C-1*5) 

READ  IN  DATA  FOR  EACH  DASH 

NUM'NUMBER  OF  DASHES  READ  AT  THIS  STATION 

NDASH-NUMBER  OF  DASH - NWT-WEIGHT(OUALITY)  OF  THIS  DASH 

XP'ACTUAL  X  READING  OF  THIS  DASH 
THREE  DASHES  PER  CARD 
read  input  tape  2*21*NUM(J) 

MUM«NUN(J) 
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READ  INPUT  TAPE  2«28> ( NDASH( Kt J) »NWT (K> J) tXPC J) tK-1 tMUH) 
WRITE  OUTPUT  TAPE  3>38. ( NOASHI J ) iNWT (K >J) tXP (X» J ) tK-1 tMUHl 


READ  IN  DATA  USED  TO  COMPUTE  RELATIVE  TIME 
XO»VQ-ODSERVEO  COORDINATES  OF  CENTER  OF  SHUTTER  ROTATION 
XC.YC-COORDINATES  OF  PROJECTION  CENTER 
F2-F0CAL  LENGTH  OF  CAMERA  SQUARED 

IFPLUS-+1  FOR  COUNTERCLOCKWISE  ROTATION - IFPLUS— 1  FOR 

CLOCKWISE  ROTATION - IFPLUS-0  FOR  NON  SUPER  SCHMIDT 

P-PERIOD  OF  REVOLUTION  OF  SHUTTER - OCCULT-NUMBER  OF 

OCCULTATIONS  PER  REVOLUTION 

IFEXT-+1  IF  CAROS  ARE  TO  BE  MADE - IFEXT-0  IF  NO  CAROS 

ARE  TO  BE  MADE 

301  READ  INPUT  TAPE  2»2S»X0IJ) •YQ( J) tXCI J) (YCI J) •F2( J) tPI J) tOCCULTI J) 
WRITE  OUTPUT  TAPE3.3S.XO< Jl • YQ( J ) tXC( Jl . YC ( J) t F2 I J I tP( J > .OCCULT ( J > 
READ  INPUT  TAPE  2.21. IFEXTI J> .IFPLUSI Jl 
WRITE  OUTPUT  TAPE  3.38, IFPLUSI J) 

READ  IN  DATA  FOR  STATION  2  WHEN  STATION  1  IS  DONE 


100  CONTINUE 


DATA  ALL  IN  COMPUTER - START  CALCULATIONS 

WRITE  TITLE  ON  NEW  PAGE 
WRITE  OUTPUT  TAPE  3, 20, TITLE 

COMPUTE  RELATIVE  COORDINATES  FOR  THE  TWO  STATIONS 
CALL  RELCO 

COMPUTE  THE  POLE  OF  THE  EVENT  FOR  EACH  STATION 
DO  101  J-1.2 
J-J 

CALL  POLE 
101  CONTINUE 

COMPUTE  THE  RADIANT 

CALL  RAO 

PERFORM  LOOP  FOR  EACH  STATION 
DO  102  J-1.2 
J-J 

COMPUTE  RANGES  AND  HEIGHTS  FOR  THE  FIVE  POINTS  OF 
SPECIAL  INTEREST  FOR  CURRENT  STATION 
CALL  RANDH 

COMPUTE  RANGES. HEIGHTS. COSINE  OF  ZENITH  DIRECTION  AND 
DISTANCE  ALONG  THE  TRAIL  FOR  EACH  DASH  OF  CURRENT  STATION 
CALL  OASHPT 

CHECK  INDICATORS  AT  THIS  POINT - IN  CASE  PROJECTION  AND 

SHUTTER  ROTATION  CENTERS  VERY  CLOSE  WILL  SET  DIVIDE  CHECK 

IN  NEXT  SUBROUTINE - THIS  WILL  NOT  HARM  RESULTS 

IF  ACCUMULATOR  OVERFLOW  12.13 

12  WRITE  OUTPUT  TAPE  3.29 
IF  DIVIDE  CHECK  14.15 

14  WRITE  OUTPUT  TAPE  3.30 
GO  TO  15 

13  IF  DIVIDE  CHECK  16.17 

16  WRITE  OUTPUT  TAPE  3.30 
GO  TO  15 

17  WRITE  OUTPUT  TAPE  3,31 
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COMPUTE  RELATIVE  TIMES  FOR  THE  CURRENT  STATION  FOR  EACH  DASH 
DIRECTION  COSINES  IN  ALTITUDE  AZIMUTH  SYSTEM  ALSO 
COMPUTED  FOR  EACH  POINT  IN  SUBROUTINE  RELTIM 
IS  CALL  RELTIM 

COMPUTE  FOR  STATION  2  OR  END  PROGRAM  HERE 

102  CONTINUE 

CHECK  INDICATORS  AT  END  OF  PROGRAM 
IF  ACCUMULATOR  OVERFLOW  200»201 

200  WRITE  OUTPUT  TAPE  3.32 

201  IF  DIVIDE  CHECK  202.203 

202  WRITE  OUTPUT  TAPE  3.33 

203  GO  TO  1 
END 
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*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  TRAIL 

CTRAIL  LEAST  SQUARES  LINE  THROUGH  TRAIL  MAY  5.1964 
C 

C  DIMENSION  STATEMENTS 

DIMENSION  CAX(2)«CAY(2).CBX(2).CBYI2).CCX(2).CCY(2).COSZI250.2t .CO 
lSZl(5.2t.CPHII2I.CPHIP(2).CTHI2I.O(250.2).DLONG(2).EL(5.2I.ELAZ(25 

20.2) .ELC(2).ELETA(2I.ELEV(2)»ELEXI(2).ELON(2I.ELONHI2).ELONM(2).EL 
30NS(2I.ELPI2>.ELZ(2I.EM(5.2).EMAZ<250.2) .EHC(2t .EMETAI2).EMEXII2) . 
4EHP(2>.EMZI2> .ENI 5 .2) .ENAZ ( 250.2 1  .ENCI 2 ) .ENETAI 2 ) .ENP 1 2  I .ENZ 1 2 1 .ET 
5AP(250.2).EXIP(250.2t.F2(2).CTHR(2|.GTM(2) .GTS(2t.Hl5.2).HEI6HT(25 

60.2)  .IFEXT(2) .IFPLUS(2).N0ASH(250.2> .NUM(2(.NWT(250.2) .0CCULT(2) .P 
7(2).PEL<250.2).PEM<250.2).PEN<250.2).PH0(2).PHH2).PHIPC2).PHM(2). 
aPHS(2>.R(5.2) .RAN6E(250.2) .RHO( 2 1 .SI NL(2 ) .SLOPE ( 2 ) .SPHI (2) .SPHIPI2 
9).STH<2I.THETA(2).TIME(250.2).T1TLE(12) .TITLE1(12.2) 

DIMENSION  X(5.2).XBAR(2) .XC(2) .XP(250.2) .XO I  2 ) .XTR ( 30 > . Y( 5 .2 ) . YBAR 
1(2>.YC(2).YP(250.2).YO(2).YTRI30).ZETAP(250.2).ZTRAIL(2) 

C 

C 

C  COMMON  STATEMENTS 

COMMON  CAX.CAY.CBX.CBY.CCX.CCY.CDL.COSZ.COSZl.CPHI .CPHIP.CTH.D.DL. 
IDLONG.DTRAD.EL.ELAZ.ELC.ELETA.ELEV.ELEXI .ELON.ELONH.ELONM.ELONS.EL 
2P.ELR.ELZ.EM.EMAZ.EMC.EMETA.EMEXI.EMP.EMR.EMZ.EN.ENAZ.ENC.ENETA.EN 
3P.ENR.ENZ.ETAAB.ETAAB0.ETAP.EXIAB.EXIAB0.EXIP.F2.6THR.6TM.GTS.H.HE 
4IGHT.IFEXT.IFPLOS.J.X.N.NOASH.NUM.NMT.OCCULT.P.PEL.PEM.PEN.PHD.PHI 
5.PHIP.PHM.PHS.P1.R.RANGE.RHO.SOL.SINL.SINO.SLOPE.SPHI.SPHIP.STH.TC 
60RR. THETA. THETOH.THETOM.THETOS. time. TITLE.T I TLEl.X.XBAR.XC.XP.XO.X 
TTR.XO.Y.YBAR.YC.YP.YO.YTR.YO.ZETAAB.ZETAP.ZTRAIL 
C 

C  ERROR  MESSAGE  FORMAT 

26  FORMATIAIHOOENOMINATOR  LESS  THAN  l.E>10  IN  TRAIL  EO> 

C 

C  START  COMPUTATION 

C  N-NUMBER  OF  POINTS  ON  TRAIL 

EN-N 

C  SET  SUMS  EQUAL  TO  ZERO  AT  START 

SUMX-0. 

SUMY-0. 

SUMXY-0. 

SUMX2-0* 

C 

C  COMPUTE  SUMS  OF  POINTS  ON  TRAIL 

DO  100  JJ-l.N 
C  SUM  OF  X 

SUMX-SUMX+XTRC JJ) 

C  SUM  OF  Y 

SUMY-SUMY+YTR< JJ) 

C  SUM  OF  X  TIMES  Y 

SUMXY-SUMXY+XTR ( JJ )*YTR  <  JJ ) 

C  SOM  OF  X  SQUARED 

SUMX2-SUMX24XTR IJ J ) •XTR  <  J J ) 

100  CONTINUE 
C 

C  COMPUTE  AVERAGE  X.AVERAGE  Y  AND  SLOPE 

XBAR(J)-SUMX/EN 
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YBAR(J)-SUMy/EN 

COMPUTE  SLOPE - EO.  26 

SLOPE ( J ) -SUMXY-EN*XBAR ( J ) «YBAR ( J ) 

EMO-SUMX2-EN»XBAR ( J ) ♦»2 

IF  DENOMINATOR  TOO  SMALL  WRITE  ERROR  MESSAGE 
IF(EM0-l.E-10)92f91f9l 
92  WRITE  OUTPUT  TAPE  3 #26 
91  SLOPE<J)-SLOPE(J)/EMD 

VALUE  OF  AVERAGE  X  IN  MM 
XBARI  J)»XBAR(  J)<^X0 

value  OF  AVERAGE  Y  IN  MM 
102  YBAR(J)»YBAR( J)-fY0 

COMPUTE  Y  INTERCEPT - EG.  25 

ZTRAILf J)-YBAR( J)-SLOPE{ J)*XBAR(J) 

RETURN  TO  MAIN  PROGRAM 

RETURN 

END 
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*  LIST 

*  LABEL 

*  SYMBOL  table 

SUBROUTINE  RELCO 

RC  RELATIVE  COORDINATES 

DIMENSION  STATEMENTS 

DIMENSION  CAX(2) tCAY(2).CBX(2) tCBY(2l .CCX(2) •CCY(2) •COSZI250>2I tCO 
ISZK  5*2I«CPHI (2)»CPHIP(2) «CTH(2) fOIZ^OtZ >  tDLONGI 2 > .EL ( S «2 ) tELAZ ( 25 
20.2>  tELC(2  >*ELETA(2) .ELEVI 2 ) .ELEX I ( 2  I .ELONt 2 ) .ELONH( 2 ) .ELONM ( 2) .EL 
30NS(2).ELP(2I •ELZ(2) .EM ( 5.2 > .EMAZ ( 250.2 ) .EMC! 2 1 .EMETA ( 2 ) .EMEXI ( 2 ) . 
4EHP(2I .EMZ(2I .EN(5.2).ENAZ(250.2).ENC(2I .ENETAI2I.ENP(2I.ENZ(2) .ET 
5AP(250.2).EXIP(250.2) .F2(2) .GTHR(2| .GTHI 2) .GTSI2I .H(5.2|.HEIGHT(25 
60.2).IFEXT(2I .IFPLUS(2>.NDASH(250.2) .NUM(2).NWT(250.2I.OCCULT(2) .P 
7I2).PEL(250.2).PEM(250.2).PEN(250.2).PH0(2>.PHI(2).PHIP(2).PHH(2). 
aPHS(2> .R(5.2) .RANGE(250.2).RH0(2) .SI NL (2 ) .SLOPE ( 2  I .SPHI (2) .SPHIP(2 
9i.STH(2i.THETA(2>.TIHE(250.2) .TITLE(12) .TITLE1(12.2> 

DIMENSION  X(5.2I  .XBAR(2) *XC(2) .XP<250.2) .XO(2) .XTROO) .Y(5.2I .YBAR 
IC2) .YCC2).YP(250.2>.Y0(2).YTR(30) .ZETAP(250.2» .ZTRAIL(2 I 

COMMON  STATEMENTS 

COMMON  CAX.CAY.CBX.CBY.CCX.CCY.COL.COSZ.COSZI.CPHI.CPHIP.CTH.D.DL. 
IDLONG.OTRAD.EL.ELAZ.ELC.ELETA.ELEV.ELEXI .ELON .ELONH.ELONH.ELONS.EL 
2P.ELR.ELZ.EM.EMAZ.EMC.EMETA.EMEXI.EMP.EMR.EMZ.EN.ENAZ.ENC.ENETA.EN 
3P.ENR.ENZ.ETAAB.ETAAB0.ETAP.EXIAB.EXIAB0.EXIP.F2.GTHR.GTH.GTS.H.HE 
4IGHT.IFEXT.IFPLUS.J.X.N.NDASH.NUM.NWT.0CCULT.P.PEL.PEM.PEN.PHD.PHI 
5.PHIP.PHM.PHS.PI.R.RANGE.RHO.SDL.SINL.SINO.SLOPE.SPHI .SPHIP.STH.TC 
6ORR.THETA.THET0H.THET0M.THET0S.TIME.TITLE.TITLE1.X.XBAR.XC.XP.X0.X 
TTR.XO.Y.YBAR.YC.YP.YO.YTR.YO.ZETAAB.ZETAP.ZTRAIL 

FORMAT  STATEMENTS 
20  FORMAT! 12A6) 

24  F0RMATI6H  TIME ( 1 1 .2H) «F4*0. 1HHF4.0. 1HMFB.4 .2HS  /5H  PHK II .2H)«F4.0 
1.1H0F4.0.1HMF8.4.1HS/6H  LONG! I 1.2H)«F4.0.1HHF4.0.1HMFe.4.1HS/6H  EL 
2EV(n.2HI«F11.7.4H  KFTI 

26  FORMATI7H  THETA! I1.2HI-1PE15.7.4H  0EG/5H  ELZ ( 1 1 .2H I -1PC15.7. lOH 

1  EMZ!  I1.2H)-1PE15.7.10H  ENZ ( 1 1 .2H) »1PE15,7/ lOH  PHIPRIMEdl. 

22H)>1PE15*7.4H  OEG/SH  RHO! 1 1 .2HI -IPEIS. 7 .4H  KFT) 

27  FORMATiaH0EXlA80«lPE15.7.14H  KFT  ETAABO-IPE 15«7,4H  XFT/7H  EXIAB- 
11PE15.7.13H  tCFT  ETAAB-1PE15.7.14H  XFT  ZETAAB«1PE15.7.4H  KFT/7H 

2  ERR0R»1PE15.7.4H  KFT) 

COMPUTE  SIDEREAL  TIMES  AND  GEOGRAPHIC  COORDINATES 
FOR  CURRENT  STATION 
DO  101  J-1.2 

JaJ 

WRITE  STATION  TITLE  AND  INPUT  COORDINATES 
WRITE  OUTPUT  TAPE  3 .20. ( T I TLEl < JL. J I . JL«1 . 12> 

WRITE  OUTPUT  TAPE  3 .24 . J.GTHR! J ) .GTM! J) .GTS! J ) . J.PHD! J I .PHMC J) .PHS 
1 ! J) .J.ELONH! J) .ELONM! Jl .ELONS! Jl .J.ELEV! J) 

COMPUTE  UNIVERSAL  TIME  IN  HOURS 
UT*GTHR!J)+!6TM! JI+GTS! Jl/60.)/60. 

C  COMPUTE  CORRECTION  IN  THETA —  EO.  9 

0THETS-236.555*!UT/24. I 
DTHETS-DTHETS/3600. 

C  COMPUTE  THETAO  IN  HOURS 
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TTHET0«THET0H+CTHET0M*THET0S/60.>/60. 

C  COMPUTE  LONGITUDE  OF  STATION  IN  HOURS 

ELONIJ)°ELONH( J>4^IEL0NM( JI+ELONSi J)/60> t/60« 

COMPUTE  SIDEREAL  TIME  IN  DEGREES  AND  RADIANS— EO.  10 
THET*15.»<TTHETO+UT-ELON<J)+OTHETS) 

THETA! J>»THET»DTRAD 
CTH<J»-COSF!THETA<JJ I 
STH! J»»SINF(THETA( JI I 

LATITUDE  CALCULATIONS 

PHI (J»»(PH0I JI+(PHM<J)+PHS<J)/60.I/60.»«0TRAD 
SPHI(J»«SINF(PHI(JJ) 

CPHI  <J)>C0SF(PHI  ( Jn 

COMPUTE  DIRECTION  COSINES  OF  ZENITH  DIRECT  ION— EO.  31 
ELZ(  J)«CPHH  J)»CTH<  JJ 
EM2(JJ»CPHI ( J)»STH(J» 

ENZ(J)«SPHI ( J) 

COMPOTE  GEOCENTRIC  LATITUDE  OF  STATION— EO,  28 

69  PHIPI J J— 695.6635»SINFC2.*PHI ( J) )+1.1731*SINFIA,*PHI I Jl I -.0026*$ IN 
1F(6.*PHHJ)I 

70  PHIP<J)«PHS( JI+PHIPI JI 
PHIP0>PHD( J)*IPHM(J)*PHIP(J>/60. >/60. 

PHIPIJ»«PHIPD*OTRAD 

COMPUTE  EARTHS  RADIOS  OF  STATION— EO.  28 

71  RHO(JI«0. 99832005+1. 683*94E-3*COSF(2. •PHI! JI l-3.549E-6*COSF< A.»PHI 
ICJ) )+8.E-9*COSF<6.*PHl <J> ) 

72  RHOI J)«20.92642796E+3*RHO( J) 

73  RHOI J)»RHO< JI+ELEV( J) 

CPHIP(J)«COSF(PHIP< Jl ) 

SPHIPI JI«SINF(PHIPI JI I 

LONGITUDE  CALCULATIONS 

ELONI JI«15.*<ELONH< JI+IELONM( JI+ELONS< JI/60.I/60.I*OTRAO 

WRITE  OUT  SIDEREAL  TIME .DIRECT  ION  COSINES  OF  ZENITH 
DIRECTION. GEOCENTRIC  LATITUDE  AND  EARTHS  RADIUS 
101  WRITE  OUTPUT  TAPE  3.26. J.THET. J.ELZ I JI .J.EMZI JI .J.ENZI JI .J.PHIPO. J 
l.RHOt JI 


CALCULATE  LONGITUDE  DIFFERENCE  BETWEEN  STATIONS 
LONGITUDE  MEASURED  POSITIVE  TOWARD  WEST 
74  DL«ELON(ll-ELON(2l 

IF  VERY  SMALL  USE  SERIES  TO  CALCULATE  COS  AND  SIN 
I F(ABSF(DL-0.1 1  1107.113.113 

107  CDL-COSFIELONI II l*COSF<ELON<2l l+SlNF(ELON(ll l•SINF«ELON(2I I 
SDL-S INF  I ELON (111 *COSF ( ELON  <  2 1 1 -COSF  <  ELON ( 1 1 1 *S I NF  C  ELON (211 
GO  TO  114 

113  IF(OL}108.109.110 

108  DL—DL 
SDL— SINF(OLI 
GO  TO  111 

109  SDL«0. 

GO  TO  111 

110  SDL-SINF(OLI 

111  CDL-COSFIDLl 


CALCULATE  DIRECTION  COMPONENTS  OF  VECTOR  FROM  STATION  1 
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C  TO  STATION  2  AT  SIDEREAL  TIME  ZERO  AT  STATION  1— EO.  29 

11*  EXIA80-RHO12)*CPHIPt2)«C0L-RHO(l»«CPHIP(l) 

ETAABO«RHO( 2 ) *CPHI P ( 2 » *SDL 
2ETAAB»RH0(2)*SPHIP(2)-RH0(1)*SPHIP(1> 

C  RANGE  BETWEEN  STATION  1  AND  STATION  2--EO.  29 

RAB-SQRTF(EXIAB0»»2+ETAAB0**2+ZETAAB**2) 

COMPUTE  DIRECTION  COMPONENTS  OF  VECTOR  FROM  STATION  1 
TO  STATION  2  AT  TIME  OF  EVENT--EQ.  30 
EXIAB»EXIAB0»CTH( 1 I-ETAAB0»STH( 1 ) 

ETAAB«ETAAB0*CTH(1 )+EXIABO*STH( 1 ) 
ERR0R»S0RTF<EXIA6**2+ETAAB«*2+2ETAAB**2)-RAB 

WRITE  OUT  DIRECTION  COMPONENTS  OF  THESE  VECTORS 
112  WRITE  OUTPUT  TAPE  3>27.EXIAaOtETAABO»EXIABtETAAB>ZETAABtERROR 
C  RETURN  TO  MAIN  PROGRAM 

RETURN 
END 
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LIST 

LABEL 

SYMBOL  TABLE 
SUBROUTINE  POLE 

POLE  CALCULATION  OF  POLE  OF  TRAIL 
DIMENSION  STATEMENTS 

DIMENSION  CAX<2)«CAYI2I •CBX(2)fCBy(2) .CCX(2>tCCr(2)fCOSZ<2SOt2) tCO 
1SZ1<5»2)<CPHI(2)>CPH|P(2) »CTH(2)i0(250<2) fDLON6(2) •EL(5.2>tELAZ<25 
20<2) •ELC(2I>ELETA(2) (ELEVCZ) .ELEXI (2) tELON( 2 ) tELONHI 2 ) tELONMI 2 >  tEL 
30NS<2)»ELP<2) tELZ(2!).EM(S.2) tEMAZ(2S0»2) tEMC(2) .EMETA(21>EMEXI(2) , 
4EMPI2)>EMZ(2I (ENI S >2) t ENAZ ( 250t 2 ) tENCC 2 ) tENETAI 2 ) >ENP< 2 ) tENZI 2 ) tET 
5AP(250t2) tEXIPI2S0t2) tF2(2) tGTHR(2) .GTM(2I iGTSI 2 ) tH( 5t2 ) tHEIGHT425 
60«2)«IFEXTI2) • IFPLUSC 2 > tN0ASH( 250 t 2 ) tNUMI 2 ) •NWT ( 2S0t ^) fOCCULT(2)fP 
7(2)  •PCL(2S0t2) •PEM(250t2) tPEN(250t2) •PH0(2) (PHI (2) •PHIP(2) •PHM(2) t 
8PHS(2)>R(5>2)fRANGE(250»2)>RHO(2)«SINL(2>tSLOPE(2).SPHI(2)fSPHIP<2 
9) tSTHCZ) •THETA(2) •TIME(250>2>  tTITLE(12) •TITLEl{12t2) 

DIMENSION  X(5t2) tXBAR(2) (XCI 2 ) fXPi 250 tZ I >X0 ( 2 ) tXTRC 30) t Y( 5f 2 ) tYBAR 
1(2) >70(2) tYP(250«2) >70(2) .YTR( 30) 'ZETAPI 250>2 ) tZTRAI L (2  ) 

COMMON  STATEMENTS 

COMMON  CAX>CAY«CBX>CBY>CCX>CCY>COL»COSZtCOSZl«CPHI tCPHIPtCTHtOtOL* 
lDLONCtDTRAO<ELtELAZ>ELC>ELETA>ELEVfELEXI tELON.ELONHtELONMtELONStEL 
2PtELRtELZ«EMtEMAZ»EMCtEMETAtEHEXI<EMP.EMR«EMZ<EN<ENAZ«ENC>ENETAtEN 
3P<ENRtENZ>ETAABtETAAB0tETAP.EXIABtEXIAB0>EXIP<F2«GTHR.GTM.GTS«HtHE 
4ICHT. IFEXTtlFPLUStJfX.N.NDASH.NUM.NWT. OCCULT. PiPEL.PEM. PEN. PHO. PHI 
S.PHIPtPHM.PHS.PI.R.RANGE.RHO.SDL.SINL.SINQ.SLOPE.SPHl .SPHIP.STH.TC 
60RR. THETA. THETOH.THETOM.THETOS. time. TITLE.T I TLEI.X.XBAR.XC.XP.XO.X 
TTR.XO.Y.YBAR.YC.YP.YO.YTR.YO.ZETAAB.ZETAP.ZTRAIL 

FORMAT  STATEMENTS 

20  FORMAT (12A6) 

29  FORMAT(5HOCAX-1PE15.7.7H  CBX>1PE15.7.7H  CCX-1PE15.7/5M  CAY-IPE 

115. 7. 7h  CBV-1PE15.7.7H  CCY-1PE15.7/5H  ELC-1PE15.7.9H  ELEXI-1 

2PE15.7.9H  ELETA-1PE15.7/5H  EMC-1PE15.7 .9H  EMEXI-1PE15.7,9H  E 
3META-1PE15.7/5H  ENC-1PE15.7.9H  ENETA-1PE15. 7 ) 

31  F0RMAT(3H0Y-1PE15.7.5H  X  +1PE15.7) 

32  FORMAT (9eH0  X  Y  EL 

lEM  EN  ERROR) 

33  FORMATdH  F10.5.F15.5.1P4E20.7) 

34  F0RMAT(19H  SINL  LESS  THAN  E-5) 

35  F0RHAT(5H0ELP«1PE15.7.10H  EMP-1PE15.7,10H  EMP-1PE15.7/4H 

1  SINL-1PE15.7/14H  TRAIL  LENGTH-1PE15.7.4H  OEG/7HOERROR-1PE15.7.3H 

2  1PE15.7.3H  1PE15.7) 

WRITE  OUT  TITLE.  INVERSE  PLATE  CONSTANTS.  AND  DIRECTION 
COSINES  OF  PLATE  CENTER  AND  COORDINATE  AXES 
WRITE  OUTPUT  TAPE  3. 20. TITLE 

WRITE  OUTPUT  TAPE  3  .20. ( T ITLEl ( JL . J ) . JL>1 .12 ) 

63  WRITE  OUTPUT  TAPE  3.29 .CAX (J) .CBX CJ ) .CCX (J) .CAY (J) .CBY ( J) .CCY( J) .E 
1LCIJ).ELEXI(J).ELETA(J).EMC(J).EMEXI<J).EMETA(J).ENC<J).ENETA(J) 
WRITE  OUT  TRAIL  EQUATION 
WRITE  OUTPUT  TAPE  3 .31 .SLOPE ( J) .ZTRAIL ( J) 

WRITE  HEADS  OF  COLUMNS 
WRITE  OUTPUT  TAPE  3.32 

CALCULATION  OF  STANDARD  COORDINATES  AND  DIRECTION 
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COSINES  FOR  THE  FIVE  POINTS 
PERFORM  LOOP  FOR  THE  FIVE  POINTS 

J-NUHBER  OF  STATION! I. E.  1  OR  21 - K-NUMBER  OF  POINT 

READ  AT  THIS  STATION 
00  201  K>1»9 

COMPUTE  Y - EQ.  25 

Y(K.J|-SLOPE(  Jl^XIKtJI-fZTRAILIJI 

COMPUTE  XIBAR  AND  ETABAR  FOR  GIVEN  POINT—EO.  23 
EXIBAR-CAX! J)*X(XtJI+CBX(J)«Y!K.J)+CCXCJ) 

ETABAR-CAYC J)»XIX.J)+CBYC J>*Y(K*J»+CCY( J» 

COMPUTE  DIRECTION  COSINES  OF  THIS  POINT — EO,  24 
DENOM-SORTF! l.+EXI6AR**2+ETABAR**2| 

EL  I X  t  J )  ■  ( ELC  (Jl  l-ELETA!  J )  »ETABAR4'ELEX  IC  J I  •EX  IBAR I  /DENOM 
EM(K.J»-<EMCCJ)+EMETAIJ»  *ET ABAR*EMEX I ( J I *EX I B AR » /DENOM 
EN  ( K  •  J I  - 1 ENC I J  >  -t-ENETA  I  J  I  •ETABAR )  /DENOM 

SUM  OF  DIRECTION  COSINES  SQUARED  MUST  BE  ONE 
ERRl«ECCX*J)^*2+EMtK»Jl^^2+EN<K»J)^^2-l. 

WRITE  OUT  XtY  AND  DIRECTION  COSINES  OF  THE  POINT 

201  WRITE  OUTPUT  TAPE  3»33 >X(K t Jl • YCK tU > .EL! Xt J ) tEHIK. J I .EN(K t J) tERRl 

COMPUTE  POLE  OF  TRAIL—EQ.  27 
C1»EM!1.J)^EN!5.J)-EN! l.J)^EM!5.J» 

C2-EN!1.J)<EL!5.J)-EL! l.J)^EN!5.J) 

C3-EL(1.JI^EM!5.J)-EM! ltJ)^EL!5.J) 

SINL!J)-S0RTF!C1^^2+C2^^2+C3^^2» 

ASINL-A8SF!SINL!JI I 
IF !ASINL~l.E-5 >202.203.203 

202  WRITE  OUTPUT  TAPE  3.34 

203  IFIASINL-1. >204.204.205 
205  SINL!J>*1. 

204  ELP!J>-C1/SINL! J> 

EMP!U>-C2/SINL!J> 

ENP!J>-C3/SINL! J) 

COMPUTE  ERRORS 

ERRP-ELP ! J  > ••2*EMP ! J  > ••2+ENP ! J  > ••2-1 , 
ERRP1-ELP!J>^EL!1.J>+EMP!J>^EM! 1. J>fENP! J>^EN! 1 .J> 
ERRP2«ELP!J)^EL!5.J>*EMP! J>^EH!5.J>+ENP! J>^EM!5.J> 

COMPUTE  TRAIL  LENGTH 
OLONG!J>>ASINF!SINL!J> I/DTRAD 

WRITE  OUT  DIRECTION  COSINES  OF  POLE. SIN  TRAIL  LENGTH 
AND  LENGTH  OF  TRAIL 

WRITE  OUTPUT  TAPE  3.35 .ELP! J > .EMP! J > .ENP ! J> .SINL ! J> .DLON6! J > .ERRP • 
1ERRP1.ERRP2 
C  RETURN  TO  MAIN  PROGRAM 

RETURN 
END 
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LIST 
LABEL 

SYMBOL  TABLE 
SUBROUTINE  RAO 

RAO  CALCULATION  OF  RAOIANT 

OIMENSION  STATEMENTS 

OIMENSION  CAX(2) fCAYIZ ) t CBXI 2  1 1 CBY ( 2 ) •CCX(2 I tCCY 1 2  I •COSZ(2SO«2 ) tCO 
1SZ1(5«2 ) tCPHI (2) tCPHIPIZ) tCTH(2) tD(2S0t2l tOLONG (2  I t EL (5 .2  I *ELAZI 25 
20tZ>«ELC(2>>ELETA(2l fELEVI2) fELEXl(2l tELONIZ) tELONH(2)«ELONM(2) *EL 
30NS(2I iELP(2) •ELZ(2) •EM(5*2I •EMAZ(250«2l t EMC( 2 ) tEMETAI 2 ) t EMEXI ( 2) t 
AEMP(2l tEMZ(2l •EN(5t2)«ENAZ(250t2l tENCIZI t ENETAI 2 ) .ENP ( 2 t »ENZ( 2 ) tET 
5API250t2l fEXIP(250t2l •F2l2t •GTHR(2I .GTM(2)tGTS(2l tH(5«2l«HEIGHT(25 
60t2l tIFEXT(2l • IFPLUS(2 I tNOASHI 250*2 ) tNUMIZ ) tNWT (250*2 ) tOCCULT ( 2 > tP 
7I2)*PEL(250*2)*PEM(250*2I *PEN(250*2I*PH0(2) *PHI(2) *PHIP(2) *PHM(2| * 
8PHS(2) *R(5*2) *RANGE(250*2I *RH0( 2  I *S INL ( 2 ) *SLOPE( 2 ) «SPHI (21 *SPHIP(2 
91  *STH(2 t  *THETA(2I •TIME(250*2) •TITLE( 12) *TITLE1(12.2) 

OIMENSION  X(5*2I *XBAR(2) *XC(2) *XP(250*2I *XQ ( 2 ) *XTR( 30 ) * Y( 5*2 > • YBAR 
1(2)*YC(2)*YP(250*2) *Y0(2)*YTR(30) *ZETAP ( 250*2 ) .ZTRAIL ( 2 ) 

COMMON  STATEMENTS 

COMMON  CAX*CAY*CBX*CBY*CCX*CCY*CDL*COSZ*COSZl*CPHI *CPHIP*CTH«0*0L* 

iolong*dtrao*el*elaz*elc*eleta*elev*elexi *elon.elonh*elonm.elons*el 

2P*ELR*ELZ*EM*EMAZ*EMC*EMETA*EMEXI*EMP*EMR*EMZ*EN*ENAZ*ENC*ENETA.EN 
3P*ENR*ENZ*ETAAB*ETAABO*ETAP»EXIAB*EXIABO*EXIP*F2*GTHR*GTM*GTS*H*HE 
4IGHT*IFEXT. IFPLUS*J*X*N*NDASH*NUM*NWT*0CCULT*P.PEL*PEM*PEN.PH0*PHI 
5*PHIP,PHM*PHS*PI*R.RANGE*RHO*SOL*SINL*SINO*SLOPE*SPHI *SPHIP.STH»TC 
60RR*THETA*THETOH.THETOM.THETOS*TIME.TITLE.TITLE1*X*XBAR*XC*XP*XQ,X 
7TR*X0*Y.YBAR*YC*YP*YQ*YTR*Y0*ZETAAB*ZETAP*ZTRAIL 

FORMAT  STATEMENTS 
20  FORMAT (12A6) 

37  F0RMAT(8H0RA0IANTI 

38  FORMAT! 5H0ELR«1PE15.7*7H  EMR-1PE15.7*7H  ENR-1PE15 .7/6H  SINO-IP 

1E15.7/7H  ERROR-1PE15.7) 

WRITE  OUT  TITLE 
WRITE  OUTPUT  TAPE  3.20*TITLE 

COMPOTE  RAOIANT - EO.  32 

CA-EMP(1)*ENP(2)-ENP(1 )*EMP(2) 

C5-ENP ( 1 ) tELP ( 2 1 -ELP ( 1 ) *ENP ( 2 ) 

C6-ELP(1 l«EMP(2)-EMP( 1)»ELP(2) 

COMPOTE  SIN  Q - E(3.  32 

SIN0»SQRTF(CA»«2+C5«*2+C6**2) 

ELR-C4/SINQ 
EMR-C5/SINQ 
50  ENR>C6/SIN0 

ERROR-ELR**2+EMR«*2+ENR**2-l.+ELR*ELP(l)+EMR*EMp(l)'fENR*ENP( D+ELR 
l«ELP(  2 )  •••EMR*EMP  ( 2 )  ♦ENR^ENP  ( 2 ) 

IS  RAOIANT  ABOVE  THE  HORIZON 
I F ( ELZ ( II ‘ELR+ENZ ( 1 ) *EMR+ENZ ( 1 1 *ENR ) 7 1 8  *  8 
C  NO*  CHANGE  SIGNS  SO  THAT  IT  IS  ABOVE  HORIZON 

7  ELR—ELR 
EMR— EMR 
ENR—ENR 
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e  WRITE  OUTPUT  TAPE  3t3T 

C  WRITE  OUT  DIRECTION  COSINES  OF  RADIANT  AND  SIN  0 

9  WRITE  OUTPUT  TAPE  3f38tELR>EMRtENR.SINO. ERROR 
C  RETURN  TO  MAIN  PROGRAM 

RETURN 
END 
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SYMBOL  TABLE 
SUBROUTINE  RANOH 

CRANOH  RANGES  AND  HEIGHTS  FOR  THE  FIVE  POINTS 

DIMENSION  STATEMENTS 

DIMENSION  CAXI2)>CAY(2).CBX(2l.C8Y(2).CCX(2)tCCY(2>.COSZ(2SO«2>  tCO 
lSZl(S«2i«CPHII2itCPHIP(2l tCTH(2* tO(2SO>2)>OLONGI2) ,EL(»t2>iELAZI2S 
20>2)»ELC(2I.ELETA(2>>ELEV(2)tELEXI(2)«ELON(2>«ELONH(2)>ELONM(2l tEL 
30NS(2I tELP(2l tELZ(2>>EM(S«2i iEMAZ(250i2t iEMCI2) tEMETAI 2  I tEMEXI ( 2 >  t 
4EMPI2)<EMZ(2I •EN(5«2I »ENAZ ( 250»2 1 tENC (2 1  tENETAI 2 ) tENPl 2 ) tENZ ( 2 ) tET 
5AP(250*2) tEXIP(250»2)iF2(2) tGTHR ( 2 1  tOTMI 2 ) •CTS( 2 ) tHI S tZ ) (HEIGHT ( 25 
60(2) (IFEXT(2I t IFPLUSC 2  I (NDASHI 250(2 ) (NUM(2 ) (NWT ( 250(2 ) (OCCULT! 2 1 (P 
712) (PEL (250(2 ) (PEM(250 (2) (PEN! 250(2 ) (PH0(2) (PHH2) (PHIP(2) (PHMt2) ( 
8PHS(2) (R(5(2I •RANGE(250(2) (RHO( 2 ) (SINLi 2 ) (SLOPE ( 2) (SPHI (2) (SPHIP(2 
9) (STH (2) (THETA (2) (TIME (250(2 ) (TITLE (12) (TITLEK 12(2) 

DIMENSION  X(5(2I (XBAR(2) (XC(2) (XP(250(2l (X0(2) (XTROO) (Y(5(2) (YBAR 
1(2)(YC(2)(YP(250(2I(YQ(2)(YTR(30)(ZETAP(250(2)(ZTRAIL(2) 

COMMON  STATEMENTS 

COMMON  CAX(CAY(CBX(CBY(CCX(CCY(CDL(COSZ(COSZl(CPHI .CPHIP(CTH(D(DL( 
10LONG(DTRAO(EL(ELAZ(ELC(ELETA(ELEV(ELEXI (ELON(ELONH(ELONM(ELONS(EL 
2P(ELR(ELZ(EH(EMAZ(EMC(EMETA(EMEXI tEHP(EMR(EMZ(EN(ENAZ(ENC(ENETA(EN 
3P(ENR(ENZ(ETAAB(ETAABO(ETAP(EXIA8(EXIABO(EXIP(F2(GTHR(GTM(GTS(H(HE 
4IGHT( IFEXT( IFPLUS (J(K( N(NDASH(NUM(NWT (OCCULT (P(PEL (PEM(PEN(PHO(PHI 
S(PHIP(PHM(PHS(PI(R(RANGE(RH0(SDL(SINL(SINQ(SL0PE(SPHI(SPHIP(STH(TC 
60RR(THETA(THETOH(THETOM(THETOS(TIME(TITLE(TITLE1(X(XBAR(XC(XP(XQ(X 
7TR(X0tY(YBAR(YC(YP(Y0(YTR(Y0(ZETAAB(2ETAP(ZTRAIL 

FORMAT  STATEMENTS 
20  FORMAT (12A6) 

33  FORMATdH  F10.5(F15.5(  1P3E20(7( 

41  FORMAT! 77H0  X  Y  R 

IH  DH) 

WRITE  OUT  TITLES  AND  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3(20(TITLE 

WRITE  OUTPUT  TAPE  3(20 ( ( T I TLEl ( JL ( J ) ( JL-1(12 ) 

62  WRITE  OUTPUT  TAPE  3(41 

COMPUTE  DISTANCE  FROM  STATION  TO  PLANE  OF  EVENT 
IS  THIS  STATION  1 
IF( J-l)206(207(206 

NO  COMPUTE  DISTANCE  FROM  STATION  2  TO  PLANE  DETERMINED 
BY  TRAIL  OF  EVENT  AND  STATION  1 - EO(  33 

206  L-1 

SST»-(EXIAfl*ELP(L)4ETAAB«EMP(L)+ZETAAB«ENP(L)) 

GO  TO  2oa 

YES  COMPOTE  DISTANCE  FROM  STATION  1  TO  PLANE  DETERMINED 
BY  TRAIL  OF  EVENT  AND  STATION  2 - EQ.  33 

207  L“2 

SST-EXIAB«ELP(L)'fETAAB*EMP(L)'«-ZETAAB*ENP(L) 

COMPUTE  RANGES(HEIGHTS  AND  COSINE  OF  ZENITH  FOR  EACH 
POINT 

J-NUMBER  OF  STATION — IC-NUMBER  OF  POINT  AT  STATION 
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208  DO  205  K-l»5 

COMPUTE  COS  EPSILON - EO<  34 

SK«EL(KtJ»*ELP(L»+EMIK.J)«EMP(L»+EN<ICtJ»*ENPtLI 
C  COMPUTE  RANGE  OF  POINT  FROM  STATION— EO.  34 

RtlC.J)=SST/SIC 

C  COMPUTE  COSINE  OF  ZENITH  DISTANCE  OF  POINT— EO.  35 

COSZl(KtJI«ELZ(  JI«EL(K»J|4-EM2(  JI*EM(IC.JI-l-ENZ(  J|*EN<K»  Jt 

C  COMPUTE  HEIGHT  ABOVE  SEA  LEVEL  OF  POINT - EO.  38 

Hl»R«K.J»«COSZHK.J» 

210  OH«SORTFtRHOC  JI«RHO(  JI+RttC.J>*R<HtJI'*-2.»RH0I  JJ*H1»-(RH0(J»*H1) 
H(K.J»«H1+DH+ELEV( J» 

WRITE  OUT  RANGE.HEIGHT.  AND  DH  FOR  EACH  POINT 
205  WRITE  OUTPUT  TAPE  3.33 .X (K. J I t Y(K • J> .RIK . J> tH(K., J) .OH 

RETURN  TO  MAIN  PROGRAM 

RETURN 
END 
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LIST 
LABEL 

SYMBOL  TABLE 
subroutine  OASHPT 

COASHPT  RANGES  AND  HEIGHTS  FOR  THE  DASHES 

DIMENSION  STATEMENTS 

DIMENSION  CAX(2) •CAY(2) •CBX(2) •CBY(2) <CCX(2)»CCr<21 •COSZ(2S0t2) tCO 
lSZl(St2) fCPHI (2 }>CPHIP(2)>CTH(2) t0(2S0>2) tOLONGIZ) tELlStZ) tELAZ(2S 
20>2) tELC(2) tELETA(2) tELEV(2) (ELEXI (2 ) *ELON< 2 ) tELONHI 2 ) .ELONMI 2) tEL 
3ONS(2)fELP<2)fELZ(2)>EM<St2)tEMAZ(2S0>2)»EMC(2).EMETA<2)fEMEXI(2). 
AEMP(2)<EMZ(2> •EN(5t2)tENAZ(2$0<2)<ENC(2)>ENETA<2)tENP<2) tENZ(2) tET 
SAP(2S0<2)<EXIP(2S0t2).F2(2)fGTHRI2)fGTM(2)>GTS(2).H(S«2)«HEIGHT<25 
60t2).IFEXT(2)<IFPLUS(2)fN0ASH(2S0«2).NUM<2)tNWT(2S0t2)>OCCULT<2)»P 
7«2  ).PEL«250.2).PEM« 250.2), PEN! 250.2 ».PHD(2JtPHI (2) ,PHIP<2) .PHM(2) . 
8PHS(2 ) .R(5.2) .RANGE<250.2) .RHO(2) .S INL ( 2 ) .SLOPE! 2 ) .SPHI ( 2 ) .SPHIPI 2 
9).STH<2}.THETA(2) .TIME(250.2>.TITLE(12).TITLE1(12.2) 

DIMENSION  X(5.2) .XBAR(2} .XC(2) ,XPi250.2) .XOI 2 ) .XTR! 30) . Y C 5 .2 ) . YBAR 
1(2),YC(2).YPI250.2) .YO(2).YTRI30).ZETAP!250.2) .ZTRAIL(2) 

COMMON  STATEMENTS 

COMMON  CAX.CAY.CBX.CBY.CCX.CCY.CDL.COSZ.COSZl.CPHI.CPHIP.CTH.D.OL, 
IDLONG.DTRAD.EL.ELAZ.ELC.ELETA.ELEV.ELEXI.ELON.ELONH.ELONM.ELONS.EL 
2P.ELR.ELZ.EM.EMAZ.EMC.EMETA.EMEXI .EMP.EMR.EMZ.EN.ENAZ.ENC.ENETA.EN 
3P.ENR.ENZ.ETAAB.ETAA60.ETAP.EXIAB.EXIAB0.EXIP.F2.GTHR.GTM,GTS.H.HE 
AIGHT.IFEXT.IFPLUS.J.K.N.NDASH.NUM.NWT.OCCULT.P.PEL.PEM.PEN.PHO.PHI 
5.PHIP,PHM.PHS.PI.R.RANGE.RHO.SOL.SINL.SINO.SLOPE.SPHI.SPHIP,STH,TC 
eORR.THETA.THETOH.THETOM.THETOS.TIME.TITLE.TITLEl.X.XBAR.XC.XP.XO.X 
TTR.XO.Y.YBAR.YC.YP.YQ.YTR.YO.ZETAAB.ZETAP.ZTRAIL 

FORMAT  STATEMENTS 
20  FORMAT! 12A6) 

24  F0RMAT!1H  1 5 . 2FI5.5 . 1P4E20.7 ) 

25  FORMAT! 106H0  DASH  X  Y  R 

1  COSZ  H  D) 

MUM-NUMBER  OF  POINTS  OR  DASHES  TO  PROCESS 
MUM«NUM!J) 

WRITE  OUT  TITLE  AND  COLUMN  HEADINGS 
WRITE  OUTPUT  TAPE  3.20. TITLE 

WRITE  OUTPUT  TAPE  3 .20. ! T ITLEl ! JL. J ) . JL> 1.12 ) 

62  WRITE  OUTPUT  TAPE  3.25 

IS  THIS  STATION  1 
IF!J-1)102,103.102 

NO  COMPUTE  DISTANCE  FROM  STATION  2  TO  PLANE  DETERMINED 
BY  TRAIL  OF  EVENT  AND  STATION  1— EO.  33 

102  L-1 

SST— !EX1AB*ELP!L)+ETAAB*EMP!L)+ZETAAB*ENP!L) ) 

GO  TO  104 

YES  COMPUTE  DISTANCE  FROM  STATION  1  TO  PLANE  DETERMINED 
BY  TRAIL  OF  EVENT  AND  STATION  2 — EQ.  33 

103  L-2 

SST-EXIAB*ELP!L)+ETAAB»EMP!L)+ZETAAB*ENP!L) 

C  SET  LINE  COUNT  TO  1 

104  M-1 
C 
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COMPUTE  RANGEfHEIGHT. DIRECTION  COSINES  AND  DISTANCES 

along  the  trail  for  each  dash 

J«NUMBER  OF  STATION - K-NUMBER  OF  POINT  AT  STATION 

DO  101  K*1.MUM 

COMPUTE  Y - EQ.  2S 

YP(IC*J ) -SLOPE (J)«XP(IC*J)+ZTRAIL(J) 

COMPUTE  ETA8AR  AND  XIBAR - EG.  23 

EXIBAR«CAX(  J)»XP<XtJ)^CBX(J)*YP(K.JI-»-CCX(J) 

ETABAR-CAYI J)*XP(K*J)+CBY(J)*YP(KtJ)4CCYt J) 

COMPUTE  DIRECTION  COSINES - EG.  24 

DENOM-SQRTFa  .-»-EX  IBAR**2-^ETABAR**2 ) 

PELtX.J)*(ELC(  J)^^ELETA(J)*ETA8AR^ELEXn  J)*EXIBAR)/DEN0M 
PEM(X*J)«(EMC< J )+EMETA( J)*ETABAR+EMEXI ( J ) *EX I  BAR ) /DENOM 
PENCXf J)»(£NCC J)+ENETA( J)*ETABAR»/DENOM 

COMPUTE  COS  EPSILON - EG.  34 

SK»PEL(X»  J)*ELP(L  )+PEMCX*J)^^EMPCL)-^PEN(KtJ)^^ENP(L» 

COMPUTE  RANGE  FROM  STATION - EG,  34 

RANGE(X*J)-SST/SX 

COMPUTE  ZENITH  DIRECTION  COSINE - EG.  35 

COSZ(X.J)«ELZ(JI*PELtX.J)+EMZ( J)«PEM(Kf J)+ENZ(J)*PENCX.J) 
COMPUTE  HEIGHT  ABOVE  SEA  LEVEL - EG.  38 


Hl-RANGEU.  J)*COSZ<IC*J) 

DH-SGRTF CRHO<J)*RHO(J) ♦RANGE CK.J)*RAN6EU»J)>2.*RH0(J)*H1)-CRH0(J) 
l+Hl) 

HEIGHT! KtJ»»Hl+DH+ELEVfJ) 

COMPUTE  DIRECTION  COMPONENTS  OF  POINT - EG.  39 

EXlP(Kf J)«RANGE(X*J)«PEL(X.J) 

ETAP!K.J)-RANGE<X»J)*PEM(X.JI 
ZETAP(ICtJI-RANGE(Xf  J)*PEN(X.J) 

111  0D1«(EXIP(K. JI-EXIP(l.J) )**2^(ETAP(X.J)-ETAP( 1*J) )*«2^(2ETAP(K»J)- 


IZETAPd  .J)  )*»2 
DD1»SGRTF<DD1 ) 


COMPUTE  DISTANCE  OF  POINT  FROM  FIRST  POINT  ON  TRAIL — EG.  40 
IS  THIS  THE  FIRST  POINT 
IF<X-l)l05f 106.105 

YES. DISTANCE  OF  FIRST  POINT 

106  Dl-DDl 

D(X.J).D1 
GO  TO  107 

NO  COMPUTE  DISTANCE  FROM  FIRST  POINT 
105  D(X.J)-DD1-D1 


C 

c 

c 


c 

c 


WRITE  OUT  DASH  NUMBER. X.Y. RANGE. COSINE  ZENITH  DIRECTION. 
HEIGHT.  AND  DISTANCE  ALONG  THE  TRAIL 
107  WRITE  OUTPUT  TAPE  3.24. NDASHI X.J ) .XP ( K .J ) .YP( X.J ). RANGE ( X .J ) .COSZ ( 
lX.J) .HEIGHT! XfJ) .D(X.J) 

INCREMENT  LINE  COUNT 


M*M^1 

IS  PAGE  FINISHED 
IF!M-50)101. 101.108 

YES  START  NEW  PAGE  WITH  PROPER  HEADINGS 
108  WRITE  OUTPUT  TAPE  3. 20. TITLE 

WRITE  OUTPUT  TAPE  3 .20 . ! T I TLE 1 ! JL . J ) . JL- I . 12 ) 
WRITE  OUTPUT  TAPE  3.25 

SET  LINE  COUNT  TO  1  FOR  NEW  PAGE 


M*1 


NO  60  ON  TO  NEXT  POINT 
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101  CONTINUE 


RETURN  TO  MAIN  PROGRAM 

RETURN 

END 
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♦  LIST 

♦  LABEL 

♦  SYMBOL  TABLE 

SUBROUTINE  RELTIM 

CTIME  CALCULATION  OF  TIME  FOR  THE  DASHES 
DIMENSION  STATEMENTS 

DIMENSION  CAX(2) •CAY(2) *CBX(2) •CBY(2) fCCX(2) •CCY(2) fCOSZ(250*2J  fCO 
lSZl(5f2) •CPHI (2) fCPHIP(2»  fCTH( 2 ) .D(250t2 ) .DLONG(2) tEL ( 5 1 2 ) t ELAZ( 25 

20.2) .ELC(2).ELETA(2).ELEV(2).ELEXI(2) •ELON(2l .EL0NH(2).EL0NM(2) .EL 
30NS(2) .ELP(2) .ELZ ( 2 ) .EM( 5 .2 ) . EMAZ ( 250 .2 J .EMC( 2 ) .EMETA (2  )  .EMEXI ( 2 )  . 
4EMP(2)»EMZ(2) tEN(5.2) .ENAZ(250.2) •ENC(2) .ENETAI 2 ) . ENP ( 2 ) .ENZ ( 2 ) .ET 
5APC  250. 2) fEXlP(250.2) fF2(2) . GTHR ( 2 ) . 6TM ( 2 ) .GTS ( 2 ) .H ( 5 . 2 ) .HE  I GHT ( 25 

60.2) . IFEXT(2) . I FPLUSI 2 ) .NDASHI 250 .2 ) .NUM ( 2 ) .NWT ( 250 . 2 ) .OCCULT ( 2 ) .P 
7(2) .PEL (250.2 ) .PEM( 250.2 ) .PEN( 250.2) .PHD( 2) .PHKZ) .PHIP(2) .PHM(2) . 
8PHS(2).R(5.2) .RANGE(250.2) .RHO(2) »S I NL ( 2 ) .SLOPE ( 2 ) .SPHI (2) .SPHIP(2 
9) .STH(2).THETA(2) . T IME ( 250 . 2 ) . T I TLE ( 12 ) . T I TLE 1 ( 1 2 . 2 ) 

DIMENSION  X(5.2 ) .X8AR(2 ) .XC ( 2 ) .XP ( 250 .2 ) .XO ( 2 ) .XTR ( 30 ) .Y(5.2) .YBAR 
1(2) .YC(2) .YPI 250.2 ) .YO(2) .YTR(30) *ZETAP(250.2) .ZTRAIL(2  ) 

COMMON  STATEMENTS 

COMMON  CAX. CAY. CBX.CBY.CCX.CCY. COL. COSZ.COSZl. CPHI .CPHIP.CTH.D.OL . 
IDLONG.DTRAO.EL.ELAZ.ELC.ELETA.ELEV.ELEXI .ELON.ELONH.ELONM.ELONS.EL 
2P.ELR.ELZ.EM.EMAZ.EMC.EMETA.EMEXI .EMP.EMR.EMZ.EN.ENAZ.ENC.ENETA.EN 
3P.ENR.ENZ.ETAAB.ETAABO.ETAP.EXIAB.EXIABO.EXIP.F2.GTHR.GTM.GTS.H.HE 
4IGHT. IFEXT. IFPLUS. J.K.N.NDASH.NUM.NWT.OCCULT.P.PEL.PEM.PEN.PHD.PHI 
5.PHIP.PHM.PHS.PI . R. RANGE. RHO. SDL. SINL.SINO. SLOPE .SPHI .SPHIP.STH.TC 
60RR.THETA.THETOH.THETOM.THETOS.TIME.TITLE.TITLE1.X.XBAR.XC.XP.XO.X 
7TR. XO. Y. YBAR. YC.YP.YO.YTR.YO.ZETAAB.ZETAP.Z TRAIL 

FORMAT  STATEMENTS 

20  FORMAT! 12A6) 

21  FORMAT(4HOXO*1PE15.7.7H  YO- IPE 1 5 • 7 . 7H  XC-1PE15.7.7H  YC*lP 

1E15.7/4H  F2-1PE15.7.7H  P-1PE15.7. IIH  OCCULT-1PE15.7 ) 

22  FORMAT (127H0  DASH  X  TIME  RANGE 

1  HEIGHT  ALPHA  BETA 

2  GAMMA) 

23  FORMATdH  I5.F10.3.1P6E19.7 ) 

24  F0RMAT(22H0DIVIDE  CHECK  AT  OY/DX ) 

26  FORMAT! 15. I3.2F10.6.F9.3.F8.3.3F9#6) 

USE  LONG  FORM  WHEN  PROJECTION  AND  SHUTTER  ROTATION 
CENTERS  MORE  THAN  A  CM  APART 
LONG  FORM  FOR  CALCULATING  OMEGA — EG.  44 
OMEGF(XMEAS)-ATANF(Z*( W»(XMEAS-X01 ) -DYTR ) / ( W*0YTR+ ( XMEAS-XOl ) ) ) 

SHORT  FORM  FOR  COMPUTING  OMEGA - EG.  42 

SOMEGF(XMEAS)-ATANF( ( XMEAS-XGl ) /DYTR ) 

COMPUTE  DT  BY  EG.  43 

DT IF (OMEGA) -0.1591 5494 *SI6N*(OMEGA-OMEGAO)*PER 

CALCULATE  CONSTANTS 
CON-P( J)/OCCULT( J) 

TWOPI-2.*PI 

MUM-NUM(J) 

REDEFINITIONS  OF  QUANTITIES  TO  USE  WITH  FUNCTIONS 
X01-X0( J) 


C 


on  non  on  on  n  non  noon  non 


PER»P< J) 

C  DISTANCE  BETWEEN  CENTER  OF  ROTATION  AND  PROJECTION 

DX  =  XO( J)-XC(  J) 

DY*YO<J)-YC(  J) 

W“DY/DX 

MAY  GIVE  DIVIDE  CHECX  IF  DX  VERY  SMALL — NEED  NOT  AFFECT 
RESULTS  IF  SOMEGF  USED 
IF  DIVIDE  CHECK  6*7 

6  WRITE  OUTPUT  TAPE  3.24 

7  D2*DX**2+DY**2 
Z«S0RTFI1.+D2/F2(J)) 


CALCULATION  OF  OMEGAO 

DISTANCE  OF  POINT  FROM  Y  OF  ROTATION  CENTER 
DYTR*YP( 1. J)-YQ( J) 

X  COORDINATE  OF  THE  FIRST  POINT 
XMEAS-XP(l.J) 

NEED  WE  USE  THE  LONG  FORMULA 
1F(D2-100« ) 10*10*11 

NO  SET  UP  TO  USE  SHORT  FORMULA 

10  ASSIGN  50  TO  KO 

ASSIGN  51  TO  MO 

GO  TO  70 

YES  SET  UP  TO  USE  LONG  FORMULA 

11  ASSIGN  60  TO  KO 

ASSIGN  61  TO  MO 

70  GO  TO  K0*(50*60) 

COMPUTE  OMEGAO  BY  SHORT  FORMULA - EQ.  42 

50  OMEGAO*SOMEGF(XMEAS) 

SET  CODE  FOR  SHORT  FORMULA  USE 

AA1«1 • 

GO  TO  71 

COMPUTE  OMEGAO  BY  LONG  FORMULA - EO*  44 

60  OMEGAO«OMEGF(XMEAS ) 

SET  CODE  FOR  LONG  FORMULA  USE 

AA1*2. 

LOOP  FOR  EACH  POINT 

J-NUMBER  OF  STATION - K*NUMBER  OF  POINT  AT  STATION 

71  DO  100  K=1*MUM 

COMPUTE  DIRECTION  COSINES  OF  POINT  IN  ALTITUDE  AZIMUTH 
SYSTEM - EQ.  12 

ELAZ(K*J)*PEM(K*J)*CTH( J>-PEL(K*J)*STH( J) 

EMAZIK* J)*PEN(K*J)*CPHI ( J I - ( PEL ( K * J ) *CTH ( J ) +PEM ( K * J ) ♦STH ( J ) )«SPHI ( 
IJ) 

ENAZIK* J)«(PEL(K*J)*CTH( JI+PEM(K*J)«STH( J) I «CPH I < J ) +PEN ( K * J ) »SPH I ( 
IJ) 

C  X  MEASUREMENT  OF  POINT 

XMEAS«XP(K. J) 

C  DISTANCE  OF  Y  TO  PROJECTION  CENTER 

DYTR«YP(K*J)-YQ(J) 

C  COMPUTE  OMEGA  OF  POINT  BY  PROPER  FORMULA 

GO  TO  MO*(51*61) 

51  OMEGA*SOMEGF(XMEAS) 

AA2«1. 

GO  TO  72 


rL8 


onnnnn  nnn  onnn  n  n 


61  0MEGA-0ME6F(XMEAS) 

AA2-2# 

C  SIGN  OF  ROTATION  DIRECTION 

72  SIGN*IFPLUS( J) 

C  COMPUTE  OT  FOR  POINT - EQ.  43 

DT»OTIF(OMEGA) 

C  DASH  NUMBER 

75  ENDASH*NDASH(ICf J) 

C  COMPUTE  RELATIVE  TIME  OF  POINT - EO.  41 

100  TIME<KtJ)-CON«ENOASH+OT 

WRITE  OUT  TITLE  AND  CAMERA  PARAMETERS  AND  HEADINGS 
WRITE  OUTPUT  TAPE  3f20fTlTLE 

WRITE  OUTPUT  TAPE  3 t 20 t ( T I TLE 1 ( JL . J ) . JL- 1 . 1 2 ) 

WRITE  OUTPUT  TAPE3 t 21 t XO ( J ) .YO( J) t XC ( J ) t YC ( J ) fF2 C J ) *P ( J ) .OCCULT ( J ) 

103  WRITE  OUTPUT  TAPE  3.22 
SET  LINE  COUNT  TO  1  AT  START  OF  PAGE 

M»1 

WRITE  OUT  DASH  NUMBER. X  READING# RELATI VE  T IME. RANGE. HEIGHT 
AND  DIRECTION  COSINES  ALTITUDE  AZIMUTH  SYSTEM  EACH  POINT 
DO  104  K«1#MUM 

WRITE  OUTPUT  TAPE  3 .23 •NDASH CX . J ) tXP ( X . J ) .T IME ( X . J ) .RANGE ( X . J ) .HE  I 
1GHT(X#J) •ELAZCX#J) .EMAZ(X.J) .ENAZ(XtJ) 

INCREMENT  LINE  COUNT 

M-M^l 

END  OF  THIS  PAGE 
IF(M-50)l04f 104fl06 

YES  SET  UP  NEXT  PAGE 
106  WRITE  OUTPUT  TAPE  3f20#TITLE 

WRITE  OUTPUT  TAPE  3f20# ( TITLEl ( JL#J ) .JL»1 .12) 

WRITE  OUTPUT  TAPE  3.22 

SET  LINE  COUNT  TO  1 

M*1 

NO  60  ON  TO  NEXT  POINT 

104  CONTINUE 

DO  WE  WANT  CARDS 
IF<IFEXT(J))200#201.200 

YES  WRITE  DISTANCE. TIME. range. HEIGHT. DIRECTION  COSINES 
.DASH  NUMBER  AND  DASH  WEIGHT  FOR  EACH  POINT  ON  CARDS 

200  WRITE  OUTPUT  TAPE  14. 20. TITLE 
WRITE  OUTPUT  TAPE  14.20# I T1 TLEl ( JL . J ) . JL«1 . 12 ) 

DO  204  X-1#HUM 

WRITE  OUTPUT  TAPE  14#26#NDASH< X#J) #NWT(X #J) #D(X.J) .T IME(X#J) .RANGE 
1 (X.J) .HEIGHT! X.J) #ELAZ ( X # J ) #EMAZ ! X# J ) #ENAZ ( X# J ) 

204  CONTINUE 
C  RETURN  TO  MAIN  PROGRAM 

201  RETURN 
END 
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